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Batch  LCAP2  is  an  improved  version  of  LCAP,  Ref.  1>  which  was  originally 
developed  in  1966.  The  major  difference  in  usage  between  these  two  programs  is 
the  FORTRAN  callable  implementation  of  each  LCAP2  operator  which  enables  the 
analyst  to  develop  code  to  automate  the  analysis  of  complex  systems. 
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Batch  LCAP2  (Linear  Controls  Analysis  Program)  is  a  FORTRAN  program  which 
provides  the  control  analyst  with  the  capability  to  numerically  perform  clas¬ 
sical  linear  control  analysis  techniques  such  as  transfer  function 
manipulation,  transfer  function  evaluation,  frequency  response,  root  locus, 
inverse  time  response  and  sampled-data  transforms.  It  is  able  to  deal  with  con¬ 
tinuous  and  sampled-data  systems,  including  multiloop  multirate  digital 
systems,  using  s,  z  and  w  transforms. 


2.0  DESCRIPTION 


This  program  was  designed  to  provide  the  control  system  analyst  with  most 
of  the  classical  analysis  tools  needed  for  analyzing  continuous  and  sampled-data 
systems  by  transform  techniques.  A  set  of  transfer  function  and  polynomial  oper¬ 
ators  has  been  defined  in  a  fashion  similar  to  the  instruction  set  of  a 
computer.  Transfer  function  and  polynomial  arrays  ara  defined  to  be  referenced 
with  indices  so  that  they  may  be  easily  addressed  by  the  operators.  The  combina¬ 
tion  of  this  set  of  LCAP2  operators  and  the  form  of  the  data  structure  provides  a 
vary  flexible  and  easy  to  use  program. 

The  data  structure  of  the  program  includes  (1)  s,  z  and  w  plane  transfer 
functions  designated  as  SPTFi,  isl,2...,  ZPTFi,  i=l,2...,  NPTFi  i=l,2..., 
raspectivaly,  and  (2)  polynomials  designated  as  POLYi  ,i=l,2...  Operations  on 
these  transfer  functions  or  polynomials  are  specified  by  references  to  their 
indices.  An  arbitrarily  large  number  of  transfer  functions  and  polynomials  are 
available  to  the  user  since  disk  storage  is  utilized  when  the  number  becomes  too 
large. 

The  transfer  functions  are  represented  as  ratios  of  polynomials.  The  user 
can  load  data  into  a  transfer  function  using  either  the  coefficient  or  the  root 
form  representation.  Data  structures  used  for  the  transfer  functions  and  polyno¬ 
mials  require  that  the  order  of  the  polynomials  be  less  than  fifty. 

Each  of  the  LCAP2  operators  is  implemented  as  a  FORTRAN  subroutine  with  a 
minimal  number  of  arguments  used  to  specify  the  transfer  functions  or  polynomi¬ 
als  involved.  For  example,  to  add  SPTFI  to  SPTF2  and  store  the  results  into 
SPTF3,  the  FORTRAN  statement  is  CALL  SPADD(3,1,2) . 

A  typical  use  of  these  operators  for  a  simple  system  would  be  to  reduce  a 
block  diagram  to  a  single  open  or  closed  loop  transfer  function  using  the  add, 
subtract,  multiply  and  divide  operators.  Then  one  of  the  operators  used  to 
implement  the  classical  control  analysis  techniques  can  be  applied.  For  example, 
if  SPTFi  is  the  open  loop  transfer  function,  the  operator  SFREQ(i)  can  be  used 
to  compute  the  frequency  response  so  that  the  system  can  be  evaluated. 
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For  complex  continuous  systems*  the  shove  block  diegrem  reduction  method 
may  not  be  prscticel  to  apply.  Cramer's  method  for  transfer  function  evaluation 
could  then  be  Used. 

For  sempled-data  systems  z  and  w  plane  operators  are  provided.  Analysis  of 
small  order  systems  can  be  performed  in  either  the  z  or  w  plane.  The  analyst  may 
prefer  to  perform  the  analysis  in  tha  w  plane  since  this  would  allow  the  use  of 
the  Bode  design  techniques.  However*  if  the  order  of  the  system  is  high*  the 
analysis  must  be  performed  in  the  w  plane  since  z  plane  coefficients  cannot  be 
as  accurately  represented  by  tha  computer. 

Multirate  sampled-data  systems  with  integer  rate  sampling  can  be  analyzed 
by  LCAP2.  Two  types  of  operators  based  upon  Sklansky's  frequency  decomposition 
method  are  available  for  this  type  of  analysis. 


3.0  LANGUAGE  AND  HOST  COMPUTER 


The  program  is  written  entirely  in  CDC  (Control  Data  Corp.)  FORTRAN 
EXTENDED  A  with  the  exception  of  one  subroutine  which  is  written  in  assembly 
language.  The  batch  version  of  LCAP2  runs  on  the  CDC  176  under  the  SCOPE  2.1 
Operating  system.  Batch  jobs  typically  require  140-240k  words.  Description  of 
the  code  for  this  program  is  given  in  Ref.  2. 

An  interactive  version  of  LCAP2*  Ref.  3*  is  available  using  the  CDC  INTER¬ 
COM  which  runs  on  the  CDC  835  computer.  INTERCOM  can  also  be  used  to  define  and 
edit  a  batch  file  which  can  be  submitted  to  the  CDC  176. 


4.0  JOB  STRUCTURE 

The  basic  operations  of  a  batch  LCAP2  job  are> 

(1)  Creation  of  FORTRAN  program  -  main  program  and 
subroutines  (optional) 

(2)  Compilation  of  the  source  code  from  (1) 

(3)  Loading  of  routines  from  the  LCAP2  and  system  libraries 

(4)  Execution  of  the  LCAP2  program  from  (2) 


(5)  (Optional)  -  Cataloging  of  data  file  if  one  is 
created  by  LCAP2 


(6)  (Optional)  -  Loading  and  execution  of  the  HARDCPY 
program  to  produce  hardcopy  plots  created  by  LCAP2 

To  facilitate  the  development  of  the  FORTRAN  program  by  the  user,  the  CDC 
UPDATE1  program  is  utilized.  An  LCAP2  program  library  has  been  defined  so  that 
the  first  part  of  the  main  program,  which  contains  many  lines  of  COMMON  block 
and  EQUIVALENCE  statements,  need  not  be  written  by  the  user.  This  block  of  code 
is  copied  from  the  program  library  and  added  to  the  the  user's  FORTRAN  code  to 
create  the  source  code.  The  input  (card  images)  for  the  UPDATE  program  will  be 
of  the  form i 

XI DENT  idname 
XINSERT  START. 1 
XDECK  MAIN 
XCALL  LCAP2 

CALL  INITO  (initialization  of  LCAP2  parameters) 

CALL  MINITO  (initialization  of  matrix  parameters) 


(user's  FORTRAN  code) 


CALL  LEXIT  (required  if  hardcopy  plots  are  generated) 

END 

The  x  in  column  1  defines  an  UPDATE  directive.  The  first  directive,  XIDENT 
idname,  specifies  an  identification  name,  idname,  which  can  be  1  through  9  char¬ 
acters  long.  The  second  directive,  XINSERT  START. 1,  defines  the  location  where 
the  input  data  to  follow  is  to  be  inserted.  The  directive,  XCALL  LCAP2,  will 
write  the  code  in  COMDECK  LCAP2  to  the  file  COMPILE.  This  COMDECK  LCAP2  contains 
tha  main  program  statment  and  all  of  the  COMMON  block  and  EQUIVALENCE  statements 
required  by  the  main  program.  The  remaining  input  data  are  the  user's  FORTRAN 
code  which  will  be  copied  to  the  file  COMPILE  to  complete  the  creation  of  the 
main  program.  The  words  in  parenthesis  are  comments  and  are  not  part  of  the  FOR¬ 
TRAN  code. 

The  job  control  cards  for  setting  up  the  above  operations  are  given  in  the 
next  section. 


1  Tha  UPDATE  program  maintains  and  updates  source  decks  for  libraries  under 
the  SCOPE  2,  NOS  1,  and  NOS/BE  1  operating  systems. 


Two  forms  for  job  submittal  are  given  below.  The  first  will  be  an  explicit 
one  which  includes  a  complete  list  of  control  cards  required.  The  second  is  a 
shortened  form  which  attaches  and  uses  a  procedure  to  generate  the  control 
cards . 

The  first  form  is: 


(  control  cards  for  accounting  ) 
FILE,TAPE30,BT=I. 

ATTACH(TAPE30, lfn, 10= . , ST=PF6) 

ATTACH(0LDPL,8LCAP2PLX, 10=9487) 

UPDATE. 

FTN(I=C0MPILE,R=3) 

FILE,TAPE31,BT=I. 

REQUEST (TAPE31 ,  XPF) 

RETURN(OLDPL) 

ATTACH(LCAPLIB,8LCAP2LIBX, 10=9487) 
ATTACH(PLOTLIB, 3FTNPL0TLIB) 
LIBRARY(LCAPLIB, PLOTLIB) 

L  DSET ( PRESET =ZER0 ) 

LGO. 

CATAL0GCTAPE31 , 8f ilename, 10= . . ST=PF6 ) 

HARDCPY, ST=IBMD8 . 

XEOR 


(optional,  use  only  if  old 
data  is  to  be  restored) 
(attach  LCAP2  program  library) 

(compile  output  of  UPDATE) 
(optional,  only  if  LCAP2 
STORE  operator  is  to  be  used) 

(attach  LCAP2  library) 

(attach  plot  library) 


(load  and  execute  LCAP2  program) 
(optional,  use  only  if  LCAP2 
STORE  operator  was  used) 

(omit  argument  if  A3  plotter  desired) 
(end  of  record) 


(UPDATE  input  deck  as  described 
in  previous  section) 


XEOR 


(end  of  record) 


( 


control  cards  for  accounting  ) 


FIIE,TAPE38»BY*I . 

ATTACH(TAPE30, Ifn, ID* . ,ST*PF6) 

ATTACH(X#  8LCAP2CC# ID=9487 ) 

BEGIN, LCAP2CC » X . 

CATAL0G( TAPE31 , 8f i lename, ID* . . ST=PF6 ) 

HARDCPY,ST*IBMD8. 

NEOR 


(optional#  usa  only  if  old 
data  is  to  ba  restored) 

(attach  LCAP2  control  card  procedure) 

(optional#  usa  only  if  LCAP2 
STORE  operator  was  used) 

(omit  argument  if  A3  plotter  desired) 
(end  of  record) 


(UPDATE  input  deck  as  described 
in  previous  section) 


NEOR 


(end  of  record) 


In  the  second  form,  the  file  X  will  generate  the  same  control  cards  as  the 
first  form  except  for  the  (1)  FILE,TAPE30, . . . ,  (2)  ATTACH(TAPE30, .  .  . ,  (3)  CATA- 
L0GCTAPE31, . . .  ,  and  (4)  HARDCPY  statements.  It  is  recommended  that  the  second 
form  ba  used  unless  the  user  must  change  some  of  the  control  cards.  An  example 
when  this  is  necessary  is  if  the  print  limit  is  exceeded.  The  statement  LGO. 
should  be  changed  to  LG0(PL* . )  where  the  value  of  PL  is  the  number  of  print 


The  following  is  a  list  of  the  LCAP2  operators  grouped  by  type  of 
operation.  A  brief  description  of  each  operation  is  included.  A  more  detailed 
description  of  each  operator  is  given  in  Appendix  A. 


PADD(i, i.k) 


PEQUCi, i) 


PLDC(i) 


PLDR(i) 


PMPYCi, i.k) 


PPRNCi  ) 


PRTSCi) 


PSUB(i, i.k) 


-  Polynomial  Add 
POLYi  =  POLYi  +  POLYk 

-  Polynomial  Equal 
POLYi  =  POLYi 

-  Polynomial  Load.  Coefficient  Form 
POLYi  =  POLY 

-  Polynomial  Load.  Root  Form 

ROOTi  =  ROOT,  POLYi  =  PSYNTHLROOTi ) 

-  Polynomial  Multiply 
POLYi  =  POLYi  *  POLYk 

*  Polynomial  Print 
Print  out  POLYi 

-  Find  Roots  Of  Polynomial 
ROOTi  =  Roots  of  POLYi 

-  Polynomial  Subtract 
POLYi  =  POLYj  ~  POLYk 


CPYPS(i.i.k)  - 


CPYSP(i.i.k)  - 


Copy  Polynomials  Into  S-Plane  Transfer  Function 
SPTFi  =  POLYi  /  POLYk 

Copy  S-Plane  Transfer  Function  Into  Polynomials 
POLYi  =  numerator  of  SPTFi 
POLYk  =  denominator  of  SPTFi 


FREQS(FAUXl)  -  S-Plane  Frequency  Response  With  User  Supplied  Function 

f 

SELCR(i)  -  Eliminate  Common  Roots  From  S-Plane  Transfer  Function 

SPTFi  =  SPTFi  with  common  roots  eliminated 
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SFREQ(i) 


-  S-Plane  Frequency  Response 

Compute  frequency  response  of  SPTFi 

SLOCK i)  -  S-Plane  Root  Locus 

Compute  root  locus  of  SPTFi 

SNORN(i)  -  Normalize  S-Plane  Transfer  Function 

SPTFi  =  SPTFi  with  coefficients  normalized 

SPADD(i,j,k)  -  S-Plane  Transfer  Function  Add 
SPTFi  =  SPTFj  +  SPTFk 

SPDIV(i,j,k)  -  S-Plane  Transfer  Function  Divide 
SPTFi  =  SPTFj  /  SPTFk 

SPEQU(i,j)  -  S-Plane  Transfer  Function  Equal 
SPTFi  =  SPTFj 

SPLDC(i)  -  S-Plane  Transfer  Function  Load#  Coefficient  Form 
SPTFi  =  POLYN  /  POLYD 

SPLDR(i)  -  S-Plane  Transfer  Function  Load,  Root  Form 

SROOTi  =  ROOTN  /  ROOTD  ,  SPTFi  =  PSYNTH( SROOTi ) 

SPMPY(i#j#k)  -  S-Plane  Transfer  Function  Multiply 
SPFTi  =  SPTFj  X  SPTFk 

SPPRN(i)  -  S-Plane  Transfer  Function  Print 

Print  out  SPTFi 

SPRTS(i)  -  Find  Roots  Of  S-Plane  Transfer  Function 
SROOTi  =  Roots  of  SPTFi 

SPSUBCi# j,k)  -  S-Plane  Transfer  Function  Subtract 
SPTFi  =  SPTFj  -  SPTFk 

STIME(i)  -  Inverse  Laplace  Transform  and  Time  Response 

Compute  time  response  of  SPTFi  by  partial  fraction  expansion 

SHMRX(i#j)  -  S-to-H  Plane  Multirate  Transform 
(slow-to-fast  sampler) 

WPTFi  =  w  plane  multirate  transform  of  SPTFj 

SHXFM(i#j)  -  S-to-M  Plane  Transform 

WPTFi  =  z  plane  transform  of  SPTFj 

SZMRXCi# j)  -  S-to-Z  Plane  Multirate  Transform 
(slow-to-fast  sampler) 

ZPTFi  =  z  plane  multirate  transform  of  SPTFj 


SZXFM(i, j) 


S-to-Z  Plane  Transform 

ZPTFi  =  z  plane  transfrom  of  ?TFj 


CPYPZ(i,j»k)  -  Copy  Polynomials  Into  Z-Plane  Transfer  Function 
ZPTFi  =  POLYj  /  POLYk 

CPYZP(i,j,k)  -  Copy  Z-Plane  Transfer  Function  Into  Polynomials 
POLYj  =  numerator  of  ZPTFi 
POLYk  =  denominator  of  ZPTFi 

FREQZ(FAUXl)  -  Z-Plane  Frequency  Response  Nith  User  Supplied  Function 

ZELCR(i)  -  Eliminate  Common  Roots  From  Z-Plane  Transfer  Function 

ZPTFi  =  ZPTFi  with  common  roots  eliminated 

ZFREQ(i)  -  Z-Plane  Frequency  Response 

Compute  frequency  response  of  ZPTFi 

ZLOCI(i)  -  Z-Plane  Root  Locus 

Compute  root  locus  of  ZPTFi 

ZMRFQ(i)  -  Z-Plane  Multirate  Frequency  Response 

Compute  multi rate  rate  frequency  response  of  ZPTFi 
by  application  of  frequency  decomposition  method 

ZMRXFM(i,j)  Z-Plane  Multirate  Transform  By  Frequency  Decomposition 
(fast-to-slow  sampler) 

ZPTFi  =  Multirate  transform  of  ZPTFj  by  application 
of  frequency  decomposition  method. 

ZNORM(i)  -  Normalize  Z-Plane  Transfer  Function 

ZPTFi  =  ZPTFj  with  coefficients  normalized 

ZPADD(i,j,k)  -  Z-Plane  Transfer  Function  Add 
ZPTFi  =  ZPTFj  +  ZPTFk 

ZPDIV(i,j,k)  -  Z-Plane  Transfer  Function  Divide 
ZPTFi  =  ZPTFj  /  ZPTFk 

ZPEQU(i>j)  -  Z-Plane  Transfer  Function  Equal 
ZPTFi  =  ZPTFj 

ZPLDC(i)  -  Z-Plane  Transfer  Function  Load*  Coefficient  Form 

ZPTFi  =  POLYN  /  POLYD 


ZPLDR(i) 


-  Z-Plana  Transfer  Function  Load,  Root  Form 
ZROOTi  =  ROOTN/  ROOTD,  ZPTFi  =  PSYNTHtZROOTi ) 

ZPMPY<i,j,k)  -  Z-Plana  Transfer  Function  Multiply 
ZPTFi  »  ZPTFj  *  ZPTFk 

ZPPRN(i)  -  Z-Plana  Transfer  Function  Print 

Print  out  ZPTFi 

ZPRTS(i)  -  Find  Roots  Of  Z-Plane  Transfer  Function 

ZROOTi  =  Roots  of  ZPTFi 

ZPSUB(i,j,k)  -  Z-Plane  Transfer  Function  Subtract 
ZPTFi  =  ZPTFj  -  ZPTFk 

ZSXFM(i)  -  Z-to-S  Root  Transformation 

Compute  "s  plane  equivalent"  of  roots  of  ZPTFi 

ZTINE(i)  -  Inverse  Z-Transform  and  Time  Response 

Compute  time  response  of  ZPTFi 

ZVCNG(i,j,n)  -  Z-to-ZN  Transform 

ZPTFi  =  ZPTFj  with  z  replaced  with  z**n 

ZMXFMCi, j)  -  Z-to-H  Plane  Transform 

HPTFi  =  Bilinear  transform  of  ZPTFi 


Copy  Polynomials  Into  W-Plane  Transfer  Function 
HPTFi  =  POLYj  /  POLYk 

CPYWP(i,j»k)  -  Copy  M-Plane  Transfer  Function  Into  Polynomials 
POLYj  -  numerator  of  WPTFi 
POLYk  -  denominator  of  HPTFi 

FREQZ(FAUXl)  -  Z-Plane  Frequency  Response  With  User  Supplied  Function 

HELCR(i)  -  Eliminate  Common  Roots  From  H-Plane  Transfer  Function 

HPTFi  =  HPTFi  with  common  roots  eliminated 

HFREQ(i)  -  H-Plane  Frequency  Response 

Compute  frequency  response  of  HPTFi 

HLOCI(i)  -  H-Plane  Root  Locus 

Compute  root  locus  of  HPTFi 

HMRFQ(i)  -  H-Plane  Multirate  Frequency  Response 

Compute  multirate  frequency  response  of  WPTFi  by 
application  of  frequency  decomposition  method 


HMRXFM(i * j ) 


NNORM(i) 


HPADD(i,j,k)  - 


WPDIV(i*j,k)  - 


WPEQUCi, j) 


WPLDCCi) 


HPLDR(i) 


HPMPYCi* j*k)  - 


HPPRN(i) 


HPRTS(i) 


HPSUB(i,j*k)  - 


WSXFM(i) 


HZXFMCi, j) 


DTERMti, j) 


DETRM(i) 


H-Plane  Multirate  Transform  By  Fraquency  Dacomposition 
(fast-to-slow  samplar) 

WPTFi  =  Multirate  transform  of  WPTFj  by  application 
frequency  decomposition  method 

Normalize  W-Plane  Transfer  Function 
WPTFi  =  WPTFi  with  coefficients  normalized 

W-Plane  Transfer  Function  Add 
HPTFi  =  WPTFj  +  NPTFk 

W-Plane  Transfer  Function  Divide 
WPTFi  =  WPTFj  /  WPTFk 

W-Plane  Transfer  Function  Equal 
WPTFi  =  WPTFj 

W-Plane  Transfer  Function  Load*  Coefficient  Form 
WPTFi  =  POLYN  /  POLYD 

W-Plane  Transfer  Function  Load*  Root  Form 
HROOTi  =  ROOTN  /  ROOTD*  WPTFi  =  PSYNTHCHROOTi) 

W-Plane  Transfer  Function  Multiply 
WPTFi  =  WPTFj  *  WPTFk 

W-Plane  Transfer  Function  Print 
Print  out  WPTFi 

Find  Roots  Of  H-Plane  Transfer  Function 
WROOTi  =  Roots  of  WPTFi 

W-Plane  Transfer  Function  Subtract 
WPTFi  =  WPTFj  -  WPTFk 

W-to-S  Root  Transformation 

Find  "s  plane  equivalent"  of  roots  of  HPTFi 

W-to-Z  Plane  Transform 

ZPTFi  -  Bilinear  transform  of  WPTFj 


Determinant  Of  Matrix  H(s)  With  Substitution  Of  Vector  £(s) 
For  Use  In  Transfer  Function  Evaluation  Via  Cramer's  Method 
POLYi  =  det  H(s)  with  column  j  replaced  with  £(s) 

Old  Version  Of  Operator  DTERM 
(No  substitution  of  £(s)  ) 

POLYi  =  det  M(s) 


V 

V 


STORE(i)  -  Store  Data  From  Current  Batch  Job 
(printout  suppressed  if  i.EQ.O) 

RESTORE(i)  -  Restore  Data  From  Old  Batch  Job  Or  Interactive  Session 
(printout  suppressed  if  i.EQ.O) 


Examples  are  presented  to  demonstrate  the  utility  of  typical  LCAP2  opera¬ 
tors.  Examples  1  through  12  were  prepared  to  be  executed  sequentially  in  one 
batch  job.  Example  12  demonstrates  the  use  of  the  STORE  operator  used  to  save 
transfer  function*  polynomial  and  matrix  data  for  a  restart  capability.  Example 
13  demonstrates  the  RESTORE  operator  which  restores  data  stored  from  a  previous 
batch  job. 

Each  of  the  examples  begins  with  a  statement  of  the  problem  and  is  followed 
by  the  user's  FORTRAN  code  and  the  printed  output  of  the  program.  To  differen¬ 
tiate  between  FORTRAN  code  and  comments*  the  FORTRAN  code  is  in  upper  case  and 
the  comments  are  in  lower  case. 

In  Examples  1*5*6*7,9  and  10*  which  have  requests  for  both  printer  and 
hardcopy  plots*  only  the  printer  plots  are  shown.  The  hardcopy  plot  file  genei — 
ated  by  these  examples  is  processed  by  the  HARDCPY  program  after  the  LCAP2  pro¬ 
gram  has  been  executed.  The  hardcopy  plots  from  these  examples  are  presented  in 
Appendix  F. 

The  examples  in  this  section  are  essentially  the  same  as  those  presented  in 
the  Interactive  LCAP2  User's  Guide,  Ref.  3,  so  that  differences  in  usage  between 
both  versions  of  LCAP2  can  be  compared.  It  should  be  noted  that  Batch  LCAP2  can 
include  FORTRAN  expressions  as  part  of  the  data  entry,  but  was  not  used  in  these 
examples  since  the  problems  all  used  numeric  data. 


Pr obi  mii  Coo  put*  froquoncy  rosponso  of 


25 


2 

•  ♦  5s  ♦  25 

between  .10  to  100.  ra d/*oc. 


Data  will  ba  loadad  in  coafficiant  form  using  tha  operator  SPLDC.  Polynomial 
arrays  POLYN  and  P01YD  are  used  with  SPLDC. 


Tha  FORTRAN  coda  for  this  example  isi 


EXAMPLE  1 

LOAD  IN  DATA  USINO  LCAP2  OPERATOR  SPLDC  AND  POLYNOMIAL  ARRAYS 
POLYN  AND  POLYD. 


P0LYN(1)=0. 

P0LYNC2)=25. 

P0LYD(1)S2. 

P0LYD(2)S25. 

P0LYD(3)S5. 

P0LYD(A)*1. 

CALL  SPLDC(l) 


"dag .  of  num . " 

"coaff.  of  sXXO* 

"dag.  of  denom." 

"coaff.  of  sXXO" 

"coaff.  of  sxxi" 

"coaff.  of  sxx2" 

"load  coafficiant  data  into  SPTF1" 


ENTER  FREQUENCY  RESPONSE  PARAMETERS  FOR  USE  HITH  SFREQ 


FAUT0*1 

RAD=1 

N0MEG=3 

0MEGA(  1 )  ='.  1 

0MEGA(2)  =  1 . 

OMEGA(3)-100. 

FDLAY»0 

FNIC0*1 

FB0DE=1 

CYCLE= 0 


FNYQS=1 

N9DB=1 

GRAFP=1 

FILM=1 


".NE.O  (preset=l)  for  auto.  fraq.  selection  mode" 
".NE.O  (praset=l)  for  rad/sec*  otherwise  hz" 
"number  of  values  of  OMEGA  to  be  entered" 
"0MEGA(l)=first  frequency  value  to  be  used" 

"user  specified  frequency  value" 

"0MEGA(N0MEG)=last  frequency  value  to  be  used" 
"time  delay  (preset=0)" 

".NE.O  (presets0)  for  Nichols  plot" 

".NE.O  (preset=l)  for  Bode  plot" 

".EQ.O  for  auto,  selection  of  2  or  3  cycle  for 
Bode  plots  (1  cycle  not  available)" 

".NE.O  (preset=0)  for  Nyquist  plot" 

".NE.O  Cpreset=0)  for  hardcopy  Nyquist  plot  in  db" 
".NE.O  (preset=l)  for  printer  plot" 

".NE.O  (presets0)  for  hardcopy  (high  resolution) 


uneuA 
RAD/SEC 
.1000 
.1200 
.  1400 
.1728 
.2074 
.2488 
.2986 
.3583 
.4300 
.5160 
.6192 
.7430 
.8916 
L.070 
L  .284 
!.084 
'.724 
(.364 
(.844 
1.324 
1.804 
i.000 
i .  480 
.120 
.920 
.880 
.840 
0.12 

1.40 
3.00 
4.92 

6.84 

9.40 

1.96 

5.16 
9.00 

2.84 

7.96 
3.08 
9.48 

7.16 

4.84 
5.08 
5.32 
1.12 
10.0 


REAL 

IMAGINARY 

OB 

PHASE 

.  100E+01 

- . 200E-01 

.002 

-1.15 

.  100E+01 

-.240E-01 

.003 

-1.38 

.  100E+01 

- . 288E-01 

.004 

-1.65 

. 100E+01 

- . 346E-01 

.005 

-1.98 

. 100E+01 

-.415E-01 

.007 

-2.38 

. 100E+01 

-.499E-01 

.011 

-2.86 

. 100E+01 

-.599E-01 

.015 

-3.43 

. 100E+01 

-.720E-01 

.022 

-4.12 

. 100E+01 

-.866E-01 

.032 

-4.95 

. 100E+01 

-. 104E+00 

.046 

-5.95 

. 100E+01 

- . 126E+00 

.066 

-7.17 

. 100E+01 

- . 152E+00 

.095 

-8.64 

•999E+00 

- . 184E+00 

.136 

-10.44 

.998E+00 

“ . 224E+00 

.194 

-12.64 

. 995E+00 

- . 274E+00 

.276 

-15.37 

. 965E+00 

- . 487E+00 

.673 

-26.77 

. 889E+00 

-.688E+00 

1.017 

-37.77 

.728E+00 

~ . 894E+00 

1.236 

-50.87 

. 539E+00 

- .  101E+01 

1.202 

-61.99 

. 311E+00 

“ . 107E+01 

.907 

-73.74 

. 828E-01 

- . 103E+01 

.320 

-85.42 

. 225E-04 

- . 100E+01 

.000 

-90.00 

- • 162E+00 

-.883E+00 

-.940 

-100.40 

- . 285E+00 

-.701E+00 

-2.421 

-112.15 

- .  332E+00 

-.503E+00 

-4.399 

-123.48 

- .  317E+00 

“ . 336E+00 

-6.707 

-133.27 

- .  273E+00 

~ . 231E+00 

-8.834 

-140.25 

- .  226E+00 

- .  148E+00 

-11.362 

-146.83 

- .  184E+00 

-.999E-01 

-13.584 

-151.50 

- .  144E+00 

-.651E-01 

-16.014 

-155.71 

- .  111E-00 

-.418E-01 

-18.536 

-159.32 

-.874E-01 

- . 285E-01 

-20.731 

-161.96 

-.661E-01 

-. 183E-Q1 

-23.275 

-164.57 

“ . 517E-01 

-.124E-01 

-25.488 

-166.50 

-.394E-01 

-.816E-02 

-27.902 

-168.31 

- .  297E-01 

-.528E-02 

-30.410 

-169.92 

- . 232E-01 

-.361E-02 

-32.598 

-171.14 

- .  173E-01 

-.232E-02 

-35.140 

-172.37 

-.135E-01 

-.158E-02 

-37.354 

-173.29 

-.102E-01 

-.104E-02 

-39.774 

-174.17 

-.765E-02 

- . 674E-Q2 

-42.292 

-174.96 

-.595E-02 

-. 461E-03 

-44.489 

-175.56 

-.443E-02 

-.297E-03 

-47.043 

-176.17 

-.343E-02 

-.202E-03 

-49.268 

-176.63 

-.260E-02 

-.133E-03 

-51.700 

-177.08 

-.250E-02 

-.125E-03 

-52.030 

-177.13 

PHASE 

MARGIN 

178.85 
178.62 

178.35 
178.02 

178.62 
177.14 

176.57 
175.88 
175.05 
174.05 
172.83 

171.36 
169.56 

167.36 

164.63 

153.23 

142.23 
129.13 
118.01 
106.26 

94.58 

90.00 

79.60 

67.85 
56.52 
46.73 
39.75 
33.17 

28.57 
24.29 
20.68 
18.04 
15.43 
13.50 
11.69 
10.08 

8.86 

7.63 

6.71 

5.83 
5.04 
4.44 

3.83 
3.37 
2.92 
2.87 


NICHOLS  PLOT  CMAGN.  VS  PHASE) 

EXAMPLE  1  S  PLANE  FREQUENCY  RESPONSE 


10/31/83 


16.0 
14 . 0 
12.0 
10.0 
8.00 
6.00 
4.00 
2.00 

0  - 

-2.00 

-4.00 

-6.00 

-8.00 

-10.0 

-12.0 

-14.0 

-16.0 

-18.0 

-20.0 

-22.0 

-24.0 

-26.0 

-28.0 

-30.0 

-32.0 

-34.0 

-36.0 

-38.0 

-40.0 

-42.0 

-44.0 

-46.0 

-48.0 

-50.0 

-52.0 

-54.0 

-56.0 

-58.0 

-60.0 

-62.0 

-64.0 

-66.0 

-68.0 

I  I  I  I  I 
-360.000  -288.000  -216.000 
FILM  PLOT  COMPLETED 


I 

I 

I 

I 


I 

I 

I 

XXXXXXX  I 

- xxxxxx - xxxxxxx 

XXX  I 

XX  I 

XX  I 

X  I 

X  I 

XX  I 

X  I 

X  I 

XX  I 

X  I 

X  I 

XX  I 

X  I 

X  I 

X  I 

XX  I 

X  I 

X  I 

X  I 

X  I 

X  I 

X  I 

X  I 

X  I 

X  I 

X  I 


I 

I 

I 

I 

I 

I 

I 

I 

I  I  I  I  I  I 
-144.000  -72.000  0.000 


BODE  PLOT  (HAGN.  VS  FREQ.) 

EXAMPLE  1  S  PLANE  FREQUENCY  RESPONSE 


10/31/83 


16.0 

14.0 

12.0 

10.0 

6.00 

4.00 

2.00 

0. 

-2.00 

-4.00 

-6.00 

-8.00 

-10.0 

-12.0 

-14.0 

-16.0 

-18.0 

-20.0 

-22.0 

-24.0 

-26.0 

-28.0 

-30.0 

-32.0 

-34.0 

-36.0 

-38.0 

-40.0 

-42.0 

-44.0 

-46.0 

-48.0 

-50.0 

-52.0 

-54.0 

-56.0 

-58.0 

-60.0 

-62.0 

-64.0 

-66.0 

-68.0 


I  xxxx 

XX-XX-XX-XXX-XX-XXXXXXXXXX - XX - 

I  X 

I  X 

I  XX 

I  XX 

I  XX 

I  XX 

I  X 

I  X 

I  X 

I  XX 

I  XX 

1  XX 

I  XX 

I  XX 

I  X 

I  X 

I  X 

I  XX 

I  XX 

I  XX 

I  XX 

I  XX 

I  X 

I  XX 

I  X 

I  X 


I  I 
. 100E+00 


I  I  I 

1.00 


III 

10.0 


I  I  I 

100. 


*«  'v 


BODE  PLOT  (HA6N.  VS  FREQ.) 

EXAMPLE  1  S  PLANE  FREQUENCY  RESPONSE 


30.0 

20.0 

10.0 

0. 


xx-xx-xx-xxx 


-10.0 

-20.0 

-30.0 

-40.0 

-50.0 

-60.0 

-70.0 

-80.0 

-90.0 

-100. 

-110. 

-120. 

-130. 

-140. 

-150. 

-160. 

-170. 

-180. 

-190. 

-200. 

-210. 

-220. 

-230. 

-240. 

-250. 

-260. 

-270. 

-280. 

-290. 

-300. 

-310. 

-320. 

-330. 

-340. 

-350. 

-360. 

-370. 

-380. 

-390. 


xx  xxxx 

xxxx 

XX 

XXX 

XX 

XX 

XX 

X 

XX 

XX 

X 

X 

XX 

XXX 

XXX 

xxxx 

xxxxxxxxxx 

xxxxx 


I  I  I  I  I  I  I  I  I  I 
. 100E+00  1.00  10.0  100 

FILM  PLOT  COMPLETED 


10/31/83 


I 


x 


10  VO  sr  N 


NYQUIST  PLOT 

EXAHPLE  1  S  PLANE  FREQUENCY  RESPONSE 


10/31/83 


iMPLE  2  ROOTS  OF  A  POLYNOMIAL 


Problem:  Compute  the  roots  of  the  polynomial 

3  2  1 

x  +  13x  +  38x  +  34 

The  FORTRAN  code  for  this  example  is: 


EXAMPLE  2 

LOAD  IN  DATA 

USING  LCAP2  OPERATOR  PRTS  AND  POLYNOMIAL  ARRAY  POLY 

P0LY(1)=3 

"deg.  of  polynomial" 

POLY ( 2) =34 . 

"coeff.  of  xXXO" 

P0LY(2)=38. 

"coeff.  of  xxxi" 

P0LY(2)=13. 

"coeff.  of  xXX2" 

P0LY(2)=1. 

"coeff.  of  xXX3" 

CALL  PLDC(l) 

"load  coefficient  data  into  P0LY1" 

CALL  PRTS(l) 

"find  roots  of  P0LY1" 

The  printer  output  for  this  case  is: 


DEGREE  OF  POLY  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 
34.  38.  13.  1. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X  PLDC  -  POLYNOMIAL  LOAD  IN  COEFFICIENT  FORM  * 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DEGREE  OF  P0LY1  IS  3 
34.  38.  13.  1. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


DEGREE  OF  P0LY1  IS  3 
34.  38.  13.  1. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
x  PRTS  -  FIND  ROOTS  OF  P0LY1  * 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


CM  W 


1  -1.8445105  -.49938380  1.9109168  .96524896 

-1.8445105  .49938380  1.9109168  .96524896 

-9.3109790  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  34.000000 


DEGREE  OF  POLY  IS  3 
34.  38.  13.  1. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


I 


I 


Problem >  Multiply  the  following  two  e-plane  transfer  functions 


30C  -  +  1  ) 

25  10 

_  and  _ 

2  s  s  s 

s  +  5s  +  25  70(s)(  -  +  1  )C  -  +  1  )(  -  +  1  ) 

7  1  +  32  1-32 

The  first  transfer  function  is  the  same  as  SPTF1  from  Example  1.  The  second 
transfer  function  is  given  in  root  form.  When  transfer  function  data  is 
expressed  in  root  form*  gains  associated  with  the  numerator  and  denominator  must 
be  specified  to  uniquely  define  the  transfer  function.  The  gains  used  by  LCAP2 
correspond  to  the  low  order  non-zero  coefficient  of  the  numerator  and  denomina¬ 
tor  polynomials  if  the  root  form  were  expanded  out.  For  this  example*  these  two 
gains  would  be  30  and  70.  (A  distinction  is  made  for  the  low  order  non-zero  coef¬ 
ficient  since  the  low  order  coefficient  can  be  zero  as  is  that  of  the 
denominator  in  this  example.) 

The  FORTRAN  code  for  this  example  ist 


EXAMPLE  3 

FIRST  TRANSFER  FUNCTION  IS  THE  SAME  AS  SPTF1  OF  EXAMPLE  1 

LOAD  IN  ROOT  DATA  FOR  THE  SECOND  TRANSFER  FUNCTION  USING 
SUBROUTINE  SPLDR  AND  POLYNOMIAL  ROOT  ARRAYS  ROOTN  AND  ROOTD 


R00TN(1)=(1 . ,30. ) 

R00TN(2)=(-10. ,0. ) 
R00TD(1)=(4. ,70. ) 

ROOTDC 2)=(-l.,2.) 
R00TD(3)=(-1 . ,-2. ) 
ROOTD<4)=(-7.,0.) 
R0OTD(5)  =  (0.  ,0.  ) 
CALL  SPLDR(2) 


"real  part  -  no.  of  roots  and  imag. 

part  =  low  non-zero  gain  of  num.*" 
"first  root" 

"real  part  =  no.  of  roots  and  imag. 

part  =  low  non-zero  gain  of  denom. 
"first  root" 

"second  root" 

"third  root" 

"fourth  root" 

"load  root  data  into  SPTF2" 


MULTIPLY  THE  TWO  TRANSFER  FUNCTIONS 

CALL  SPMPY(3,1,2)  "multiply  SPTF1  and  SPTF2  and  store 

in  SPTF3" 


The  printer  output  for  this  example  ist 


THE  ROOTS  OF  ROOTN  ARE 


NO.  REAL  IMAG.  OMEGA  ZETA 

1  -10.000000  0. 

THE  ROOTS  OF  ROOTD  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -1.0000000  2.0000000 

2  -1.0000000  -2.0000000 

3  -7.0000000  0. 

4  0.  0. 

XXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  SPLDR  -  LOAD  TRANSFER  FUNCTION  IN  ROOT  FORM  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THE  NUMERATOR  ROOTS  OF  SR00T2  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -10.000000  0. 


LON  ORDER  NON  ZERO  COEFFICIENT  =  30.000000 

THE  DENOMINATOR  ROOTS  OF  SR00T2  ARE 


NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-1.0000000 

2.0000000 

2.23606798 

.447213595 

2 

-1.0000000 

-2.0000000 

2.23606798 

.447213595 

3 

-7.0000000 

0. 

4 

0. 

0. 

LON  ORDER  NON  ZERO  COEFFICIENT  =  70.000000 

DEGREE  OF  NUMERATOR  OF  SPTF2  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

30.  3. 

DEGREE  OF  DENOMINATOR  OF  SPTF2  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  70.  38.  18.  2. 


BODE  GAIN  a 


42857143 


DEGREE  OF  NUMERATOR  OF  SPTF1  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 
25. 

DEGREE  OF  DENOMINATOR  OF  SPTF1  IS  2  (COEFFICIENTS  IN  ASCENDING  ORDER) 
25.  5.  1. 

BODE  GAIN  =  1.000000 


THE  NUMERATOR  ROOTS  OF  SR00T2  ARE 

IMAG.  OMEGA  ZETA 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  30.000000 

THE  DENOMINATOR  ROOTS  OF  SR00T2  ARE 

IMAG.  OMEGA  ZETA 

2.0000000 
-2.0000000 
0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  70.000000 

DEGREE  OF  NUMERATOR  OF  SPTF2  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

30.  3. 

DEGREE  OF  DENOMINATOR  OF  SPTF2  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  70.  38.  18.  2. 

BODE  GAIN  =  .42857143 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  SPTF3  =  SPTF1  X  SPTF2  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DEGREE  OF  NUMERATOR  OF  SPTF3  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 
750.  75. 

DEGREE  OF  DENOMINATOR  OF  SPTF3  IS  6  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  1750.  1300.  710.  178.  28.  2. 


NO.  REAL 

1  -10.000000 

NO.  REAL 

1  -1.0000000 

2  -1.0000000 

3  -7.0000000 


BODE  GAIN 


42857143 


Not*  that  if  th*  second  transfer  function  war*  expressed  in  a  different  form* 
the  numerator  and  denominator  non-zero  coefficients  would  not  necessarily 
appear  explicitly.  For  example*  if  th*  root  locus  representation  were  used* 
i  .*.  * 


3(  s  +  10  ) 


2(  s  ♦  7  >t  s  ♦  1  +  52  >(  s  +  1  -  *2) 


them 

numerator  low  order  non-zero  coefficient  =  (3)(10)  =  30 
denominator  low  order  non-zero  coefficient  =  (2)(7)(l+32)(l-;j2)  =  70 
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CALL  SPLDC(5>  "load  HCs)  data  into  SPTF5" 

C  . 

C  HULTIPLY  G(S)  AND  H(S) 

C 

CALL  REMARK* ( 40HG( S ) XH ( S )  IS  IN  SPTF6  ) 

CALL  SPMPY(6,*,5)  "G(s)xH(s)  is  in  SPTF6" 

C  . 

C  LOAD  IN  UNITY  TRANSFER  FUNCTION 

C 

POLYN( 1 )=0 
P0LYN(2)=1 
POLYDC1)=0 
POLYD(2)=l 
CALL  SPLDC(7) 

C  . 

C 

CALL  REMARK* (40 HI  +  G(S)»G(S)  IS  IN  SPTF8  ) 

CALL  SPADD(8,7,6)  "1  +  GCS)XH(S)  is  in  SPTF8" 

C  . 

C  FORM  CLOSED  LOOP  TRANSFER  FUNCTION  BY  DIVIDING  G(S)  BY  G(S)+H(S) 

C 

CALL  SPDIV(9,*,8)  "closed  loop  transfer  function  is  in  SPTF9" 

C  . 

C  ELIMINATE  COMMON  ROOTS  BETWEEN  NUMERATOR  AND  DENOMINATOR 
C 

CALL  REMARK*! 40HCLOSED  LOOP  T.  F.  IS  IN  SPTF9  ) 

CALL  SELCRC9)  "eliminate  common  roots  of  SPTF4" 

C  . 


The  printer  output  for  this  example  isi 


G(S)  IS  IN  SPTF* 


DEGREE  OF  POLYN 
15. 

IS 

0 

(COEFFICIENTS 

IN 

ASCENDING 

ORDER) 

DEGREE  OF  POLYD 
0.  5.  6.  1. 

IS 

2 

(COEFFICIENTS 

IN 

ASCENDING 

ORDER) 

"deg.  of  num." 

"coeff.  of  sXXO" 

"deg.  of  denom." 

"coeff.  of  sXXQ" 

"unity  transfer  function  is  in  SPTF7" 


XXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  SPLDC  -  LOAD  TRANSFER  FUNCTION  IN  COEFFICIENT  FORM  x 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


DEGREE  OF  NUMERATOR  OF  SPTF*  IS  0 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


H(S)  IS  IN  SPTF5 


DEGREE  OF  POLYN  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1.5  1. 

DEGREE  OF  POLYD  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

0.  1. 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  SPLDC  -  LOAD  TRANSFER  FUNCTION  IN  COEFFICIENT  FORM  * 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

DEGREE  OF  NUMERATOR  OF  SPTF5  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1.5  1. 

DEGREE  OF  DENOMINATOR  OF  SPTF5  IS  I  (COEFFICIENTS  IN  ASCENDING  ORDER) 
15.  1. 

BODE  GAIN  =  . 10000000E+00 


CP=  5.33 

G(S)*H(S)  IS  IN  SPTF6 

DEGREE  OF  NUMERATOR  OF  SPTF4  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

15. 

DEGREE  OF  DENOMINATOR  OF  SPTF4  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  5.  6.  1. 

BODE  GAIN  =  3.0000000 

DEGREE  OF  NUMERATOR  OF  SPTF5  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1.5  1. 

DEGREE  OF  DENOMINATOR  OF  SPTF5  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 
15.  1. 

BODE  GAIN  =  . 10000000E+00 


KXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  SPTF6  =  SPTF4  x  SPTF5  * 

XXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  NUMERATOR  OF  SPTF6  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 
22.5  15. 

DEGREE  OF  DENOMINATOR  OF  SPTF6  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  75.  95.  21.  1. 

10DE  GAIN  =  .30000000 


CP=  5.35 

DEGREE  OF  POLYN  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

DEGREE  OF  POLYD  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

KXXXKXKXXXKXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXX 

X  SPLDC  -  LOAD  TRANSFER  FUNCTION  IN  COEFFICIENT  FORM  X 

KXXXXKXKXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  NUMERATOR  OF  SPTF7  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

DEGREE  OF  DENOMINATOR  OF  SPTF7  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

BODE  GAIN  -  1.00000000 


CP=  5.36 

1  ♦  G(S)*H(S)  IS  IN  SPTF8 

DEGREE  OF  NUMERATOR  OF  SPTF7  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

DEGREE  OF  DENOMINATOR  OF  SPTF7  IS  0  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1. 

BODE  GAIN  »  1.00000000 


DEGREE  OF  NUMERATOR  OF  SPTF6  IS  1 
22.5  15. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


DEGREE  OF  DENOMINATOR  OF  SPTF6  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  75.  95.  21.  1. 

BODE  GAIN  -  .30000000 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  SPTF8  =  SPTF7  +  SPTF6  X 

XXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  NUMERATOR  OF  SPTF8  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
22.5  90.  95.  21.  1. 

DEGREE  OF  DENOMINATOR  OF  SPTF8  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  75.  95.  21.  1.0000000 

BODE  GAIN  =  .30000000 


CP=  5.38 


DEGREE  OF  NUMERATOR  OF  SPTF4  IS  0 
15. 

DEGREE  OF  DENOMINATOR  OF  SPTF4  IS  3 
0.  5.  6.  1. 

BODE  GAIN  =  3.0000000 

DEGREE  OF  NUMERATOR  OF  SPTF8  IS  4 
22.5  90.  95.  21.  1. 

DEGREE  OF  DENOMINATOR  OF  SPTF8  IS  4 
0.  75.  95.  21.  1.0000000 

BODE  GAIN  =  .30000000 


(COEFFICIENTS  IN  ASCENDING  ORDER) 

(COEFFICIENTS  IN  ASCENDING  ORDER) 

(COEFFICIENTS  IN  ASCENDING  ORDER) 

(COEFFICIENTS  IN  ASCENDING  ORDER) 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X  SPTF9  =  SPTF4  /  SPTF8  * 
XXXXXXXKXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXKXX 


DEGREE  OF  NUMERATOR  OF  SPTF9  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 

0.  1125.  1425.  315.  15. 

DEGREE  OF  DENOMINATOR  OF  SPTF9  IS  7  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  112.5  585.  1037.5  765.  226.  27.  1. 


BODE  GAIN 


10.000000 


CLOSED  LOOP  T.F.  IS  IN  SPTF9 


DEGREE  OF  NUMERATOR  OF  SPTF9  IS  A  (COEFFICIENTS  IN  ASCENDING  ORDER) 

0.  1125.  1425.  315.  15. 

DEGREE  OF  DENOMINATOR  OF  SPTF9  IS  7  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  112.5  585.  1037.5  765.  226.  27.  1. 

BODE  GAIN  =  10.3.0000000 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  SELCR  -  ELIMINATE  COMMON  ROOTS  OF  S-PLANE  X 

X  TRANSFER  FUNCTION  9  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


THE  COMMON  ROOTS  ELIMINATED  ARE 


NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-1.0000000 

0. 

2 

-5.0000000 

0. 

3 

0. 

0. 

XXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  SPTF9  =PSYNTH( SR00T9  )  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


THE  NUMERATOR  ROOTS  OF  SROOT9  ARE 
NO.  REAL  IMAG.  OMEGA  ZETA 

1  -15.000000  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  1125.0000 

THE  DENOMINATOR  ROOTS  OF  SROOT9  ARE 


NO.  REAL  IMAG.  OMEGA  ZETA 

1  -.755541434  0. 

2  -.416857608  0. 

3  -4.73265645  0. 

4  -15.0949445  0. 


LOW  ORDER  NON  ZERO  COEFFICIENT  =  112.50000 


DEGREE  OF  NUMERATOR  OF  SPTF9  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 
1125.  75. 

DEGREE  OF  DENOMINATOR  OF  SPTF9  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
112.5  450.  475.  105.  5. 


BODE  GAIN 


10.000000 


Problem i  Compute  tho  root  locus  for  the  following  transfer  function 

s 

1 . 8(  -  +1  ) 

1 


s<  -  ♦  1  )(  -  +  1  )(  -  ♦  1)(  -  +  1  ) 

2  3  5  -  jl  5  +  jl 

by  varying  the  nominal  gain  from  .125  to  2. 

The  FORTRAN  code  for  this  example  isi 


EXAMPLE  5 

LOAD  IN  ROOT  DATA 

FOR  TRANSFER  FUNCTION 

ROOTN( 1)=(1..1.8) 

"real  part  =  no.  of  roots  and  imag. 
part  =  low  non-zero  gain  of  nun.," 

ROOTN(2)=(-1.,0.) 

"first  root" 

R00TD(1)  =  (5. »1 .  ) 

"real  part  =  no.  of  roots  and  imag. 
part  =  low  non-zero  gain  of  denom.," 

R00TDC2)=(-5.,-l. 

)  "first  root" 

R00TD(3)=<-5.,1.) 

"second  root" 

ROOTDt  4)a<-3.»0.) 

"third  root" 

R00TD(5)=(-2. *0.  ) 

"fourth  root" 

R00TDC6 )  =  ( 0 . , 0  .  ) 

"fifth  root" 

CALL  SPLDR(IO) 

"data  stored  in  SPTF10" 

ENTER  ROOT  LOCUS 

PARAMETERS  FOR  USE  WITH  SLOCI 

NL0CI=2 

"number  of  values  of  KGAIN  to  be  entered 

KGAINC 1 )= . 125 

"KGAIN(l)-first  gain  value  to  be  used 

KGAIN(2)=2. 

"KGAINtNLOCI )=last  gain  value  to  be  used 

KFLG=0 

".EQ.O  to  increment  gain  by  multiplying  by  KDELT, 
.NE.O  to  increment  gain  by  adding  by  KDELT" 

KDELT=2 

"value  for  changing  gains  (preset=l.E4  so  that  no 
additional  gains  are  computed  by  the  program" 

GRAFP=1 

".NE.O  (preset=l)  for  printer  plot 

FILM=1 

".NE.O  (preset=0)  for  hardcopy  (high  resolution) 
plot" 

RLXMN=-9 

"min.  x  axis  for  plot" 

RLXMX=1 

"max.  x  axis  for  plot" 

"auto,  scaling  of  x  axis  if  RLXMN=RLXMX" 

RLYMN=-1 

"min.  y  axis  for  plot" 

RLYMX=9 

"max.  y  axis  for  plot" 

"auto,  scaling  of  y  axis  if  RLYMN=RLYMX" 

Ul 


CALL  HEADIN4C 1 » 40HEXAMPLE  5  S  PLANE  ROOT  LOCUS 
CALL  SLOCI(IO)  "compute  root  loci  of  SPTF10" 


The  printer  output  for  this  example  isi 


THE  ROOTS  OF  ROOTN  ARE 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

-1.0000000 

0. 

THE  ROOTS  OF  ROOTD  ARE 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

-5.0000000 

-1.0000000 

2 

-5.0000000 

-1.0000000 

3 

-3.0000000 

0. 

4 

-2.0000000 

0. 

5 

0. 

0. 

XXKXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

*  SPLDR  =  LOAD  TRANSFER  FUNCTION  IN  ROOT  FORM 

X 

XXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THE  NUMERATOR  ROOTS  OF  SROOTIO 

ARE 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

-1.0000000 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  = 

1.8000000 

THE  DENOMINATOR  ROOTS  OF  SROOTIO 

ARE 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

-5.0000000 

-1.0000000 

2 

-5.0000000 

-1.0000000 

3 

-3.0000000 

0. 

4 

-2.0000000 

0. 

5 

0. 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  = 

1.0000000 

1 

DEGREE 

OF  NUMERATOR 

OF  SPTF10  IS  1  (COEFFICIENTS 

IN 

ASCENDING  ORDER) 

BODE  GAIN  =  1.8000000 


THE  NUMERATOR  ROOTS  OF  SROOTIO  ARE 


IMAG. 


OMEGA 


-1.0000000 


LOW  ORDER  NON  ZERO  COEFFICIENT  =  1.8000000 


THE  DENOMINATOR  ROOTS  OF  SROOTIO  ARE 


IMAG. 


OMEGA 


-5.0000000 

-5.0000000 

-3.0000000 

-2.0000000 

0. 


-1.0000000 

-1.0000000 

0. 

0. 

0. 


LOW  ORDER  NON  ZERO  COEFFICIENT 


1.0000000 


DEGREE  OF  NUMERATOR  OF  SPTF10  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

1.8  1.8 

DEGREE  OF  DENOMINATOR  OF  SPTF10  IS  5  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  1.  12.17948717949  .525641025641  .09615384615385  .006410256410256 


BODE  GAIN  = 


1.8000000 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X  SLOCI  -  ROOT  LOCUS  OF  S-PLANE  TRANSFER  X 
X  FUNCTION  10  x 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THE  OPEN  LOOP  ZEROES  ARE 


IMAG. 


OMEGA 


-1.000000 


iW  r 


,C 


,1? 


VJ. 


THE  OPEN 

LOOP  POLES  ARE 

REAL 

IMAG. 

OMEGA 

5.0000000 

-1.0000000 

5.0000000 

1.0000000 

2.0000000 

0. 

3.0000000 

0. 

CLOSED 

LOOP  POLES 

FOR  GAIN  =  .1250000 

(GAIN  NO. 

1)  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-1.7947341 

-1.1077310 

2.1090611 

.85096357 

2 

-1.7947341 

1.1077310 

2.1090611 

.85096357 

3 

-5.5892217 

1.6614254 

5.8309290 

.95854737 

4 

-5.5892217 

-1.6614254 

5.8309290 

.95854737 

5 

-.23208832 

0. 

CLOSED 

LOOP  POLES  FOR 

GAIN  =  .2500000 

(GAIN  NO. 

2)  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-1.4171985 

-1.4418858 

2.0217532 

.70097500 

2 

-1.4171985 

1.4418858 

2.0217532 

.70097500 

3 

-5.8577544 

1.9606112 

6.1771582 

.94829277 

4 

-5.8577544 

-1.9606112 

6.1771582 

.94829277 

5 

-.45009415 

0. 

CLOSED 

LOOP  POLES 

FOR  GAIN  =  .5000000 

(GAIN  NO. 

3)  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-.94735061 

-1.9000842 

2.1231564 

.44619915 

2 

-.94735061 

1.9000842 

2.1231564 

.44619915 

3 

-6.1975818 

-2.3343777 

6.6226384 

.93581763 

4 

-6.1975818 

2.3343777 

6.6226384 

.93581763 

5 

-.71013518 

0. 

7  W  '?.' 

CLOSED 

LOOP  POLES  FOR 

GAIN  -  1.0000000 

(GAIN 

NO. 

4)  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-.44922549 

-2.4680974 

2.5086468 

.17907084 

2 

-.44922549 

2.4680974 

2.5086468 

.17907084 

3 

-6.6183282 

-2.7904455 

7.1825382 

.92144699 

4 

-6.6183282 

2.7904455 

7.1825382 

.92144699 

5 

-.86489263 

0. 

CLOSED 

LOOP  POLES  FOR 

GAIN  -  2.0000000 

(GAIN 

NO. 

5)  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

.10025104 

-3.1088677 

3.1104837 

-.32230049E-01 

2 

.10025104 

3.1088677 

3.1104837 

-.32230049E-01 

3 

-7.1323100 

-3.3395771 

7.8754442 

.90563907 

4 

-7.1323100 

3.3395771 

7.8754442 

.90563907 

5 

-.93588209 

0. 

FILM  PLOT  COMPLETED 

ROOT  LOCUS 

EXAMPLE  5  S  PLANE  ROOT  LOCUS 


10/31/85 


10.2  I 

9.90  I 

9.60  I 

9.30  I 

9.00  I 

8.70  I 

8.90  1 

8.10  I 

7.80  I 

7.50  I 

7.20  I 

6.90  I 

6.60  I 

6.30  I 

6.00  I 

5.70  I 

5.40  I 

5.10  I 

4.80  I 

4.50  I 

4.20  I 

3.90  I 

3.60  I 

3.30  *  I 

3.00  IK 

2.70  K  I 

2.40  K  *1 

2.10  X  I 

1.80  K  K  I 

1.50  X  I 

1.20  X  I 

.900  I 

.600  I 

.300  I 

0.  XXXX-XI  — 

-.500  I 

-1.00  I 

I  I  I  I  I  I  I  I  I  I 
-9.00  -8.00  -7.00  -6.00  -5.00  -4.00  -3.00  -2.00  -1.00  0.00 


.1 


r 


Problem!  Find  the  inverse  Laplace  transform  and  the  step  response 
of  the  following  transfer  function 

432s  +  4320 

-  -  2  T~ 

s  +  35s  +  345s  +  1008s  +  2160 

Plot  the  response  in  increments  of  .05  seconds  from  0  to  5  seconds. 

The  inverse  Laplace  transform  is  computed  by  the  partial  fraction  method.  The 
algorithm  used  to  compute  the  partial  fraction  expansion  requires  that  (1)  there 
are  no  multiple  poles  other  than  those  at  the  origin  and  (2)  the  degree  of  the 
numerator  must  not  be  greater  than  the  number  of  non-zero  poles  of  the  denomina¬ 
tor.  Up  to  5  poles  at  the  origin  are  allowed. 

The  FORTRAN  code  for  this  example  is> 


EXAMPLE  6 

LOAD  COEFFICIENT  DATA  FOR  TRANSFER  FUNCTION 


POLYNt 1 ) =1 . 
P0LYN(2) =4320 . 
P0LYN(3)=432. 
P0LYD( 1 )=4 . 
P0LYD(2)=2160. 
P0LYD(3)=1008 . 
P0LYD(4)=345. 
P0LYD(5)=35. 
P0LYD(5)  =  1 . 
CALL  SPLDC(ll) 


"deg.  of  num." 
"coeff.  of  sXXO" 
"coeff.  of  sXXl" 
"deg.  of  denom." 
"coeff.  of  sX*0" 


"coeff. 

"coeff. 


of  sxxi" 
of  SXK2" 


"coeff.  of  sXX3" 

"coeff.  of  sx*4" 

"load  coefficient  data  into  SPTF11" 


ENTER  PARAMETERS  FOR  TIME  RESPONSE 


TSTEP=1 

TMAGN=1 . 

TZERO=0 . 

TEND=5. 

TDELT=.05 

GRAFP=1 

FILM=1 


".NE.O  (preset=l)  for  step  response( 

. EQ.O  for  impulse  response" 

"magnitude  of  the  input" 

"starting  time  for  evaluating  the  response" 

"end  time  for  evaluating  the  response" 

"time  increment  for  evaluating  the  response" 
".NE.O  (preset=l)  for  printer  plot 
".NE.O  Cpreset=0)  for  higher  resolution  electro 
static  plot" 


CALL  HEADIN5( 1 , 

+50HEXAMPLE  6  INVERSE  LAPLACE 
CALL  HEADINK6*  10HNSE  ) 


rRANSFORM  AND  TIME  RESPO 


The  printer  output  for  this  example  ist 


DEGREE  OF  POLYN  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

4320.  432. 

DEGREE  OF  POLYD  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 

2160.  1008.  345.  35.  1. 

XXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
x  SPLDC  -  LOAD  TRANSFER  FUNCTION  IN  COEFFICIENT  FORM  * 

XXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  NUMERATOR  OF  SPTF11  IS  1  (COEFFICIENTS  IN  ASCENDi'lG  ORDER) 
4320.  432. 

DEGREE  OF  DENOMINATOR  OF  SPTF11  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
2160.  1008.  345.  35.  1. 

BODE  GAIN  =  2.0000000 


CP=  5.58 

THE  NUMERATOR  ROOTS  OF  SR00T11  ARE 
NO.  REAL  IMAG.  OMEGA  ZETA 

1  -10.000000  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  4320.0000 

THE  DENOMINATOR  ROOTS  OF  SR00T11  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 


1 

2 

3 

4 


-1.5000000 

-1.5000000 

-12.000000 

-20.000000 


-2.5980762 

2.5980762 

0. 

0. 


3.00000000 

3.00000000 


LOW  ORDER  NON  ZERO  COEFFICIENT  = 


.500000000 

.500000000 


2160.0000 


DEGREE  OF  NUMERATOR  OF  SPTF11  IS  1 
4320.  432. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


DEGREE  OF  DENOMINATOR  OF  SPTF11  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
2160 .  1008.  345 .  35.  1. 


BODE  GAIN  *  2.0000000 


KXKXXXXXKKXXXKXXXXXXKXXKKXKXXXXXXXKXXXXKXXXXXXXXXXKXXXXXXKXX 

X  STIME  -  TIME  RESPONSE  OF  S-PLANE  TRANSFER  * 

X  FUNCTION  11  X 

XKXXXKKXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXKXXX 


COMPUTE  STEP  RESPONSE  IF  TSTEP.NE.O.,  (TSTEP  =  1.0  ) 

OTHERNISE  COMPUTE  IMPULSE  RESPONSE. 

SCALE  OUTPUT  BY  TMAGN,  (TMAGN  =  1.00000  ) 


NO. 


ROOT 


PARTIAL  FRACTION  COEFFICIENT 


1  -1.5000000 

2  -1.5000000 

3  -12.000000 

4  -20.000000 

5  0. 


-2.5980762 

2.5980762 

0. 

0. 

0. 


-.99977959 
-.99977959 
.76923077E-01 
- . 76923077E-01 
20.000000 


-.69735056 

.69735056 

0. 

0. 

0. 


ANALYTICAL  SOLUTION  IS  THE  SUMMATION  OF  THE  FOLLOWING  4  TERMS 


(  2.0000  )  x  txx  0 

((-2.000  )x( cos(  2.60  xT))+(-1.395  )x(sin(  2.60  XT)))xExx(- 

(  .76923E-01)  x  Exx(  -12.000  x 

(  -.77364E-01 )  X  Exx(  -20.000  x 


1.500 

T) 

T) 


XXX  TIME  RESPONSE  xxx 


TIME 

VALUE 

TIME 

VALUE 

TIME 

VALUE 

0. 

0. 

1.7000 

2.1497 

3.4000 

2.0054 

.50000E-01 

.66945E-02 

1.7500 

2.1235 

3.4500 

2.0066 

.10000 

.44037E-01 

1.8000 

2.0985 

3.5000 

2.0075 

.15000 

.10879 

1.8500 

2.0749 

3.5500 

2.0082 

.20000 

.20673 

1.9000 

2.0530 

3.6000 

2.0086 

.25000 

.32913 

1.9500 

2.0330 

3.6500 

2.0087 

.30000 

.46993 

2.0000 

2.0152 

3.7000 

2.0086 

.35000 

.62345 

2.0500 

1.9995 

3.7500 

2.0084 

.40000 

.78455 

2.1000 

1.9860 

3.8000 

2.0084 

.45000 

.94871 

2.1500 

1.9747 

3.8500 

2.0076 

.50000 

1.1120 

2.2000 

1.9655 

3.9000 

2.0070 

.55000 

1.2711 

2.2500 

1.9582 

3.9500 

2.0063 

.60000 

1.4233 

2.3000 

1.9529 

4.0000 

2.0057 

.65000 

1.5664 

2.3500 

1.9493 

4.0500 

2.0050 

.70000 

1.6985 

2.4000 

1.9472 

4.1000 

2.0042 

.75000 

1.8186 

2.4500 

1.9466 

4.1500 

2.0035 

.80000 

1.9257 

2.5000 

1.9471 

4.2000 

2.0029 

.85000 

2.0194 

2.5500 

1.9487 

4.2500 

2.0022 

.90000 

2.0997 

2.6000 

1.9511 

4.3000 

2.0016 

.95000 

2.1667 

2.6500 

1.9542 

4.3500 

2.0011 

1.0000 

2.2209 

2.7000 

1.9578 

4.4000 

2.0006 

1.0500 

2.2209 

2.7500 

1.9618 

4.4500 

2.0001 

1.1000 

2.2937 

2.8000 

1.9661 

4.5000 

1.9998 

1.1500 

2.3140 

2.8500 

1.9704 

4.5500 

1.9994 

1.2000 

2.3249 

2.9000 

1.9748 

4.6000 

1.9992 

1.2500 

2.3276 

2.9500 

1.9791 

4.6500 

1.9990 

1.3000 

2.3230 

3.0000 

1.9832 

4.7000 

1.9988 

1.3500 

2.3123 

3.0500 

1.9871 

4.7500 

1.9987 

1.4000 

2.2966 

3.1000 

1.9907 

4.8000 

1.9986 

1.4500 

2.2769 

3.1500 

1.9940 

4.8500 

1.9986 

1.5000 

2.2542 

3.2000 

1.9970 

4.9000 

1.9986 

1.5500 

2.2293 

3.2500 

1.9996 

4.9500 

1.9986 

1.6000 

2.2032 

3.3000 

2.0019 

5.0000 

1.9987 

1.6500 

2.1764 

3.3500 

2.0038 

TIME  RESPONSE 

EXAMPLE  6  INVERSE  LAPLACE  TRANSFORM  AND  TIME  RESPONSE 
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3.30 

3.20 

3.10 
3.00 

2.90 
2.80 

2.70 
2.60 

2.50 
2.  AO 

2.30 

2.20 

2.10 
2.00 

1.90 
1.80 

1.70 
1.60 

1.50 
1 .  AO 

1.30 
1.20 
1.10 
1.00 
.900 
.800 
.700 
.600 
.500 
.  AOO 
.300 
.200 

. 100E+00 
0. 

-. 100E+00 
-.200 
-.390 
-.AOO 
-.500 
-.600 
-.700 


xxxxxx 

XX  XX 

X  xxxx 

X  xxxxxx  xxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

X  xxxx 

X 

X 

X 

X 

X 


X 

X 


X 

X 


-.800 

-.900 


* 

I  I  I  1 

I 

I 

I 

I 

I 

I 

I 

0.00  .50  1.00  1.50 

FILM  PLOT  COMPLETED 

2.00 

2.50 

3.00 

3.50 

A. 00 

A. 50 

5.00 

P0LYD<1>*4. 

POL YD<2>=. 243334776 
P0LYDC3)=-1. 58617798 
P0LYD(4)*3. 43416283 
P0LYD<5)=-3. 09127983 
P0LYD(6)S1 . 

CALL  ZPLDC(l) 


wdeg.  of  danom." 

"coeff.  of  z**0" 

"coeff.  of  z**l" 

"coeff.  of  zXX2" 

"coeff.  of  zXX3" 

"coeff.  of  zXXA" 

"load  transfer  function  into  ZPTF11' 


ENTER  TIME  RESPONSE  PARAMETERS  FOR  USE  HITH  ZTIME 


TSTEPsl 


".NE.O  (preset=l)  for  step  response, 

. EQ.O  for  impulse  response” 

TMAGN=1  "magnitude  of  the  input” 

TENDS2.0  "end  time  for  evaluating  the  response” 

SAMPT* . 04  "sampling  period" 

CALL  HEADIN5(1, 

♦50HEXAMPLE  7  INVERSE  Z  TRANSFROM  AND  TIME  RESPONSE  ) 

CALL  HEADIN1(6#10H  ) 

CALL  ZTIMEC1)  "compute  inverse  z  transform  of  ZPTF1" 


DEGREE  OF  DENOMINATOR  OF  ZPTF1  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
.243334776  -1.58617798  3.43416283  -3.09127983  1. 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
X  ZTIME  -  TIME  RESPONSE  OF  Z-PLANE  TRANSFER  * 

X  FUNCTION  1  x 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

COMPUTE  STEP  RESPONSE  IF  TSTEP.NE.O.,  (TSTEP  =  1.0  ) 

OTHERWISE  COMPUTE  IMPULSE  RESPONSE. 

SCALE  OUTPUT  BY  TMAGN,  (TMAGN  =  1.00000  ) 

FINAL  VALUE  =  .999909 

xxx  TIME  RESPONSE  xxx 


TIME 

VALUE 

TIME 

VALUE 

TIME 

VALUE 

. 20135E-03 

.68000 

.13117 

1.3600 

.25310 

. 40000E-01 

. 12321E-02 

.72000 

.14013 

1.4000 

.25856 

. 80000E-01 

. 37411E-02 

.76000 

.14886 

1.4400 

.26388 

.12000 

.78894E-02 

.80000 

.15736 

1.4800 

.26907 

.16000 

. 13486E-01 

.84000 

.16562 

1.5200 

.27413 

.20000 

. 20269E-01 

.88000 

.17364 

1.5600 

.27907 

.24000 

. 27989E-01 

.92000 

.18142 

1.6000 

.28391 

.28000 

. 36424E-01 

.96000 

.18897 

1.6400 

.28864 

.32000 

. 45388E-01 

1.0000 

.19629 

1.6800 

.22927 

.36000 

. 54723E-01 

1.0400 

.20339 

1.7200 

.29781 

.40000 

. 64299E-01 

1.0800 

.21027 

1.7600 

.30227 

.44000 

. 74009E-01 

1.1200 

.21695 

1.8000 

.30664 

.48000 

. 83764E-01 

1.1600 

.22342 

1.8400 

.31093 

.52000 

. 93493E-01 

1.2000 

.22970 

1.8800 

.31516 

.56000 

.10314 

1.2400 

.23581 

1.9200 

.31931 

.60000 

.11266 

1.2800 

.24173 

1.9600 

.32340 

.64000 

.12201 

1.3200 

.24750 

2.0000 

.32743 

AT  THE  LAST  TIME  POINT  THE  SINGLE  PRECISION  VALUE  DIFFERS  FROM  THE 
DOUBLE  PRECISION  VALUE  BY  .16E-08  PERCENT. 
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EXAMPLE  7 

INVERSE  Z  TRANSFORM  AND  TIME 

RESPONSE 

.370 

I 

.360 

I 

.350 

I 

.340 

I 

.330 

I 

X 

.320 

I 

X  XX 

.310 

I 

XX 

.300 

I 

XX 

.290 

I 

XX 

.280 

I 

XX 

.270 

I 

X  X 

.260 

I 

XX 

.250 

I 

XX 

.240 

I 

XX 

.230 

I 

X 

.220 

I 

XX 

.210 

I 

X 

.200 

I 

XX 

.190 

I 

X 

.180 

I 

X 

.170 

I  XX 

.160 

I  X 

.150 

I  X 

.140 

I  X 

.130 

I  X 

.120 

I  X 

.110 

I  X 

. 100E+00 

I  X 

. 900E-01 

I  X 

. 800E-01 

I  X 

. 700E-01 

I  X 

. 600E-01 

I  X 

.500E-01 

I  XX 

. 400E-01 

I  * 

. 300E-01 

I  * 

. 200E-01 

I  X 

. 100E-01 

I  XX 

U  . 

- . 100E-01 

I 

-.200E-01 

I 

- . 300E-01 

I 

-.400E-01 

I 

-.500E-01 

I 

■ 

I  I  I  I  I 

I  I 

I  I 

I  I 

0 

.00  .20  .40  .60  .80 

1.00  1.20 

1.40  1.60 

1.80  2.00 

FILM  PLOT  COMPLETED 

IVlV 


RANSFORM  OF  AN 


Problem)  Compute  the  z  transform  of  the  following  transfer  function  with  a  zero 
order  hold  and  a  delay  of  .008  seconds.  The  sampling  period  is  .08  sec¬ 
onds. 


2900s  +  2900 


s  +  As  +  124s  +  363s 

The  partial  fraction  expansion  method  is  used  to  compute  the  z  transform.  The 
algorithm  used  to  compute  the  partial  fraction  expansion  requires  that  (1)  there 
are  no  multiple  poles  other  than  those  at  the  origin  and  (2)  the  degree  of  the 
numerator  must  not  be  greater  than  the  number  of  non-zero  poles  of  the  denomina¬ 
tor.  Up  to  5  poles  at  the  origin  are  allowed  (this  includes  the  one  from  the  zero 
order  hold  if  applicable). 

The  FORTRAN  code  for  this  example  is> 


DEGREE  OF  POLYN  IS  1 
2900.  2900. 

DEGREE  OF  POLYD  IS  4 
0.  0.  363.  124.  4.  1. 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  SPLDC  -  LOAD  TRANSFER  FUNCTION  IN  COEFFICIENT  FORM  X 

KXXXXXXXXXXXXXXXXXXKXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  NUMERATOR  OF  SPTF12  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 
2900.  2900. 

DEGREE  OF  DENOMINATOR  OF  SPTF12  IS  4  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  0.  363.  124.  4.  1. 

BODE  GAIN  =  7.9889807 


THE  NUMERATOR  ROOTS 

OF  SR00T12  ARE 

NO. 

REAL 

I  MAG. 

OMEGA 

ZETA 

1 

-1.0000000 

0. 

LOW  ORDER  NON  ZERO 

COEFFICIENT  = 

2900.0000 

THE  DENOMINATOR  ROOTS  OF  SR00T12 

ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-.50000000 

-10.988630 

11.0000000 

. 45454545E-01 

2 

-.50000000 

10.988630 

11.0000000 

. 45454545E-01 

3 

-3.0000000 

0. 

4 

0. 

0. 

5 

0. 

0. 

LOW  ORDER  NON  ZERO 

COEFFICIENT  = 

363.00000 

DEGREE  OF  NUMERATOR 

OF  SPTF12  IS  1 

(COEFFICIENTS 

IN 

ASCENDING  ORDER) 

2900. 

.  2900. 

DEGREE  OF  DENOMINATOR  OF  SPTF12  IS  5 

(COEFFICIENTS 

IN 

ASCENDING  ORDER) 

0.  0. 

.  363.  124.  4. 

1. 

BODE 

GAIN  =  7.9889807 

XXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 

X  ZR00T2  =  SZXFM  OF  SR00T12  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXX 

X  ZPTF2  =  SZXFM  OF  SPTF12  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx: 


SAMPLING  PERIOD,  SAMPT  =  .0800  DELAY  =  .0080  ZOH  =  1 

NO.  ROOT  PARTIAL  FRACTION  EXPANSION  COEFFICIENT 


1 

-.50000000 

-10.988630 

.  96684606E-01 

-  .40456452E-02 

2 

-.50000000 

10.988630 

-. 966846 06E-01 

.  40456452E-02 

3 

-3.0000000 

0. 

1.6914552 

0. 

4 

0. 

-1.8848244 

0. 

5 

0. 

0. 

5.2599625 

0. 

6 

0. 

0. 

7.9889807 

0. 

THE 

NUMERATOR  ROOTS 

OF  ZR00T2  ARE 

NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-12.268933 

0. 

2 

- .67135130E 

-04 

0. 

3 

.92311634 

0. 

4 

-.13469859 

0. 

5 

-1.1567519 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  -.23106189E-04 

THE  DENOMINATOR  ROOTS  OF  ZR00T2  ARE 
NO.  REAL  IMAG.  OMEGA  ZETA 

1  .61284138  -.73996066  .96078944  -.63785191 

2  .61284138  .73996006  .96078944  -.63785191 

3  .78662786  0. 

4  1.00000000  0. 

5  1.00000000  0. 

6  0.  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  -46.473538 

DEGREE  OF  NUMERATOR  OF  ZPTF2  IS  5  (COEFFICIENTS  IN  ASCENDING  ORDER) 
-.00002310618894132  -.3443427037346  -2.507846683653  .6794170629569 
2.464719122914  .1950347342772 


1m! 


L 


i 


T 


Ikl 


Problem!  Compute  the  frequency  response  of  the  following  function 


\  T/n  j2*pi*k/n 
/  G  (z  e  ) 


where  T=.24  seconds,  n=3  and  G  (z  )  is  the  z-transform  computed  in 

n 

Example  9. 

The  above  function  is  Sklansky's  frequency  decomposition  method  for  expressing 
the  output  transform  of  a  fast  to  slow  sampler  in  terms  of  the  faster  input 
transfer  function.  An  example  of  how  this  frequency  decomposition  method  can  be 
applied  to  the  stability  analysis  of  a  multiloop  multirate  control  system  is 
given  in  Example  2  of  Ref.  4.  The  above  function  can  represent  the  open  loop 
transfer  function  of  Eq  (6.6)  in  Ref.  4. 

The  operator  ZMRFQ  evaluates  the  frequency  response  of  the  above  function  by 

T/n 

using  only  the  transform  of  the  faster  sampled  signal  G  (z  ).  The  response  is 

n 

computed  by  the  indicated  summation  with  shifted  values  of  z  .  This  operation 

T/n  n 

yields  only  the  multirate  frequency  response  of  G  (z  ).  No  z-transform  at 

n 

the  slower  sampling  rate  is  computed.  In  the  next  example  though,  an  explicit 
form  for  the  implementing  the  frequency  decomposition  method  is  described. 

Since  one-half  of  the  sampling  frequency  of  the  slower  output  sampler  is  13.09 
rad/sec,  the  frequency  range  to  be  used  for  this  example  will  be  1.0  to  13.0 
rad/sec . 

The  FORTRAN  code  for  this  example  isi 


EXAMPLE  9 

TRANSFER  FUNCTION  IS  ZPTF2  FROM  PREVIOUS  CASE 
ENTER  FREQUENCY  RESPONSE  PARAMETERS  FOR  ZMRFQ 


N0MEG=2 
OMEGA ( 1 )  =1 . 
0MEGA(2)=13. 
SAMPT=.24 


"number  of  values  of  OMEGA  to  be  entered" 
"0MEGA( 1 ) =first  frequency  value  to  be  used" 
"0MEGA(N0MEG)slast  frequency  value  to  ba  used" 
"sampling  period  of  slower  output  sampler" 


pit 


CALL  REMARK5( 

♦50HZPTF2  IS  AT  THE  FASTER  SAMPLING  RATE,  SAMPT=.08  ) 

INTEGER  RATIO  OF  OUTPUT/INPUT  SAMPLING  PERIODS  IS  ENTERED  AS  THE 
SECOND  ARGUMENT  OF  ZMRFQ(I,M) 

CALL  REMARK5( 

+50HFREQ.  RESP.  MILL  BE  AT  THE  SLOWER  SAMPLING  RATE,  ) 

CALL  REMARK2<  20HSAMPT® . 08  ) 

CALL  HEADIN5Q, 

+50HEXAMPLE  9  MULTIRATE  FREQUENCY  RESPONSE  BY  FREQUENC) 

CALL  HEADIN2(6,20HY  DECOMPOSITION  ) 

CALL  ZMRFQ(2,3)  "compute  multirate  frequency  response  of  ZPTF2 


C 


The  printer  output  for  this  example  isi 


ZPTF2  IS  AT  THE  FASTER  SAMPLING  RATE,  SAMPT=.08 
FREQ.  RESP.  HILL  BE  AT  THE  SLOWER  SAMPLING  RATE, 
SAMPT  * . 08 


THE  NUMERATOR  ROOTS  OF  ZR00T2  ARE 


NO.  REAL  IMAG.  OMEGA  ZETA 

1  -12.268933  0. 

2  -.67135130E-04  0. 

3  .92311634  0. 

4  -.13469859  0. 

5  -1.1567519  0. 


LOW  ORDER  NON  ZERO  COEFFICIENT  =  - . 23106189E-04 

THE  DENOMINATOR  ROOTS  OF  ZR00T2  ARE 


NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

.61284138 

-.73996066 

.96078944 

-.63785191 

2 

.61284138 

.73996006 

.96078944 

-.63785191 

3 

.78662786 

0. 

4 

1.00000000 

0. 

5 

1.00000000 

0. 

6 

0. 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  -46.473538 


DEGREE  OF  NUMERATOR  OF  ZP7F2  IS  5  (COEFFICIENTS  IN  ASCENDING  ORDER) 

- . 00002310618894132  - . 3443427057346  -2.507846683653  .6794170629569 
2.464719122914  .1950347342772 

DEGREE  OF  DENOMINATOR  OF  ZPTF2  IS  6  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  -46.47353837271  213.7325205421  -416.8323061287  442.3612041218 
-256.7878801625  64. 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  ZMRFQ  -  MULTIRATE  FREQUENCY  RESPONSE  (BY  X 

x  FREQUENCY  DECOMPOSITION)  OF  ZPTF2  * 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SAMPLING  PERIOD  ( SAMPT )  =  .2400 

INTEGER  RATIO  OF  OUTPUT/INPUT  SAMPLING  PERIODS,  MTGER  =  3 

AUTOMATIC  FREQUENCY  MODE  IF  FAUTO.NE.O,  FAUTO  =  1.000 
NOMEG  =  2.000  OMEGA  =1.000  ,  13.00 


OMEGA 

ZREAL 

ZIMAG 

REAL 

IMAGINARY 

DB 

PHASE 

PHASE 

RAD/SEC 

MARGIN 

1.000 

.997 

.080 

- . 331E+01 

-. 143E+01 

11.125 

-156.66 

23.34 

1.200 

.995 

.096 

- . 247E+01 

- . 113E+01 

8.664 

-155.48 

24.52 

1.400 

.994 

.112 

- .  195E+01 

- . 906E+00 

6.648 

-155.08 

24.92 

1.700 

.991 

.136 

- . 147E+01 

-.670E+00 

4.171 

-155.51 

24.49 

2.000 

.987 

.159 

- .  117E+01 

- . 504E+00 

2.125 

-156.75 

23.25 

2.369 

.982 

.188 

-.933E+00 

-.360E+00 

-.000 

-158.91 

21.09 

2.869 

.974 

.220 

-.722E+00 

-.230E+00 

-2.411 

-162.36 

17.64 

3.469 

.962 

.274 

- . 560E+00 

- . 132E+00 

-4.805 

-166.75 

13.25 

4.269 

.942 

.335 

-.423E+00 

- . 556E-01 

-7.377 

-172.53 

7.47 

5.069 

.919 

.394 

-.  339E+00 

-. 112E-01 

-9.387 

-178.11 

1.89 

5.346 

.910 

.415 

- . 317E+00 

.458E-07 

-9.966 

-180.00 

-.00 

6.946 

.850 

.528 

-.241E+00 

. 467E-01 

-12.206 

-190.98 

-10.98 

8.546 

.775 

.632 

- .  225E+00 

. 989E-01 

-12.202 

-203.76 

-23.76 

9.546 

.722 

.692 

-.250E+00 

. 181E+00 

-10.213 

-215.89 

-35.39 

10.05 

.694 

.720 

-.270E+00 

. 286E+00 

-8.096 

-226.62 

-46.62 

10.35 

.677 

.736 

- . 262E+00 

. 409E+00 

-6.287 

-237.50 

-57.50 

10.55 

.655 

.747 

- .  208E+00 

. 530E+00 

-4.893 

-248.60 

-68.60 

10.70 

.656 

.755 

- .  110E+00 

. 630E+00 

-3.887 

-260.08 

-80.08 

10.82 

.648 

.762 

.  250E-01 

. 688E+00 

-3.247 

-272.08 

-92.08 

10.92 

.642 

.767 

.  160E+00 

. 690E+00 

-2.994 

-283.03 

-103.03 

11.02 

.636 

.772 

.  292E+00 

. 640E+00 

-3.054 

-294.53 

-114.53 

11.12 

.630 

.777 

. 393E+00 

. 546E+00 

-3.442 

,-305.70 

-125.70 

11.22 

.623 

.782 

. 447E+00 

. 435E+00 

-4.097 

-315.81 

-135.81 

11.37 

.614 

.789 

.  458E+00 

. 283E+00 

-5.376 

-328.22 

-148.22 

11.57 

.601 

.799 

.  409E+00 

. 150E+00 

-7.222 

-339.91 

-159.91 

335E+00 

.659E-01 

-9.326 

-398.97 

-168.97 

270E+00 

. 236E-01 

-11.393 

-355.00 

-175.00 

212E+00 

. 339E-02 

-13.960 

-359.09 

-179.09 

198E+00 

.  333E-03 

-19.077 

-359.90 

-179.90 

The  z-plane  frequencies  in  tha  tabulation  of  the  response  are  at  the  faster  sam¬ 
pling  rata.  Note  that  at  13.0  rad/sac  tha  z-plane  frequency  is  (.506  +  j.852)  or 
approximately  60  degrees  on  the  unit  circle.  At  the  slower  sampling  rate  this 
point  would  be  at  approximately  the  (-1.  +  j  0.)  point,  which  is  3  times  60 
degrees  on  the  unit  circle. 


Problem  Compute  the  frequency  response  of  the  function 

n~l 
1  \ 

\  T/n  j2*pi*k/n 
n  /  G  (z  e  ) 


using  the  explicit  form  of  Sklansky's  frequency  decomposition  method 

T/n 

where  T= .  24  seconds *  n= 3  and  G  (z  )  is  the  z-transform  computed  in 

n 

Example  9. 

In  Example  9  the  frequency  response  was  computed  using  the  operator  ZMRFQ 
which  numerically  evaluated  the  function  over  a  range  of  frequencies  for  z  . 

n 

Until  recently*  this  wss  the  only  means  in  LCAP2  of  applying  the  frequency 
decomposition  method.  A  new  operator*  ZMRXFM,  has  been  implemented  which  will 
yield  a  rational  representation  of  this  frequency  decomposition  method.  This 
operator  computes  the  output  transform  of  a  fast  to  slow  sampler.  This  output 
transform  is  a  transfer  function  represented  in  rational  form  at  the  slower  sam¬ 
pling  rate.  After  this  transfer  function  has  been  computed*  the  single  rate 
z-plane  frequency  response  operator*  ZFRE9*  can  be  used  to  evaluate  the 
response . 

The  numericel  technique  for  computing  the  rational  form  of  the  frequency  decom¬ 
position  method  will  be  documented  in  the  near  future. 

The  FORTRAN  code  for  this  example  is> 


EXAMPLE  10 

G(Z)  IS  ZPTF2  FROM  THE  PREVIOUS  EXAMPLE 
ENTER  PARAMETERS  FOR  ZMRXFM 


SAMPT=.24  "Sampling  period  of  slower  output  sampler" 

NTGER=3  "integer  ratio  of  output/input  sampling  periods" 

CALL  REMARK5C 

♦50ZPTF2  IS  AT  THE  FASTER  SAMPLING  RATE,  SAMPT=.0S 
CALL  REMARK5( 

+50HZPTF3  HILL  BE  AT  THE  SLOWER  SAMPLING  RATE* 

CALL  REMARK2(20HSAMPT= . 24  ) 


c 

c 

c 


CALL  ZMRXFM(3,2)  "compute  multirate  transform  of  ZPTF2  and  store 

in  ZPTF3" 


ZPTF3  IS  AT  THE  SLOWER  SAMPLING  RATE 


C  ENTER  FREQUENCY  RESPONSE  PARAMETERS  FOR  USE  WITH  ZFREQ 
C 

N0MEG=2  "number  of  values  of  OMEGA  to  be  entered" 

0MEGA(1)=1.  "OMEGA( 1 )=first  frequency  value  to  be  used" 

0MEGA(2)=13.  "0MEGA( NOMEG) =last  frequency  value  to  be  used" 

CALL  HEADIN5C1, 

♦50HEXAMPLE  10  RATIONAL  REPRESENTATION  OF  FREQUENCY  DE) 


CALL  HEADIN2C6#  20HCOMPOSITION  METHOD  ) 

CALL  ZFREQ(3)  "compute  frequency  response 

of  ZPTF3" 

The  printer  output  for 

this  example  is: 

ZPTF2 

IS  AT  THE  FASTER 

SAMPLING  RATE,  SAMPT=.08 

ZPTF3 

SAMPT  * 

HILL  BE  AT  THE  SLOWER  SAMPLING  RATE# 

1 . 24 

THE  NUMERATOR  ROOTS  OF  ZR00T2  ARE 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

-12.268933 

0. 

2 

-.67135I30E-04 

0. 

3 

.92311634 

0. 

4 

-.13469859 

0. 

5 

-1.1567519 

0. 

LOW 

THE 

ORDER  NON  ZERO  COEFFICIENT  = 

DENOMINATOR  ROOTS  OF  ZROOT2  ARE 

■ . 23106189E-04 

NO. 

REAL 

IMAG.  OMEGA 

ZETA 

1 

.61284138 

-.73996066  .96078944 

-.63785191 

2 

.61284138 

.73996006  .96078944 

-.63785191 

3 

.78662786 

0. 

4 

1 .00000000 

0. 

5 

1.00000000 

0. 

6 

0. 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  = 


-46.473538 


DEGREE  OF  NUMERATOR  OF  ZPTF2  IS  5  (COEFFICIENTS  IN  ASCENDING  ORDER) 
-.00002310618894132  -.3443427037346  -2.507846683653  .6794170629569 
2.464719122914  .1950347342772 

DEGREE  OF  DENOMINATOR  OF  ZPTF2  IS  6  (COEFFICIENTS  IN  ASCENDING  ORDER) 
0.  -46.47353837271  213.7325205421  -416.8323061287  442.3612041218 
-256.7878801625  64. 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
*  ZR00T3  =  ZMRXFM  OF  ZR00T2  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
x  ZPTF3  =PSYNTH(ZR00T3  )  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


THE  NUMERATOR  ROOTS  OF  ZR00T3  ARE 


NO.  REAL  IMAG.  OMEGA  ZETA 

1  - . 10647750E-01  0. 

2  -.59161620  0. 

3  .78658382  0. 

4  -2.9655282  0. 


LOW  ORDER  NON  ZERO  COEFFICIENT  =  -.18162533 

THE  DENOMINATOR  ROOTS  OF  ZR00T3  ARE 


NO. 

REAL 

IMAG. 

OMEGA 

ZETA 

1 

-.77650116 

-.42857182 

.88692044 

.87550262 

2 

-.77650116 

.42857182 

.88692044 

.87550262 

3 

.48675226 

0. 

4 

1.00000000 

0. 

5 

1.00000000 

0. 

LOW  ORDER  NON  ZERO 

COEFFICIENT  = 

-24.505145 

DEGREE 

OF  NUMERATOR 

OF  ZPTF3  IS  4 

(COEFFICIENTS 

IN 

ASCENDING  ORDER) 

-.1816253320627  -17 

.19496629374  -12.53387821671  34.37676332934 

12.3603697263 

DEGREE 

OF  DENOMINATOR  OF  ZPTF3  IS  5 

(COEFFICIENTS 

IN 

ASCENDING  ORDER) 

-24.5051447024  50.97511989617  39.80519865903  -70.51517819695 
-59.75999565584  64. 


KXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXX 
*  ZFREQ  -  FREQUENCY  RESPONSE  OF  Z-PLANE  TRANSFER  X 

x  FUNCTION  3  x 

KXXXKXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

SAMPLING  PERIOD  (SAMPT)  =  .2400 

AUTOMATIC  FREQUENCY  MODE  IF  FAUTO.NE.O,  FAUTO  =  1.000 
NOMEG  *  2.000  OMEGA  <  1.000  ,  13.00 


OMEGA 

RAD/SEC 

ZREAL 

ZIMAG 

REAL 

IMAGINARY 

DB 

PHASE 

PHASE 

MARGIN 

1.000 

.971 

.238 

- . 331E+01 

- . 143E+01 

11.125 

-156.66 

23.34 

1.200 

.959 

.284 

- . 247E+01 

-.113E+01 

8.664 

-155. 48T 

24.52 

1.400 

.944 

.330 

- . 195E+01 

- . 906E+00 

6.648 

-155.08 

24.92 

1.700 

.918 

.397 

-. 147E+01 

-.670E+00 

4.171 

-155.51 

24.49 

2.000 

.887 

.462 

-.  117E+01 

- . 504E+00 

2.125 

-156.75 

23.25 

2.369 

.843 

.538 

- .  933E+00 

-.360E+00 

-.000 

-158.91 

21.09 

2.869 

.772 

.635 

- .  722E+00 

- . 230E+00 

-2.411 

-162.36 

17.64 

3.469 

.673 

.740 

-.560E+00 

-. 132E+00 

-4.805 

-166.75 

13.25 

4.269 

.519 

.854 

-.423E+00 

- . 556E-01 

-7.377 

-172.53 

7.47 

5.069 

.347 

.938 

-.339E+00 

- . 112E-01 

-9.387 

-178.11 

1.89 

5.346 

.284 

.959 

- . 317E+00 

.458E-07 

-9.966 

-180.00 

-.00 

6.946 

-.096 

.995 

- . 241E+00 

. 467E-01 

-12.206 

-190.98 

-10.98 

8.546 

-.462 

.887 

- . 225E+00 

.989E-01 

-12.202 

-203.76 

-23.76 

9.546 

-.660 

.752 

- . 250E+00 

. 181E+00 

-10.213 

-215.89 

-35.39 

10.05 

-.745 

.667 

-.270E+00 

. 286E+00 

-8.096 

-226.62 

-46.62 

10.35 

-.791 

.612 

-.262E+00 

. 409E+00 

-6 . 287 

-237.50 

-57.50 

10.55 

-.819 

.573 

-.208E+00 

. 530E+00 

-4.893 

-248.60 

-68.60 

10.70 

-.839 

.544 

-.  110E+00 

. 630E+00 

-3.887 

-260.08 

-80.08 

10.82 

-.855 

.518 

. 250E-01 

.  688E+00 

-3.247 

-272.08 

-92.08 

10.92 

-.868 

.497 

.  160E+00 

.690E+00 

-2.994 

-283.03 

-103.03 

11.02 

-.879 

.476 

.292E+00 

•640E+00 

-3.054 

-294.53 

-114.53 

11.12 

-.890 

.455 

. 393E+00 

.546E+00 

-3.442 

-305.70 

-125.70 

11.22 

-.901 

.434 

•447E+00 

. 435E+00 

-4.097 

-315.81 

-135.81 

11.37 

-.916 

.401 

.  458E+00 

. 283E+00 

-5.376 

-328.22 

-148.22 

11.57 

-.934 

.357 

. 409E+00 

. 150E+00 

-7.222 

-339.91 

-159.91 

11.82 

-.954 

.300 

.  335E+00 

•654E-01 

-9.326 

-348.97 

-168.97 

12.12 

-.973 

.231 

. 270E+00 

. 236E-01 

-11.343 

-355.00 

-175.00 

12.62 

-.994 

.112 

. 212E+00 

. 339E-02 

-13.460 

-359.09 

-179.09 

13.00 

-1.000 

.022 

.  198E+00 

. 333E-03 

-14.077 

-359.90 

-179.90 

'l. 
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NICHOLS  PLOT  (MAGN.  VS  PHASE)  10/31/83 

EXAMPLE  10  RATIONAL  REPRESENTATION  OF  FREQUENCY  DECOMPOSITION  METHOD 


xxxxxxx 

X 

XXX  XX 

X 

XX  XX 

X 

XX 

X 

XX 

XX 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

XX 

X 

xxxxx 

X 

X 

X 

X 

I  I  I  I 

I 

I  I  I  ,  I 

360.000  -288.000 

FILM  PLOT  COMPLETED 

-216.000 

-144.000  -72.000 

I 


I 


IJ>J 


Problem i  Given 


s  +  2 


-(  s  +  3  ) 


M(s)  = 


.01  s  +  .15  s  +1 


fits)  = 


0 

x  (s) 

1 

and  X(s>  = 

1 

x  (s) 

find  the  transfer  function  x  (s)/u(s)  by  application  of  Cramer's  method. 

1 


H(s)  = 


x  (s)  det  H  (s) 
1  1 


det  M(s) 


where  M  (s)  is  equal  to  {j(s)  with  column  1  replaced  by  £(s) . 

1 

Three  operations  are  required  for  this  example.  The  first  two  compute  the  deter 
minants  and  store  the  results  in  polynomials  and  the  last  one  copies  the  polyno 
mials  into  an  s-plane  transfer  function. 

The  FORTRAN  code  for  this  example  isi 


EXAMPLE  11 

MUST  INITIALIZE  MATRIX  PARAMETERS  FOR  FIRST  TIME  USE 
CALL  MINITO 

NOTE  -  DO  NOT  CALL  MINITO  AGAIN 


ENTER  PARAMETERS  FOR  DIMENSION  AND  DEGREE  OF  MATRIX 


MXM=2  "dimension  of  matrices" 

MDEGS2  "highest  degree  of  polynomial  element" 

ENTER  MATRIX  DATA  FOR  COEFFICIENTS  OF  S**0 

M0(l,l)=2 


M0(l,2)=-3 

M0(2,2)=l 

ENTER  MATRIX  DATA  FOR  COEFFICIENTS  OF  S**l 
Ml(l,l)=l 
Ml(  1 , 2)=-l 
Ml (2,2)=. 15 

ENTER  MATRIX  DATA  FOR  COEFFICIENTS  OF  S**2 
M2 (2,2)=. 01 
ENTER  B  VECTOR 
B0(2)=l . 

COMPUTE  DETERMINANT  OF  THE  DENOMINATOR  OF  H(S) 

CALL  REMARK5( 

+50HDETERMINANT  OF  DENOMINATOR  MILL  BE  IN  POLY2  AND 
CALL  REMARK1 ( 10HROOT2  ) 

CALL  DTERM(2,0)  "compute  determinant  of  M(s)" 

"if  2nd  argument  is  zero,  no 
column  substitution  is  made) 


COMPUTE  DETERMINANT  OF  THE  NUMERATOR  OF  H(S) 

CALL  REMARK5( 

+50HDETRMINANT  OF  NUMERATOR  WILL  BE  IN  P0LY3  AND 
CALL  REMARKK10HROOT3  ) 

CALL  DTERM(3, 1 )  "compute  determinant  of  M  (s)" 

1 

"2nd  argument  is  column  number 
where  B(s)  is  substituted" 


COPY  POLYNOMIALS  INTO  S  PLANE  TRANSFER  FUNCTION 


CALL  CPYPS(I3,3,2) 


The  printer  output  for  this  example  is> 


DETERMINANT  OF  DENOMINATOR  MILL  BE  IN  P0LY2  AND 


R00T2 


MATRIX  M(S)  is 


ROM  COL 


1  1  . 20000000E+01  . 1 0000000E+01 

1  2  - . 30000000E+01  -. 10000000E+01 

2  2  . 10000000E+01  . 15000000E+00  . 1 0000000E-01 


V  v  v 

>w 


•v>; 

•  *  » 

1-3. 


•  v  .•* 

:v.-> 


B  VECTOR  IS 


12  3 

S  S  S 

2  . 10000000E+01 


ROW 


XXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  DETRM  -  FIND  DETERMINANT  OF  MATRIX  X 
X  WITH  NO  COLUMN  REPLACED  BY  B  VECTOR  X 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


XXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  R00T2  =  ROOTS  OF  DETERMINANT  X 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXX 

THE  ROOTS  OF  R00T2  ARE 


NO.  REAL  I MAG. 


OMEGA 


ZETA 


1  -7.5000000 

2  -7.5000000 

3  -2.0000000 


6 .61437828 
-6.61437828 
0. 


10.0000000 

10.0000000 


.750000000 

750000000 


LOW  ORDER  NON  ZERO  COEFFICIENT  -  2.0000000 


XKXXXXXXXXXXXXXXXXKXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  P0LY2  =  COEFFICIENTS  OF  DETERMINANT  POLYNOMIAL  x 

XXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


DEGREE  OF  P0LY2  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 

2.  1.3  .17  .01 


DETERMINANT 

R00T3 

MATRIX  M(S) 

OF  NUMERATOR  WILL 

is 

0 

BE  IN  POLY3  AND 

1 

2 

ROW  COL 

1  1 

S 

.20000000 E+ 01 

S 

.  10000000E+01 

S 

1  2 

-.30000000E+01 

- .  1 0000000E+01 

t 

2  2 

.  10000000E+01 

.  15000000E+00 

. 10000000E-01 

B  VECTOR  IS 


0  12  3 

RON  S  S  S  S 

2  . 10000000E+01 


XXXXXXKXXXXXXXKXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  DETRM  -  FIND  DETERMINANT  OF  MATRIX  X 

X  WITH  COLUMN  1  REPLACED  BY  B  VECTOR  * 

XKXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

XXXXXXXKKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  R00T3  =  ROOTS  OF  DETERMINANT  X 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THE  ROOTS  OF  ROOT 3  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -3.0000000  0. 

LON  ORDER  NON  ZERO  COEFFICIENT  =  000000 

XXXXXXXXXXXXXXXXXXXXXXKXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  P0LY3  =  COEFFICIENTS  OF  DETERMINANT  OF  POLYNOMIAL  x 

XXXXXXXXXXXXXXXXXXXXKXXXXXXKKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

DEGREE  OF  P0LY3  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

3.  1. 


THE  ROOTS  OF  R00T3  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -3.0000000  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  3.0000000 


DEGREE  OF  POLY3  IS  1 


(COEFFICIENTS  IN  ASCENDING  ORDER) 


THE  ROOTS  OF  R00T2  ARE 


NO.  REAL  IMAG.  OMEGA  ZETA 

1  -7.5000000  6 .61437828  10.0000000  .750000000 

2  -7.5000000  -6.61437828  10.0000000  750000000 

3  -2.0000000  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  2.0000000 

DEGREE  OF  P0LY2  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 

2.  1.3  .17  .01 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
X  SR00T13  =  R00T3  /  R00T2  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
x  SPTF13  =  P0LY3  /  P0LY2  X 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

THE  NUMERATOR  ROOTS  OF  SR00T13  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -3.0000000  0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  3.0000000 

THE  DENOMINATOR  ROOTS  OF  SR00T13  ARE 


REAL 

IMAG. 

OMEGA 

ZETA 

-7.5000000 

6.61437828 

10.0000000 

.750000000 

-7.5000000 

-6.61437828 

10.0000000 

750000000 

-2.0000000 

0. 

LOW  ORDER  NON  ZERO  COEFFICIENT  =  2.0000000 

DEGREE  OF  NUMERATOR  OF  SPTF13  IS  1  (COEFFICIENTS  IN  ASCENDING  ORDER) 

3.  1. 

DEGREE  OF  DENOMINATOR  OF  SPTF13  is  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 
2.  1.3  .17  .01 

BODE  GAIN  =  1.5000000 


Problem  Store  all  polynomial*  transfer  function,  and  matrix  data  so  that  it  can 
bo  used  in  a  subsequent  batch  or  interactive  LCAP2  job.  Assume  that 
Examples  1  through  11  were  all  run  on  the  same  batch  job. 

The  LCAP2  STORE  operator  will  store  all  pertinent  data  from  a  current  job  so 
that  on  subsequent  batch  or  interactive  jobs  this  data  can  be  restored  with  the 
RESTORE  operator.  This  data  is  written  on  TAPE31.  At  the  completion  of  the 
batch  job  the  user  must  cetalog  this  file. 

To  identify  the  data  on  TAPES1,  the  first  record  will  contain  70  characters  of 
alphanumeric  data  copied  from  words  HEADC64)  through  HEAD(70).  The  user  must 
enter  this  data  before  calling  STORE. 

The  FORTRAN  code  for  this  example  isi 


C  EXAMPLE  12 

C  ENTER  INFORMATION  TO  IDENTIFY  THE  TAPE  BEING  WRITTEN  ON 
CALL  HEADIN5(64, 

♦50HEXAMPLES  FOR  BATCH  LCAP2  USERS  MANUAL  ) 

CALL  HEADIN2(69, 20H0CT0BER  31,1983  ) 

C 

C  THE  ABOVE  CODE  IS  THE  SAME  AS 

C  HEAD(64)=10HEXAMPLES  F 

C  HEAD(65)=10HOR  BATCH  L 

C 
C 

C  HEADC69) =10HOCTOBER  31 

C  HEAD(70)=10H, 1983 

C 


CALL  STORE(l) 
CALL  LEXIT 
END 


if  argument  is  0  printout  will  be  suppressed” 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  STORE  -  STORE  POLYNOMIALS,  TRANSFER  FUNCTIONS,  X 
X  AND  MATRIX  DATA  FOR  A  FUTURE  JOB  X 
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


TAPE31  IDENTIFIER  IS 

EXAMPLES  FOR  BATCH  LCAP2  USERS  MANUAL  OCTOBER  31,1983 


THE  ROOTS 

OF 

R00T1  ARE 

NO. 

REAL 

I  MAG. 

OMEGA 

ZETA 

1 

-1.8445105 

-.49938380 

1.9109168 

.96524896 

2 

-1.8445105 

.49938380 

1.9109168 

.96524896 

3 

-9.3109790 

0. 

DEGREE 

OF  P0LY1  IS 

3 

(COEFFICIENTS  IN  ASCENDING  ORDER) 

34.  38. 

.  13.  1. 

THE  ROOTS 

OF 

R00T2  ARE 

NO. 

REAL 

I  MAG. 

OMEGA 

ZETA 

1 

-7.5000000 

6.6143783 

10.000000 

.75000000 

2 

-7.5000000 

-6.6143783 

10.000000 

.75000000 

3 

-2.0000000 

0. 

DEGREE 

OF  P0LY2  IS 

3 

(COEFFICIENTS  IN  ASCENDING  ORDER) 

2.  1.3 

.17  .01 

(data  for  polynomial  3) 


(data  for  s-plane  transfer  functions  1,...,5) 


(data  for  z-plane  transfer  functions  1,2,3) 


(data  for  s-plane  transfer  functions  6,..., 13) 


MXM  =  2 
MDEG  =  2 


MATRIX  M(S)  is 


0  12  3 

ROM  COL  S  S  S  S 

1  1  . 20000000E+01  . 10000000E+01 

1  2  - . 30000000E+01  - . 10000000E+01 

2  2  . 10000000E+01  . 15000000E+00  . 10000000E-01 

B  VECTOR  IS 

0  12  3 

ROM  S  S  S  S 

2  . 10000000E+01 


Note  that  the  s-plane  transfer  functions  are  not  listed  continuously  from  1 
through  13.  This  is  due  to  the  method  in  which  polynomials  and  transfer  func¬ 
tions  are  stored  in  the  program.  All  polynomials  and  transfer  functions  with 
indices  1  through  5  are  saved  in  SCM  (small  core  memory).  All  polynomials  and 
transfer  functions  with  indices  greater  than  5  are  saved  on  disk  storage.  In 
printing  out  the  data  stored  by  the  STORE  operator,  the  data  which  is  saved  in 
SCM  is  printed  out  before  the  data  which  is  saved  on  disk  storage. 

After  execution  of  Batch  LCAP2  has  been  completed,  the  job  will  proceed 
according  to  the  operations  specified  in  the  control  card  input  file  in  Section 
5.0.  First,  file  TAPE31  will  be  cataloged.  For  this  example  the  file  name  used 
is  8LCAP2BATCHDATA  with  ID=9487.  This  file  will  be  used  by  Example  13.  Then,  the 
hardcopy  plots  will  be  processed  by  the  program  HARDCPY. 


Problems  Restore  data  stored  in  Example  12 

In  Example  12,  the  STORE  operator  was  used  to  save  data  on  file  TAPE31.  Upon  com¬ 
pletion  of  the  batch  job  this  file  was  stored  as  file  8LCAP2BATCKDATA,  ID=9487. 
To  restore  this  data  for  a  new  batch  job  this  file  must  first  be  attached  as  file 
TAPE30  before  LCAP2  is  executed. 

The  job  structure  for  this  example  iss 


(control  cards  for  accounting) 

FILE,TAPE30,BT=I. 

ATTACH! TAPE30 , 8 LCAP2BATCHDATA, ID=9487 , ST=PF6 ) 
ATTACH(X, 8LCAP2CC, ID=9487 ) 

BEGIN, LCAP2CC,X. 

*EOR 

XIDENT  idname 
XDECK  MAIN 
XCALL  LCAP2 
C 

CALL  INITO 
CALL  MINITO 


(user's  FORTRAN  code) 


CALL  LEXIT 


The  user's  FORTRAN  code  for  this  example  iss 


RESTORE  DATA  STORED  IN  EXAMPLE  12 


CALL  RESTORE! 1 ) 


"if  argument  is  0  printout  will  be  suppressed" 


KXXXXXKXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

X  RESTORE  -  RESTORE  POLYNOMIALS,  TRANSFER  FUNCTIONS,  x 
X  AND  MATRIX  DATA  FROM  A  PREVIOUS  JOB  X 

XXXXXXXKXKXXXXXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

TAPE30  IDENTIFIER  IS 

EXAMPLES  FOR  BATCH  LCAP2  USERS  MANUAL  OCTOBER  31,1983 

THE  ROOTS  OF  ROOT1  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -1 . 8445105  -.49938380  1.9109168  .96524896 

2  -1.8445105  .49938380  1.9109168  .96524896 

3  -9.3109790  0. 

DEGREE  OF  POLY1  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 

34.  38.  13.  1. 

THE  ROOTS  OF  R00T2  ARE 

NO.  REAL  IMAG.  OMEGA  ZETA 

1  -7.5000000  6.6143783  10.000000  .75000000 

2  -7.5000000  -6.6143783  10.000000  .75000000 

3  -2.0000000  0. 

DEGREE  OF  P0LY2  IS  3  (COEFFICIENTS  IN  ASCENDING  ORDER) 

2.  1.3  .17  .01 

(data  for  polynomial  3) 


(data  for  s-plane  transfer  functions  1,...,5) 


(data  for  z-plane  transfer  functions  1,2,3) 


(data  for  s-plane  transfer  functions  6,..., 13) 


1  Since  the  batch  job  represented  by  Examples  1  through  11  generated  nineteen 
different  data  sets,  for  brevity  only  a  representative  sample  of  the  output 
for  this  example  is  given  here. 


MATRIX  M(S)  is 


RON  COL 


0 


S 


12  3 

S  S  S 


1  1 
1  2 
2  2 


. 20000000E+01 
-.30000000E+01 
.i0000000E+01 


. 10000000E+01 
10000000E+01 
. 15000000E+00 


. 1 0000000E-01 


B  VECTOR  IS 


ROM 


0 


s 


1 


s 


2 


S 


3 


S 


2  . 10000000E+01 


9.0  MORE  EXAMPLES 


In  Section  8  examples  were  presented  for  some  commonly  used  LCAP2 
operators.  The  FORTRAN  code  and  the  printer  output  for  each  example  were  pre¬ 
sented.  In  this  section  more  examples  covering  advanced  operations  will  be  pre¬ 
sented.  However,  no  printer  output  will  be  included. 

Two  examples  are  included.  The  first  one  describe  the  use  of  Cramer's  meth¬ 
od  for  transfer  function  evaluation.  The  second  one  describes  the  use  of  a 
user-supplied  function  to  compute  the  s-plane  frequency  response  of  a 
non-rational  function.  This  technique  can  also  be  utilized  for  rational  func¬ 
tions  which  exceed  LCAP2's  restriction  that  transfer  functions  be  less  than 
degree  50.  The  frequency  response  of  two  transfer  functions,  each  whose  degree 
is  less  that  50,  but  whose  sum  is  greater  than  50,  can  be  computed  in  a  similar 
manner  as  that  described  in  Example  15. 


Problem  Compute  the  closed  loop  transfer  function  y/u  and  the  open  loop  trans¬ 
fer  function  v/u  of  the  following  system  using  Cramer's  method. 


For  stability  analysis  the  loop  is  to  be  broken  after  the  summation  of  the 
two  feedback  loops.  Defind  a  flag  KFIG»  whose  value  is  to  be  0  or  1>  which  will 
be  used  to  open  or  close  this  feedback  loop.  Redraw  the  above  figure  with  his 
flag  and  with  state  variables  xl»  through  x6  to  obtain 


and  1(a)  is  given  by 


T 

£(s)  -  (A, 0,0, 0,0,0) 


Tha  closed  loop  transfer  function  y/u  is  given  by 


y(s)  x  (s)  det  M  (s) 

3  3 


u(s)  u(s)  det  £(s> 

where  M  (s)  is  equal  to  {](s)  with  column  1  replaced  by  £(s)  and  the  value 
3 

of  KFLG  =  1. 


The  open  loop  transfer  function  v/u  is  given  by 


v(s)  x  (s)  det  M  (s) 

6  6 


uCs)  u(s)  det  {]Cs) 

where  H  (s)  is  equal  to  {](s)  with  column  6  replaced  by  J)(s)  and  the  value 
6 

of  KFLG  *  0. 

To  use  LCAP2  to  evaluate  the  above  transfer  functions*  the  polynomial 
matrix  M(s)  is  represented  by  the  input  matrices  MU*  Ml,  M2*  Ml*  and  M4  as 

A  3  2  1 

M(s>  =  M  s  +  Mi  s  +  M2s  +  Ml  s  +  MSI 

The  Mi*  i  =  l*  •  •*«  input  matrices  should  not  be  confused  with  the  M(s)  matrix 

i 

used  in  computing  the  numerator  of  a  transfer  function. 


The  FORTRAN  code  for  this  example  isi 


REAL  K1,K2,KFLG 

MUST  INITIALIZE  MATRIX  PARAMETERS  FOR  FIRST  TIME  USE 

CALL  INITO 
CALL  MINITO 

NOTE  -  DO  NOT  CALL  MINITO  AGAIN 

ENTER  PARAMETERS  FOR  ORDER  AND  DEGREE  OF  MATRIX 

MXM=6 
MDEG=2 

ENTER  DATA 

A= . 

B= . 

C= . 

D= . 

E= . 

F= . 

G= . 

H= . 

Kl= . 

K2= . 

ENTER  DATA  FOR  CLOSED  LOOP  CASE  FIRST 
M0( 1 , 1 ) =C 

M0C1,6)=KFLG=1.  "set  KFLG  for  closed  loop" 

M0(2,l)=-E 

M0( 3»  2) =-l 

M0(4,2)=-K1 

M0(4,4)=G 

M0(5,3)=-K2 

M0(5,5)=H 

M0(6,4)=-l 

M0(6,5)=-l 

M0(6,6)=l 

M1(1,1)=B 

Ml(2,l)=-D 

M1C2,2)=F 

Ml(3,3)=l 

M1(A,A)=1 

Ml(5,5)=l 

M2(l,l)=l 

M2(2,2)=l 

BOC 1 )=A 

COMPUTE  DENOMINATOR  OF  CLOSED  LOOP 


"degree  of  matrix" 

"highest  degree  of  poly,  element1 


CALL  DTERM(1,0) 


"denom.  det.  of  closed  loop" 


c 

c 


COMPUTE  NUMERATOR  OF  CLOSED  LOOP 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 


CALL  DTERM(2,3)  "num.  det.  of  closed  loop" 

COPY  NUM.  AND  DENOM.  INTO  CLOSED  LOOP  TRANSFER  FUNCTION 
CALL  CPYPSO ,  2» 1 )  "SPTF1  is  the  closed  loop  t.f." 

CHANGE  ELEMENTS  FOR  OPEN  LOOP  CONFIGURATION 
M0(1,6)=KFLG=0  "set  KFLG  for  open  loop" 

COMPUTE  DENOMINATOR  OF  OPEN  LOOP 

CALL  DETRM(3,0)  "denom.  det.  for  open  loop" 

COMPUTE  NUMERATOR  OF  OPEN  LOOP 

CALL  DETRM(4,6)  "num.  det.  for  open  loop" 

COPY  NUM.  AND  DENOM.  INTO  OPEN  LOOP  TRANSFER  FUNCITON 
CALL  CPYPS(2,A,3)  "SPTF2  is  the  open  loop  t.f." 

DESIRED  TRANSFER  FUNCTIONS  HAVE  BEEN  COMPUTED 
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Problem:  Find  the  s-plane  transfer  function  of  the  non  rational1  function 

-ds 

A(s)  +  B(s)  e  where  A(s)  and  B(s)  are  rational  transfer  functions 
stored  in  LCAP2  as  SPTF1  and  SPTF2,  respectively. 

The  s-plane  operator  SFREQ(i)  has  a  provision  to  include  time  delay  by  cas¬ 
cading  the  response  of  SPTFi  with  time  delay  specified  by  FDLAY.  It  is  not  capa¬ 
ble  though,  of  handling  a  more  general  non  rational  function  such  as  the  one  for 
this  example. 

The  operator  FREQSC SFAUX1 ) ,  which  uses  a  user  supplied  COMPLEX  FUNCTION 
SFAUX1,  can  be  used  to  obtain  the  desired  frequency  response.  This  operator 
and  function  pair  has  been  designed  so  that  the  code  to  be  written  by  the  user 
can  be  simple.  The  usei — supplied  function  will  enable  the  user  to  (1)  compute 

-ds 

the  response  of  A(s)  and  B(s)  individually,  (2)  compute  the  response  of  e 

-ds 

and  (3)  combine  the  results  so  that  A(s)  +  B(s)  e  is  the  function  to  be 
evaluated. 

The  user  UPDATE  code  for  this  example  is: 


KIDENT  ldname 
KINSERT  START. 1 
XDECK  MAIN 
XCALL  LCAP2 

EXTERNAL  SFAUX1 

CALL  INITO 


"must  be  before  first  executable 
statement" 

"initialization  of  LCAP2  param." 


code  defining  SPTFI 
and  SPTF2 


frequency  response  parameters, 
i.e.,  FAUTO,  NOMEG,  OMEGA ( 1 ) , 
FBODE,  etc. 


CALL  FREQS(SFAUXl) 


"SFAUX1  is  the  user-supplied 
COMPLEX  FUNCTION" 


1  This  non  rational  function  could  represent  a  distributed  parameter  system. 
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CALL  LEXIT 
END 


C 

C  THE  FOLLOHING  ARE  UPDATE  CHANGES  TO  COMPLEX  FUNCTION  SFAUX1 
C 

C  COMPLEX  FUNCTION  SFAUX1  IS  IN  THE  LCAP2  SUBROUTINE  LIBRARY. 
C  USER  NEEDS  ONLY  TO  CHANGE  STATEMENT  100  TO  ADD  CODE  TO 
C  DEFINE  HIS/HER  PARTICULAR  PROBLEM. 

C 

C  USER  NEED  NOT  BE  CONCERNED  WITH  PASSING  THE  ARGUMENTS  FOR 
C  COMPUTING  THE  FREQUENCY  RESPONSE.  IT  WILL  BE  DONE  AUTO- 


C  MATICALLY. 

XDELETE  SFAUX.37  "SFAUX.37  is  the  card  ident  for 

statement  100" 

D= .  "value  of  delay  d" 

DS=-D*U  "value  of  ds" 


100  SFAUX1=SFAUX(SPTF1 )+SFAUX(SPTF2)*CEXP(CMPLX( 0 . , DS) ) 

C 

C  EXPLANATION  OF  STATEMENT  100  - 

C  SFAUXC SPTF1 )  IS  A  COMPLEX  VALUE  EQUAL  TO  THE  RESPONSE  AT 

C  COMPLEX  FREQUENCY  U.  WHERE  U  OF  COMMON/ FRQBLK/  IS  THE  S 

C  PLANE  OMEGA  FREQUENCY  BEING  VARIED  AUTOMATICALLY  BY  THE 

C  PROGRAM.  CEXP(CMPLX( 0 . >  DS) )  IS  THE  COMPLEX  VALUE  OF  THE 

C  DELAY.  SFAUX1  IS  THE  RETURNED  VALUE  OF  THE  RESPONSE  OF 
C  THE  SPECIFIED  FUNCTION. 

C 


So  that  changes  to  the  usei — supplied  COMPLEX  FUNCTION  SFAUX1  can  be  undei — 
stood,  the  code  for  this  routine  has  been  included  in  the  description  of  SFAUX1 
in  Appendix  C. 

The  code  for  COMPLEX  FUNCTION  SFAUX1  has  been  set  up  so  that  the  user  can 
access  the  first  five  s>  z  and  w  plane  transfer  functions  without  additional 
changes  to  the  code.  If  other  transfer  functions  are  to  be  used#  see  description 
of  this  routine  in  Ref.  2. 


E.  A.  Lee,  "Linear  Controls  Analysis  Program  (LCAP)  Users  Guide,"  Aerospace 
Corporation,  TOR  -  0077(2442-231-1,  5  October  1976. 

E.  A.  LEE,  "LCAP2  -  Linear  Controls  Analysis  Program,  Vol  III:  Source  Code 
Description,"  Aerospace  Corporation,  TR  -  0084(99751-1  Vol  III,  15  Novem¬ 
ber  1983. 

E.  A.  LEE,  "LCAP2  -  Linear  Controls  Analysis  Program,  Vol  II:  Interactive 
LCAP2  User's  Guide,"  Aerospace  Corporation,  TR  -  0084(99751-1  Vol  II,  15 
November  1983. 

E.  A.  Lee,  "LCAP2-Linear  Controls  Analysis  Program,"  IEEE  Control  Systems 
Magazine,  December  1982,  pp.  15-18. 


LCAP2  operators  ara  FORTRAN  subroutinea  which  ara  written  to  allow  tha  uaer 
to  aasily  porfora  control  system  analysis  operations  using  only  a  few  arguments. 
In  most  cases  tha  arguments  will  ba  indices  defining  polynomials  and/or  transfer 
functions  to  ba  operated  upon.  Descriptions  of  these  operators  are  presented  in 
alphabetical  order. 

Tha  degree  of  tha  polynomials  and  transfer  functions  referenced  and  gener¬ 
ated  by  those  LCAP2  operators  must  ba  lass  than  50. 

In  addition  to  these  operators*  tha  user  can  utilize  other  LCAP2  subrou¬ 
tines  and  functions  which  ara  described  in  Appandix  C. 


CPYPS 

Identification 

SUBROUTINE  CPYPS  -  LCAP2  Operator,  Copy  Polynomials  Into  S  Plana  Transfer 

Function 


£ 


Copy  polynomials  into  an  s  plana  transfer  functions  using  LCAP2  indices. 
For  transfer  function  evaluation  by  Cramer's  method,  this  operator  is  used  to 
define  a  transfer  function  after  two  polynomial  determinants  have  been  computed 
with  the  use  of  the  DTERM  or  DETRM  operators. 

Usaoa 

CALL  CPYPS(I, J,K) 

I  input  -  Index  of  s  plane  transfer  function  where  results  are  to 
be  stored 

J  input  -  Index  of  poly,  to  be  used  to  define  numerator  of  SPTFi,  i=I 
K  input  -  Index  of  poly,  to  ba  used  to  define  denominator  of  SPTFi,  i=I 


Method 

Copies  coefficients  of  polynomials  into  a  transfer  function.  If  the  roots 
of  tha  polynomials  ara  also  defined,  these  roots  are  also  copied  into  the  trans¬ 
fer  function. 


PREVIOUS  PAGE 

_ IS  blank 
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CPYPH 


a-. 

•VVV 


SUBROUTINE  CPYPN  -  LCAP2  Operator*  Copy  Polynomials  Into  H  Plane  Transfer 

Function 


This  operator  is  similar  to  CPYPS  except  that  it  is  for  use  with  a  w  plane 
transfer  function  instead  of  an  s  plane  transfer  function. 


CEYPZ 


SUBROUTINE  CPYPZ  -  LCAP2  Operator*  Copy  Polynomials  Into  Z  Plane  Transfer 

Function 

This  operator  is  similar  to  CPYPS  except  that  it  is  for  use  with  a  z  plan* 
transfer  function  instead  of  an  s  plane  transfer  function. 
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CPYSP 


Identification 

SUBROUTINE  CPYSP  -  LCAP2  Operator,  Copy  S  Plane  Transfer  Function  Into  Polynomials 
Purpose 

Copy  s  plane  transfer  function  into  polynomials  using  LCAP2  indices. 

Usage 

CALL  CPYSPd,  J»K) 

I  input  -  Index  of  s  plane  transfer  function  to  be  used  in  copying 

J  input  -  Index  of  polynomial  equated  with  the  numerator  of  SPTFi,  i-I 

K  input  -  Index  of  polynomial  equated  with  the  denominator  of  SPTFi,  i-I 

Method 

Copies  coefficients  of  transfer  function  into  polynomials.  If  the  roots  of 
the  transfer  function  are  available,  the  roots  are  also  stored  in  the  polynomi¬ 
als. 


cme 


Identification 

SUBROUTINE  CPYHP  -  LCAP2  Operator,  Copy  W  Plane  Transfer  Function  Into  Polynomials 

This  operator  is  similar  to  CPYSP  except  that  it  is  for  use  with  a  w  plane 
transfer  function  instead  of  an  s  plane  transfer  function. 


CPYZP 


Identification 

SUBROUTINE  CPYZP  -  LCAP2  Operator,  Copy  Z  Plane  Transfer  Function  Into  Polynomials 

This  operator  is  similar  to  CPYSP  except  that  it  is  for  use  with  a  z  plane 
transfer  function  instead  of  an  s  plane  transfer  function. 


SUBROUTINE  DTERM  -  LCAP2  Operator/  Determinant  Of  Matrix  With  Polynomial  Elements 

(New  Version) 

Purpose 

Transfer  function  evaluation  by  Cramer's  method  for  the  system  described  by 

JJts)  *ts)  =  £(s)u 

4  3  2  1 

where  ijts)  =  Ms  +  ills  +  M2  s  +  Ml  s  +  Mfi 

4  3  2  1 

Its)  =  M  s  +  Us  +  12  s  +.  II  s  +  Ifl 

MO.  Ml .  M2.  M3.  M4  are  square  matrices  of  dimension  MXM 

BO.  Bl.  B2.  B3.  B4  are  vectors  of  dimension  MXM 

£(s)  =  State  vector  of  dimension  MXM 
u  =  Scalar  input 

is  given  by 

x  (s)  det  H  (s) 

3  j 


u  det  £|(s) 


where  (s)  is  equal  to  Mts)  with  column  j  replaced  by  Its), 
j 

The  operator  DTERM  will  compute  the  determinant  of  M  ts).  Substitution  of 

3 

Its)  into  column  j  will  be  done  automatically  by  the  program. 

The  determinant  is  found  by  solving  for  its  roots  directly  and  then  com¬ 
puting  its  coefficients. 

Usage 

CALL  DTERMtI/J) 

I  input  -  Index  where  polynomial  determinant  is  to  be  stored 
J  input  -  Column  where  Its)  is  to  be  substituted  into 

tJa0  interpreted  to  mean  no  column  substitution) 


1 


••for*  DTERM  is  used*  the  Matrix  paramtars  Must  first  b*  initialized  by 
calling  H1NIT0  (only  onco) . 


2.  Matrix  paraMatars  are  in  C0MM0N/MATRIX1/.  They  era  to  ba  sat  bafora  MR00T1 
is  called.  Thasa  parasatars  ara  defined  balow> 

ParaMatars  Prasat  Description 

MXM  1  Dimansion  of  matrices  and  vectors  (1-30) 

MDEG  0  Highest  dagraa  of  polynomial  element  (0-4) 

MO  0  Matrix  for  coefficients  of  sXXO 

Ml  0  Matrix  for  coafficients  of  sxxi 

M2  0  Matrix  for  coafficiants  of  sXX2 

M3  0  Matrix  for  coafficients  of  sXX3 

M4  0  Matrix  for  coafficiants  of  sXX4 

•0  0  Vector  for  coafficiants  of  sXXO 

B1  0  Vector  for  coefficients  of  sxxi 

B2  0  Vector  for  coefficients  of  sxx2 

B3  0  Vector  for  coefficients  of  sxx3 

B4  0  Vector  for  coafficiants  of  sxx4 


Method 

If  j  is  not  zero,  £(s)  is  substituted  into  column  j  of  M(s) .  Subroutine 
MR00T1  is  than  called  to  compute  tha  determinant.  Column  j  of  H(s)  is  then 
restored  to  its  original  value. 

RtatrifitivM 

Tha  dimension  of  tha  matrix  must  not  ba  greater  than  30  x  30.  The  polynomi¬ 
al  elements  of  the  matrix  must  be  of  degree  4  or  less.  The  degree  of  the  computed 
polynomial  determinant  must  be  less  than  50. 


SUBROUTINE  FREQS  -  LCAP2  Operator*  Frequency  Response  Using  A 

User-Supplied  Function 


Purpose 

Compute  frequency  response  of  an  arbitrary  s  plane  transfer  function. 


Usage 

CALL  FREQS(FAUXl) 

FAUX1  input  -  Name  of  user  supplied  subroutine.  Must  be  declared  with  an 
EXTERNAL  statement  in  the  calling  program. 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB/ .  See  description  of 
SFREQ.  The  user  need  not  be  concerned  with  passing  the  arguments  for  com¬ 
puting  the  frequency  response.  It  is  automatically  done  by  LCAP2. 

FREQH 


Identification 

SUBROUTINE  FREQH  -  LCAP2  Operator*  H  Plane  Frequency  Response  Using  A 

User-Supplied  Function 

Purp.oat 


Compute  w  plane  frequency  r«  iponse  of  an  arbitrary  w  plane  transfer  func¬ 
tion. 

Usage 

CALL  FREQH(FAUXl) 

FAUX1  input  -  Name  of  user  supplied  subroutine.  Must  be  declared  with  an 
EXTERNAL  statement  in  the  calling  program. 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB/.  See  description  of 
SFREQ.  The  user  need  not  be  concerned  with  passing  the  arguments  for  com¬ 
puting  the  frequency  response.  It  is  automatically  done  by  LCAP2. 


FREQZ 


Identification 

SUBROUTINE  FREQZ  -  LCAP2  Operator,  Z  Plane  Frequency  Response  Using  A 

User-Supplied  Function 


EMCBffiaa 

Compute  frequency  response  of  an  arbitrary  z  plane  transfer  function. 

Usage 

CALL  FREQZ(FAUXl) 

FAUX1  input  -  Name  of  user  supplied  subroutine.  Must  be  declared  with  an 
external  statement  in  the  calling  program. 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB/ .  See  description  of 
SFREQ.  The  user  need  not  be  concerned  with  passing  the  arguments  for  com¬ 
puting  the  frequency  response.  It  is  automatically  done  by  LCAP2. 

PADD 


IdBntificfltion 

SUBROUTINE  PADD  -  LCAP2  Operator,  Polynomial  Add 
Purpose 

Add  two  polynomials  using  LCAP2  indices. 

VMflt 

CALL  PADDCI, J,K) 


I  input  -  Index  of  resultant  polynomial  sum 
J  input  -  Index  of  first  polynomial  to  be  added 
K  input  -  Index  of  second  polynomial  to  be  added 


* 


,1 


Identification 

SUBROUTINE  PEQU  -  LCAP2  Operator, 


EEfly 


Polynomial 


Equal 


■7"7  ’.y  7  r  r.  r.  r.i 


Purpose 

Equate  polynomials  using  LCAP2  indices. 


Usage 

CALL  PEQU(I, J) 

I  input  -  Index  of  resultant  polynomial 
J  input  -  Index  of  polynomial  to  be  equated  with 


PL  PC 


Identification 

SUBROUTINE  PLDC  -  LCAP2  Operator,  Polynomial  Load  In  Coefficient  Form 

Purpose 

Load  coefficients  into  polynomial  coefficient  array,  POLYi . 


Usage 


CALL  PLDC(I) 

I  input  -  Index  where  polynomial  is  to  be  stored 

1.  Polynomial  coefficients  are  entered  with  polynomial  array  POLY  (LCAP2  for¬ 
mat)  which  are  in  COMMON/HEADDB/ . 

2.  The  calling  program  must  include  COMMON/HEADDB/  and  the  appropriate  EQUIV¬ 
ALENCE  and  DIMENSION  statements  for  POLY. 

3.  The  roots  of  POLYi  will  not  be  automatically  computed.  If  this  is  desired, 
follow  this  operation  with  the  operator  PRTS(I). 
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SUBROUTINE  PL DR  -  LCAP2  Operator,  Polynomial  Load  In  Root  Fora 

Eucesst 


Load  roots  in  polynosial  root  array,  ROOT.  Aftar  tha  roots  hava  bean 
loaded*  the  coefficients  of  the  polynosial  are  computed  end  stored  in  the  poly¬ 
nosial  coefficient  array  POLYi . 


VMflt 

CALL  PLDR(I) 

I  input  -  Index  where  polynosial  is  to  be  stored 

1.  Polynosial  roots  are  entered  with  polynomial  root  array  ROOT  (LCAP2  format) 
which  is  in  COMMON/HEADDB/ . 

2.  The  calling  program  must  include  COMMON/HEADDB/  and  the  appropriate  EQUIV¬ 
ALENCE  and  DIMENSION  statements  for  ROOT. 


hse\  A  J>  I 


SUBROUTINE  PMPY  -  LCAP2  Operator*  Polynomial  Multiply 
Purpose 

Multiply  two  polynomials  using  LCAP2  indices. 

Usasut 

CALL  PMPY(I,J,K) 

I  input  -  Index  of  resultant  polynomial  product 
J  input  -  Index  of  polynomial  multiplicand 

K  input  -  Index  of  polynomial  multiplier 

Method 

If  only  the  coefficients  of  the  j-th  and  k-th  polynomials  are  available* 
the  product  is  computed  by  multiplication  of  the  coefficients.  If  the  roots  of 
the  j-th  and  k-th  polynomials  are  available*  the  product  is  computed  by  combin¬ 
ing  the  roots.  The  coefficients  of  the  product  are  then  formed  from  these  roots. 


SUBROUTINE  PPRN  -  LCAP2  Operator,  Print  Out  Polynomial 


Purpose 

Print  out  a  polynomial  using  an  LCAP2  index. 


Usage 

CALL  PPRN(I) 

I  input  -  Index  of  polynomial  to  be  printed  out 
Method 

Roots  of  the  polynomial  are  printed  out  only  if  they  are  defined  (previous¬ 
ly  computed  or  loaded  in).  Coefficients  of  the  polynomial  are  printed  out  in 
ascending  order. 


PRTS 


Identification 

SUBROUTINE  PRTS  -  LCAP2  Operator,  Find  Roots  Of  A  Polynomial 
Purpose 

Find  roots  of  a  polynomial  using  an  LCAP2  index. 


Uaeae 

CALL  PRTS(I) 

I  input  -  Index  of  polynomial 

Rertrictiena 

If  the  roots  of  POLYi  were  previously  computed  or  loaded  in,  the  program 
will  not  recompute  the  roots  from  the  coefficients.  A  message  to  this  effect 
will  be  printed. 


SUBROUTINE  PSUB  -  LCAP2  Operator,  Polynomial  Subtract 


EurpasB 

Subtract  two  polynomials  using  LCAP2  indices. 


Usaas 

CALL  PSUB(I, J,K) 

I  input  -  Index  of  resultant  polynomial  difference 

J  input  -  Index  of  first  polynomial  (minuend) 

K  input  -  Index  of  second  polynomial  (subtrahend) 


RESTORE 


Identification 

SUBROUTINE  RESTORE  -  Restore  Polynomial,  Transfer  Function  And  Matrix  Data 
E.UCP.ggfl 

Restore  polynomial,  transfer  function  and  matrix  data  from  a  previous  batch 
or  interactive  job  for  a  restart  capability  in  Batch  LCAP2. 

Usage 


CALL  RESTORE(IPRNFG) 

IPRNFG  input  -  =0  for  no  printout  of  restored  data 

File  type  for  TAPE30  must  be  declared  with  'FILE,TAPE30,BT=I . '  File  TAPE30 
must  be  attached  before  executing  LCAP2. 


SELCR 


Identification 

SUBROUTINE  SELCR  -  LCAP2  Operator,  Eliminate  Common  Roots  Of  S  Plana 

Transfer  Function 


Purpose 

Eliminate  common  roots  from  an  s  plane  transfer  function  using  an  LCAP2 
index. 

Usage 

CALL  SELCR(I) 

I  input  -  Index  of  s  plane  transfer  function 

1.  Common  root  elimination  parameters  ECRE1  (preset=2.E-4)  and  ECRE2 
(preset=l.E-8)  are  in  COMMON/HEADDB/ . 

Method 

If  a  numerator  root  nrt  and  a  denominator  root  drt  are  found  such  that 
ABS(drt/nrt  -  Cl . , 0 . )) . LT. ECRE1  for  nrt.NE.O  or  ABS(drt) . LT.ECRE2  for  nrt.EQ.O, 
roots  nrt  and  drt  are  considered  to  be  common  roots  and  will  be  eliminated  from 
the  transfer  function. 
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Identification 

SUBROUTINE  SFREQ  -  LCAP2  Operator,  S  Plana  Fraquancy  Response 

EMCBflSg 

Evaluate  s  plana  frequency  response  using  an  LCAP2  index.  Automatic  fre¬ 
quency  mode  available  to  allow  program  to  dynamically  choose  its  own  frequency 
points  to  yield  a  smooth  plot  of  the  response. 


CALL  SFREQ(I) 

I  input  -  Index  of  s  plane  transfer  function 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB .  They  are  to  be  set 
before  SFREQ  is  called.  These  parameters  are  defined  belowi 


parameter 

preset 

description 

FAUTO 

1 

.NE.O  for  automatic  frequency  mode.  Uses  NOMEG 
and  OMEGA  array. 

.EQ.O  for  user  supplied  frequency  points.  Uses 
FREQ1,  FREQ2,  . ..,  FREQS  arrays. 

NOMEG 

2 

Number  of  frequency  points  entered  in  OMEGA  array 

\ 

for  use  in  auto,  frequency  mode  (max=20) 

OMEGA 

Array  of  frequency  points  for  auto.  freq.  mode 
(Units  defined  by  RAD). 

1. 

OMEGA(l)sfirst  frequency  point  used  in  auto,  mode 

10. 

0MEGA(2)ssecond  frequency  point  used  in  auto,  mode 

0MEGA(N0MEG)=last  freq.  point  used  in  auto,  mode 

RAD 

1 

.NE.O  for  rad/sec,  .EQ.O  for  HZ 

FBODE 

1 

.NE.O  for  Bode  plot 

FNICO 

0 

•NE.O  for  Nichols  plot 

PMARG 

0 

.NE.O  for  plotting  phase  margin  instead  of  phase 
for  the  Nichols  plot 

FNYQS 

0 

.NE.O  for  Nyquist  plot 

NQDB 

0 

.NE.O  for  hardcopy  Nyquist  plot  in  db 

GRAFP 

1 

.NE.O  for  printer  (lower  resolution)  plots 

FILM 

0 

.NE.O  for  hardcopy  (higher  resolution)  plots 

FDLAY 

0 

Time  delay  (dead  time)  (s  plane  only) 

DEGMN 

-360. 

Minimum  defined  phase  in  frequency  response 
(Phase  defined  from  DEGMN  to  DEGMN+360 . ) 

CYCLE 

0 

.EQ.O  for  automatic  selection  of  2  or  3  cycle 
scale  for  Bode  plots.  (1  cycle  not  available) 
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FREQl(l)  1  Starting  fraq.  point  for  first  segment  of 

user  specified  values 

FREQK2)  10  End  freq.  point  for  second  segment  of  user 

specified  values 

FREQK3)  1  Delta  frequency  for  third  segment  of  user 

specified  values 

FREQk(l)  0  Starting  freq.  point  for  k-th  segment  X 

FREQk(2)  0  End  freq.  point  for  k-th  segment  X 

FREQkC3)  0  Delta  frequency  for  k-th  segment  x 

(k=2, 5)  (X  -  only  if  FAUTO.EQ.O) 

DBMAX  0  Maximum  db  for  plotting  frequency  response 

DBMIN  0  Minimum  db  for  plotting  frequency  response 

(Auto,  scaling  if  DBMAX=DBMIN) 

FXYDL  .5  Nyquist  plot  scale  in  units  per  inch 

(Auto,  scaling  if  FXYDL =0) 

FXYMN  -2.5  Nyquist  plot  parameter  -  minimum  real  and  imaginary 

value  plotted. 

SAMPT  1.  Sampling  period  for  use  in  w  or  z  plane  freq.  response 


Method 

If  the  automatic  frequency  mode  is  selected  (FAUT0-0),  the  program  will 
choose  frequency  values  for  evaluating  the  transfer  function  such  that  succes¬ 
sive  delta  db  and  delta  phase  values  will  be  within  specified  limits  to  yield  a 
smooth  plot.  The  program  evaluates  the  first  point  using  f  =  OMEGA(l).  Then 
choosing  deltaf  =  0MEGACD/20  initially,  the  next  frequency  to  be  used  is  com¬ 
puted  as  f  =  f  +  deltaf.  Evaluating  the  next  point  using  this  value  of  f,  the 
delta  db  and  delta  phase  is  compared  to  the  specified  limits.  If  either  is  too 
large,  deltaf  is  halved  and  the  response  is  recomputed.  If  both  are  too  small, 
deltaf  is  doubled  and  the  response  is  recomputed.  The  limits  for  the  delta  db 
response  is  EDB1/2  and  EDB1 .  The  limits  for  the  delta  phase  response  is  EDEG1/2 
and  EDEG1 .  Simultaneously  with  computing  the  next  f  to  be  used,  a  comparison  is 
made  with  the  next  value  of  OMEGA(i).  If  f  is  larger  than  OMEGA(i),  f  will  be 
replaced  with  the  value  of  OMEGA(i).  This  will  ensure  that  the  user  specified 
frequency  values  will  be  inserted  into  those  computed  by  the  program.  This  proc¬ 
ess  will  continue  until  the  last  value  OMEGAC NOMEG)  is  used. 

Since  the  plot  points  computed  to  generate  a  smooth  plot  will,  in  many  cas¬ 
es,  be  very  large,  only  a  portion  of  the  computed  response  will  be  printed  out. 
The  print  out  is  controlled  by  the  delta  db  and  delta  phase  parameters,  EDB2  and 
EDEG2,  respectively.  A  print  out  is  made  only  if  either  of  these  two  limits  are 
exceeded. 

Also  as  part  of  the  automatic  frequency  mode,  a  comparison  is  made  on  del¬ 
taf  to  keep  (deltaf/f)  within  the  limits  of  MNDH  and  MXDW.  The  lower  limit  MNDW 
is  necessary  to  prevent  an  excessive  number  of  plot  points  around  frequencies 
with  low  damping  coefficiants.  The  upper  limit  MXDW  will  ensure  enough  points  to 
yield  a  smooth  Bode  plot. 


The  abova  para  me  tars  usad  in  tha  automatic  fraquancy  mode  are  in 
COMMON/HEADDB/ .  They  can  be  changed  by  tha  user.  These  parameters  are  defind 
below t 


parameter 


preset 


description 


E0B1 

EDB2 

EDEG1 

EDEG2 

MNDW 

MXDW 

MXITF 


1.  Min.  delta  db  for  plotting 

2.  Min.  delta  db  for  printout 

4.  Max.  delta  phase  for  plotting 

10.  Max.  delta  phase  for  printout 

.0005  Min  relative  frequency  step  size 
.2  Max.  relative  frequency  step  size 

3000  Max.  no.  of  iterations  of  freq.  response 


With  either  the  automatic  or  the  non-automatic  frequency  mode,  the  program 
will  automatically  check  for  the  gain  and  phase  crossover.  When  found,  interpo¬ 
lation  is  used  to  find  the  exact  crossover  frequency  and  the  response  computed 
at  that  frequency. 

In  the  non-automatic  frequency  mode  (FAUTO=0)  the  user  can  define  up  to 
five  sets  of  frequencies  to  be  used  in  computing  the  response.  Each  of  these 
sets  are  specified  by  a  three  element  array  of  the  form  FREQk(i).  i  =  l,3 
described  above  under  Usage.  If  FREQk(l)  =  a  FREQk(2)  =  b  and  FREQk(3)  =  c.  the 
k-th  set  for  frequencies  specified  isi 


a+c,  a+2c. 


a+jc.  b 


where  j  is  the  largest  integer  such  that  (a+jc)  is  less  than  b.  Each  successive 
FREQk  array  must  define  an  increasing  set  of  frequencies  such  that  the  first 
value  of  the  segment  is  always  larger  than  the  last  value  of  the  preceding  seg¬ 
ment.  When  FREQk(3)  is  not  larger  than  FREQk(l).  as  in  tha  case  with  the  preset 
values  for  k  =  2.5.  those  segments  will  not  be  used. 
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SUBROUTINE  SLOCI  -  LCAP2  Operator,  S  Plana  Root  Locus 


P.UTRflRt 

Evaluate  s  plana  root  locus  using  an  LCAP2  index.  Automatic  gain  selection 
available  to  supplenent  user-selectod  gains. 

Uaaoa 

CALL  SLOCI(I) 

I  input  ~  Index  of  s  plana  transfer  function  to  be  evaluated 

1.  Root  locus  parameters  ara  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
SLOCI  is  called.  These  parameters  ara  defined  belowi 

parameter  preset  description 

NLOCI  2  Number  of  root  locus  gains  entered  in  array  KGAIN 

(max-25) 

KGAIN  Array  of  values  used  for  computing  root  locus  gains 

.5  KGAINC 1 )=first  user-specified  root  locus  gain 
2.  KGAIN( 2) -second  user-specified  root  locus  gain 

KGAIN( NLOCI )slast  root  locus  gain 

(Gains  computed  and  used  only  if  they  are  between 

KGAIN(l)  and  KGAIN( NLOCI)  ) 

KFLG  0  .EQ.O  to  increment  gain  by  multiplying  by  KDELT 

. NE.O  to  increment  gain  by  adding  by  KDELT 
KDELT  1.E4  Value  for  changing  gains  (preset  to  large  value  so 

that  no  additional  gains  are  computed  by  the  program) 

ITLOC  50  Max.  no.  of  different  gains  computed  for  root  locus 

GRAFP  1  .NE.O  for  printer  (low  resolution)  plot 

FILM  0  .NE.O  for  hardcopy  (high  resolution)  plot 

RLXMN  0  Minimum  x  axis  for  plot 

RLXMX  0  Maximum  x  axis  for  plot 

(Auto,  scaling  of  x  axis  if  RLXMN=RLXMX=0) 

RLYNN  0  Minimum  y  axis  for  plot 

RLYMX  0  Maximum  y  axis  for  plot 

(Auto,  scaling  of  y  axis  if  RLYMN=RlYMX=0) 

Method 

Root  locus  is  computed  by  evalueting  the  roots  of  the  polynomial  (PN  + 
GAINXPD)  where  GAIN  is  the  veried  gain  and  PN  and  PD  are  the  numerator  and  denom¬ 
inator  polynomials  of  the  transfer  function. 


SUBROUTINE  SNORM  -  LCAP2  Operator*  Normalize  S  Plana  Transfer  Function 

Puri»««t 

Normalize  s  plana  transfer  function  using  an  LCAP2  index.  Normalization  can 
be  either  with  respect  to  the  low  order  non-zero  coefficient  or  the  high  order 
coefficient  of  the  denominator. 

Osaae 


CALL  SNORN(I) 

I  input  -  Index  of  the  s  plane  transfer  function 

1.  Normalization  parametars  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
SNORM  is  called.  These  parametars  are  defined  belowi 

preset  description 

1.  Value  used  for  normalizing  tha  transfer  function 
0  If  .EQ.O,  the  low  order  non-zero  coefficient  of  the 
denominator  is  set  equal  to  the  value  of  KNORM  and 
all  other  coefficients  are  normalized  to  this  value. 

If  KN0RM=1.*  the  low  order  non-zero  coefficient  of  the 
numerator  is  the  Bode  gain. 

If  . NE.O*  the  high  order  coefficient  of  the  denominator 
is  set  equal  to  the  value  of  KNORM  and  all  other  coeffi¬ 
cients  are  normalized  to  this  value.  If  KN0RM=1.*  the 
high  order  coefficient  of  the  numerator  is  the  root 
locus  gain. 

RMtris-tinos 


parameter 

KNORM 

NRMFG 


KNORM  cannot  be  zero. 


SPADD 


Idintifisition 

SUBROUTINE  SPADD  -  LCAP2  Operator*  S  Plane  Transfer  Function  Add 
Purpose 

Add  two  s  plane  transfer  functions  usins  LCAP2  indices. 


Usage 

CALL  SPADD(I*J*K) 

I  input  -  Index  of  resultant  transfer  function  sum 
J  input  -  Index  of  first  transfer  function  to  be  added 
K  input  -  Index  of  second  transfer  function  to  be  added 


SPDIV 


Idanttfi&atian 

SUBROUTINE  SPDIV  -  LCAP2  Operator*  S  Plane  Transfer  Function  Divide 

P.urpaaa 

Divide  two  s  plane  transfer  functions  using  LCAP2  indices. 


Usage 

CALL  SPDIVCI* J*K) 


I  input  -  Index  of  resultant  transfer  function 
J  input  -  Index  of  dividend  transfer  function 
K  input  ~  Index  of  divisor  transfer  function 


SUBROUTINE  SPEQU  -  LCAP2  Operator,  S  Plane  Equal 


Eurpgss 

Equate  s  plane  transfer  functions  using  LCAP2  indices. 


Usage 

CALL  SPEQU( I *  J) 

I  input  -  Index  of  resultant  transfer  function 
J  input  -  Index  of  transfer  function  to  be  equated  with 


SPLDC 


Identification 

SUBROUTINE  SPLDC  -  LCAP2  Operator,  Load  Coefficients  Into  S  Plane  Transfer 

FUNCTION 


Purpose 

Load  coefficients  into  s  plane  transfer  function  using  an  LCAP2  index. 


Usaas 


CALL  SPLDC(I) 

I  input  -  Index  where  transfer  function  is  to  be  stored 

1.  Transfer  function  coefficients  are  entered  with  polynomial  coefficient 
arrays  POLYN  and  POLYD  (LCAP2  format)  which  are  in  CONNON/HEADDB/ .  They  are 
to  be  set  before  SPLDC  is  called. 


SUBROUTINE  SPLDR  -  LCAP2  Operator,  Load  S  Plana  Transf ar  Function  In 

Root  Form 

Purpose 

Load  roots  into  s  plana  transfar  function  using  an  LCAP2  indax. 

Usaoa 

CALL  SPLDR(I) 

I  input  -  Index  where  transfar  function  is  to  be  stored 

1.  Transfar  function  roots  are  entered  with  polynomial  coefficient  arrays 
ROOTN  and  ROOTD  (LCAP2  format)  which  are  in  COMMON/HEADDB/ .  They  are  to  ba 
sat  before  SPLDR  is  called. 


SUBROUTINE  SPMPY  -  LCAP2  Operator,  S  Plana  Transfer  Function  Multiply 
Purpose 

Multiply  two  s  plane  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  SPMPYd,  J,K) 

I  input  -  Index  of  resultant  transfer  function  product 
J  input  -  Index  of  first  transfer  function  to  be  multiplied 
K  input  -  Index  of  second  transfer  function  to  be  multiplied 

Method 

If  only  the  coefficients  of  the  j-th  and  k-th  transfer  functions  are  avail¬ 
able,  the  product  is  computed  by  multiplication  of  the  coefficients.  If  the 
roots  of  the  j-th  and  k-th  transfer  functions  are  available,  the  product  is  com¬ 
puted  by  combining  the  roots.  The  coefficients  of  the  product  are  then  formed 
from  these  roots. 


SPPRN 


Identification 

SUBROUTINE  SPPRN  -  LCAP2  Operator,  Print  Out  S  Plane  Transfer  Function 

Purpose 

Print  out  s  plane  transfer  function  using  an  LCAP2  index. 


Usage 


CALL  SPPRN(I) 

I  input  -  Index  of  transfer  function  to  be  printed  out 

tlfl.th.gd 

Roots  of  the  transfer  function  are  printed  out  only  if  they  are  defined 
(previously  computed  or  loaded  in).  Coefficients  of  the  transfer  function  are 
printed  out  in  ascending  order. 


SPRTS 


Identification 

SUBROUTINE  SPRTS  -  LCAP2  Operator,  Find  Roots  Of  S  Plane  Transfer  Function 

Purpose 

Find  roots  of  an  s  plane  transfer  function  using  an  LCAP2  index. 

Usage 

CALL  SPRTSCI) 

I  input  -  Index  of  s  plane  transfer  function 

Uilhfid 

Roots  of  the  numerator  and  denominator  are  computed  by  subroutine  PROOT . 

Riatristiwna 

If  the  roots  of  SPTFi  were  previously  computed  or  loaded  in,  the  program 
will  not  recompute  the  roots  from  the  coefficients.  A  message  to  this  effect 
will  be  printed. 


SUBROUTINE  SPSUB  -  LCAP2  Operator*  S  Plane  Transfer  Function  Subtract 
Purpose 

Subtract  two  s  plana  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  SPSUB(I* J*K) 

I  input  -  Index  of  resultant  transfer  function  difference 

J  input  -  Index  of  first  transfer  function  (minuend) 

K  input  -  Index  of  second  transfer  function  (subtrahend) 


SI  I  HE 


Idantififiitign 

SUBROUTINE  STINE  -  LCAP2  Operator/  Inverse  Laplace  Transform  And  Time  Response 

C.UM»g»t 

Compute  inverse  Laplace  transform  and  the  time  response  using  an  LCAP2 
index. 

Useoe 

CALL  STIME(I) 

I  input  -  Index  of  s  plane  transfer  function 

1.  Time  response  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
STINE  is  called.  These  parameters  are  defined  belowi 


parameter 

preset 

description 

TSTEP 

1 

.NE.O  for  step  response;  . EQ.O  for  impulse  response 

TNAGN 

1. 

Magnitude  of  input  for  time  response 

TZERO 

0 

Start  time  for  evaluating  time  response 

TEND 

1. 

End  time  for  evaluating  time  response 

TDELT 

1. 

Delta  time  for  evaluating  time  response 

TNAGN 

1. 

Magnitude  of  input  for  time  response 

GRAFP 

1 

.NE.O  for  printer  (low  resolution)  plot 

FILN 

0 

.NE.O  for  hardcopy  (high  resolution)  plot 

TXNIN 

0 

Minimum  x  axis  for  plot 

TXNAX 

0 

Maximum  x  axis  for  plot 

(Auto,  scaling  of  x  axis  if  TXMIN=TXMAX) 

TYNIN 

0 

Minimum  y  axis  for  plot 

TYNAX 

0 

Mamimum  y  axis  for  plot 

(Auto,  scaling  of  y  axis  if  TYMIN=TYMAX) 


Hathod 

The  partial  fraction  expansion  of  the  s  plane  transfer  function  times 
(1/s)/  if  the  input  is  a  step  function/  is  first  computed.  By  utilizing  the 
inverse  Laplace  transform/  the  analytical  solution  is  computed  and  printed  out. 
This  analytical  solution  is  then  evaluated  over  the  range  of  time  values  speci¬ 
fied. 

lUalcifiiiflos 

Due  to  the  algorithm  used  to  implement  the  partial  fraction  expansion/  the 
following  restrictions  on  the  form  of  the  s  plane  transfer  function  apply.  Mul¬ 
tiple  poles  are  not  allowed  except  for  those  at  the  origin.  The  poles  at  the  ori¬ 
gin  (including  the  pole  due  to  the  1/s  term  if  the  input  is  a  step)  must  be  5  or 


lass.  Also*  the  dagroa  of  the  numarator  must  not  ba  graatar  than  tha  n unbar  of 
poles  not  at  tha  origin. 


STORE 


SUBROUTINE  STORE  -  Store  Polynomial*  Transfer  Function  and  Matrix  Data 
Purpose 

Store  data  from  an  LCAPZ  batch  job  for  a  restart  capability.  This  data  can 
ba  accessed  in  a  subsequent  batch  or  interactive  job  by  using  the  RESTORE  opera¬ 
tor  . 

ilsaas 

CALL  STORE(IPRNFG) 

IPRNFG  input  -  =C  for  no  printout  of  data  to  ba  stored 

1.  To  identify  tha  data  stored*  enter  alphanumeric  information  in  HEADL64) 
through  HEADC70)  of  COMMON/HEADDB/  before  calling  STORE.  This  information 
will  be  printed  out  when  this  data  is  restored  in  a  subsequent  job. 

Method 

Data  will  be  saved  on  file  TAPE31.  The  first  record  will  be  alphanumeric 
information  copied  from  HEAD(64)  through  HEAD(70)  of  COMMON/HEADDB/.  The  other 
records  on  this  file  are  described  in  the  more  detailed  description  of  subrou¬ 
tine  STORE  in  Ref.  2. 


File  type  for  TAPE31  must  be  declared  with  'FILE,TAPE31,RT=I . *  At  tha  and 
of  tha  batch  job  the  user  must  catalog  file  TAPE31  on  tha  CDC  835  with 
* . . . *ST=PF6) *  in  tha  catalog  command. 


«N*  M  M  r 
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SUBROUTINE  SHMRX  -  LCAP2  Operator,  S  to  H  Plana  Hultirate  Transform 


Purpose 

Compute  multirate  (slow  input,  fast  output)  s  to  w  plane  transform  using 
LCAP2  indices.  The  zero  order  hold,  if  included,  is  at  the  slower  input  sampling 
rate.  The  ratio  of  the  output/input  sampling  rates  must  be  an  integer.  (note> 
the  w  is  not  the  w*  defined  by  the  Tustin's  bilinear  rule) 

Usage 

CALL  SHMRX(I, J) 

I  input  -  Index  of  w  plane  transfer  function 
J  input  -  Index  of  s  plana  transfer  function 

1 .  SHMRX  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before  SHMRX  is 
called.  These  parameters  are  defined  below: 


parameter 

DELAY 

SAMPT 

NTGER 


preset  description 

0  Time  delay,  (enter  negative  value  for  time  advance) 
1.  Sampling  period  of  the  faster  output  sampler 

1  Integer  ratio  of  output/input  sampling  rates 

(or,  input/output  sampling  periods) 

1  . NE.O  for  inclusion  of  zero  order  hold  at  the  input 


Hithsd 


Partial  fraction  expansion  of  the  s  plane  transfer  function  (including  the 
1/s  from  the  zero  order  hold  if  there  is  one)  is  computed  by  subroutine  RESDU. 
The  w  transform  (at  the  faster  output  sampling  rate)  of  each  term  of  the  expan¬ 
sion  is  then  computed.  Next,  the  terms  of  the  expansion  are  summed  and  ratio¬ 
nalized.  Roots  of  this  intermediate  transfer  function  are  then  found.  This 
result  is  multiplied  by  the  discrete  contribution  of  the  zero  order  hold  (which 
is  at  the  slower  input  sampling  rate)  to  yield  the  desired  transform. 

Bii.tristiflna 

The  algorithm  used  for  computing  the  partial  fraction  expansion  requires 
the  following  constraints  on  the  s  plane  transfer  function.  Multiple  poles  are 
not  allowed  except  for  those  at  the  origin.  The  poles  at  the  origin  (including 
the  1/s  from  the  zero  order  hold  if  there  is  one)  must  be  5  or  less.  Also,  the 
degree  of  the  numerator  must  not  be  greater  than  the  number  of  poles  not  at  the 
origin. 


*,"  v*  V 


SUBROUTINE  SHXFM  -  LCAP2  Operator,  S  to  N  Plane  Transform 


Purpose 

Compute  sampled-data  transform  from  s  to  w  plane  using  LCAP2  indices, 
(notei  this  w  is  not  the  w'  defined  by  the  Tustin's  bilinear  rule) 

Usage 


CALL  SHXFM(I*J) 

I  input  -  Index  of  w  plane  transfer  function 
J  input  -  Index  of  s  plane  transfer  function 

1.  SHXFM  parameters  are  in  COMMON/ HEADDB/ .  They  are  to  be  set  before  SHXFM  is 
called.  These  parameters  are  defined  belowi 


pa rameter 
DELAY 
SAMPT 
ZOH 


preset  description 

0  Time  delay*  (enter  negative  value  for  time  advance) 
1.  Sampling  period 

1  . NE.O  for  inclusion  of  zero  order  hold 


Method 


Partial  fraction  expansion  of  the  s  plane  transfer  function  (including  the 
1/s  from  the  zero  order  hold  if  there  is  one)  is  computed  by  subroutine  RESDU. 
The  w  transform  of  each  term  of  the  expansion  is  then  computed.  Next*  the  terms 
of  the  expansion  are  summed  and  rationalized.  Roots  of  this  intermediate  trans¬ 
fer  function  are  then  found.  This  result  is  multiplied  by  the  discrete  contrib¬ 
ution  of  the  zero  order  hold  to  yield  the  desired  transform. 


Restrictions 

The  algorithm  used  for  computing  the  partial  fraction  expansion  requires 
the  following  constraints  on  the  s  plane  transfer  function.  Multiple  poles  are 
not  allowed  except  for  those  at  the  origin.  The  poles  at  the  origin  (including 
the  1/s  from  the  zero  order  hold  if  there  is  one)  must  be  5  or  less.  Also*  the 
degree  of  the  numerator  must  not  be  greater  than  the  number  of  poles  not  at  the 
origin. 


sum 


Identification 

SUBROUTINE  SZMRX  -  l  CAP  2  Operator*  S  to  Z  Plano  Multi rata  Tranafom 


Purgflit 

Compute  aultirata  (alow  input*  fast  output)  a  to  z  plana  transfora  using 
LCAP2  indices.  Tha  zero  order  hold*  if  included*  is  at  the  slower  input  sampling 
rata.  Tha  ratio  of  tha  output/input  saapling  rates  aust  be  an  integer. 


yauut 

CALL  SZMRX(I* J) 


I  input  -  Indax  of  z  piano  transfar  function 
J  input  -  Indax  of  s  plana  transfar  function 

1 .  SZMRX  paraaotars  are  in  COMMON/HEADDB/ .  They  are  to  ba  set  before  SZMRX  is 
called.  Those  paraaeters  are  defined  belowi 


paraaatar 

DELAY 

SAMPT 

NTOER 

ZOH 


preset  description 

0  Tiae  delay*  (enter  negative  value  for  time  advance) 
1.  Sampling  period  of  the  faster  output  sampler 
1  Integer  ratio  of  output/input  sampling  rates 
(or,  input/output  sampling  periods) 

1  .NE.O  for  inclusion  of  zoro  ordor  hold  at  the  input 


Partial  fraction  expansion  of  the  s  plane  transfer  function  (including  the 
1/s  from  tha  zoro  order  hold  if  there  is  one)  is  computed  by  subroutine  RESDU. 
Tha  w  transfora  (at  tho  faster  output  saapling  rate)  of  each  term  of  the  expan¬ 
sion  is  then  computed.  Next*  the  terms  of  the  expansion  are  summed  and  ratio¬ 
nalized.  Roots  of  this  intermediate  transfer  function  are  then  found.  This 
result  is  multiplied  by  the  discrete  contribution  of  the  zero  order  hold  (which 
is  at  tha  slower  input  sampling  rata)  to  yield  the  w  plane  form  of  the  desired 
transform.  Next*  subroutine  HZTRANS  is  called  to  transform  the  w  plane  roots  to 
tha  z  plana.  Tha  coefficients  of  the  z  plane  transfer  function  are  then 
computed . 


RtBtricUflrm 

The  algorithm  use  for  computing  tha  partial  fraction  expansion  requires 
that  the  following  constraints  on  the  s  plane  transfer  function.  Multiple  poles 
are  not  allowed  except  for  those  at  the  origin.  The  poles  at  the  origin  (includ¬ 
ing  tha  1/s  from  tha  zero  order  hold  if  there  is  one)  must  be  5  or  less.  Also* 
tho  dagrao  of  tho  numerator  must  not  be  greater  than  the  number  of  poles  not  at 
tha  origin. 


.■wr;  v 
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SZXFM 

Identification 

SUBROUTINE  SZXFM  -  LCAP2  Operator,  S  To  Z  Plane  Transform 
Purpose 

Coepute  sampl ad-data  transfora  froa  a  to  z  plane  using  LCAP2  indices. 

Usage 

CALL  SZXFM( I ,  J  ) 

I  input  -  Index  of  z  plana  transfer  function 
J  input  -  Index  of  s  plana  transfer  function 

1.  SZXFN  parameters  are  in  COMMON/ HEADDB/.  Sea  description  for  SMXFM. 

Method 

Partial  fraction  expansion  of  s  plana  transfer  function  (including  tha  1/s 
froa  tha  zero  order  hold  if  there  is  one)  is  computed  by  subroutine  RESDU.  Since 
calculations  performed  in  the  w  plane  are  more  accurate  than  in  the  z  plane,  the 
w  transform  of  each  term  of  the  partial  fraction  expansion  is  computed.  Roots  of 
this  intermediate  transfer  function  are  then  found.  This  result  is  multiplied  by 
the  discrete  contribution  of  the  zero  order  hold  to  yield  the  w  plane  form  of  the 
desired  transform.  Next,  subroutine  MZTRANS  is  called  to  transform  the  w  plane 
roots  to  the  z  plane.  The  coefficients  of  the  z  plane  transfer  function  are  then 
computed. 

Rtatrictiana 

Same  restrictions  on  s  plane  transfer  functions  that  apply  for  SMXFM.  For 
higher  order  s  plane  transfer  functions,  the  w  plane  transform  will  be  more 
accurate  than  the  z  plane  transform.  To  determine  differences  in  numerical  accu¬ 
racies,  compute  both  w  and  z  plane  transforms  and  compare  frequency  responses. 


HELCR 


Idantifi cation 

SUBROUTINE  HELCR  -  LCAP2  Operator,  Eliminate  Common  Roots  Of  N  Plane 

Transfer  Function 

Purpose 

Eliminate  common  roots  from  a  w  plane  transfer  function  using  an  LCAP2 
index. 

Usage 

CALL  HELCR(I) 

I  input  -  Index  of  w  plane  transfer  function 

1.  Common  root  elimination  parameters  ECRE1  (preset=2 . E-4)  and  ECRE2 
(preset-l.E-8)  are  in  COMMON/HEADDB/ . 

Method 

If  a  numerator  root  nrt  and  a  denominator  root  drt  are  found  such  that 
ABSCdrt/nrt  -  Cl . ,0. )) . LT.ECRE1  for  nrt.NE.O  or  ABSCdrt) .LT.ECRE2  for  nrt.EQ.O, 
roots  nrt  and  drt  are  considered  to  be  common  roots  and  will  be  eliminated  from 
the  transfer  function. 


SUBROUTINE  WFREQ  -  LCAP2  Operator,  W  Plan*  Frequency  Response 


Purpose 

Evaluate  w  plane  frequency  response  using  an  LCAP2  index.  Automatic  fre¬ 
quency  mode  available  to  allow  program  to  dynamically  choose  its  own  frequency 
points  to  yield  a  smooth  plot  of  the  response. 

Usaae 

CALL  HFREQ(I) 

I  input  -  Index  of  w  plane  transfer  function 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB .  They  are  to  be  set 
before  HFREQ  is  called.  See  description  of  SFREQ  for  the  complete  list  of 
definitions  of  these  parameters.  The  parameter  SAMPT  is  described  belowi 

parameter  preset  description 

SAMPT  1  Sampling  period 


Method 

Same  as  that  described  in  detail  in  description  of  SFREQ. 

In  the  automatic  frequency  mode  (FAUTO=0)  the  program  will  avoid  using  fre¬ 
quency  values  at  or  near  the  half  sampling  frequency.  H  plane  frequencies  (im¬ 
aginary  part)  greater  than  1000.  will  not  be  used. 


’  >  T, 


HI  PCI 

Identification 

SUBROUTINE  HLOCI  -  LCAP2  Operator,  H  Plane  Root  Locus 
Purpose 

Evaluate  w  plana  root  locus  using  an  LCAP2  index.  Automatic  gain  selection 
available  to  supplement  user-selected  gains. 

Usaoa 

CALL  HLOCI(I) 

I  input  -  Index  of  w  plane  transfer  function  to  be  evaluated 

1.  Root  locus  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before  HLO¬ 
CI  is  called.  See  description  of  SLOCI  for  a  complete  list  of  definition  of 
these  parameters. 

Method 

Same  as  that  described  in  detail  in  description  of  SLOCI. 


SUBROUTINE  WMRFQ  -  LCAP2  Operator*  H  Plane  Multirate  Frequency  Response 
Purpose 

Evaluate  multirate  (fast  input*  slow  output)  frequency  response  of  a  w 
plane  transfer  function  using  an  LCAP2  index. 

Usaag 

CALL  HMRFQ(I,M) 

I  input  -  Index  of  w  plane  transfer  function 
M  input  -  Integer  ratio  of  output/input  sampling  periods 

1.  The  input  w  plane  transfer  function  is  at  tha  faster  sampling  rate. 

2.  Frequency  response  parameters  are  in  COMMON/HEADDB/ .  See  description  of 
SFREQ. 

3.  The  sampling  period*  SAMPT,  is  for  the  slower  output  sampler. 

ftettwd 

The  frequency  response  is  evaluated  by  direct  application  of  Sklansky's 
frequency  decomposition.  No  explicit  rational  representation  of  the  slower  out¬ 
put  transform  is  computed.  If  an  explicit  representation  of  the  slower  output 
transfer  function  is  desired*  see  LCAP2  operator  HMRXFM. 


HMRXFM 


IdinUficatien 

SUBROUTINE  HMRXFM  -  LCAP2  Operator,  Multirate  (fast  input,  slow  output 

sampler)  H  Transform 


Purpasa 

Compute  the  output  w  transform  of  a  fast  to  slow  sampler  using  LCAP2 
indices.  This  operation  will  yield  a  rational  transfer  function  at  the  slower 
sampling  rate. 

Usage 

CALL  WMRXFM( I , J ) 

I  input  -  Index  of  resultant  slower  output  w  plane  transfer  function 
J  input  -  Index  of  faster  input  w  plane  transfer  function 

1.  The  integer  ratio,  output/input  sampling  periods,  NTGER,  of  COMMON/HEADDB/ 
must  be  set  before  this  subroutine  is  called. 

The  faster  input  transfer  function  is  first  transformed  to  the  z  plane. 
The  output  transfer  function  of  a  fast  to  slow  rate  sampler  is  then  given  in  the 
z  plane  by  Sklansky's  frequency  decomposition  method  as, 

ail 

1  \ 

\  T/n  j2*pi*k/n 
n  /  G  (z  e  ) 

/ _  n 

k=l 

T/n 

where  G  (z  )  is  the  z  plane  transfer  function  at  the  fester  sampling  rate 
n 

and  n  is  the  integer  ratio  of  the  output/input  sampling  periods.  Using  the 
root  form  representation  of  the  input  transform,  a  rational  representation 
of  the  slower  rate  output  trensfer  function  is  computed.  This  transfer  function 
is  then  transformed  to  the  w  plane. 


SUBROUTINE  NNORM  -  LCAP2  Operator,  Normalize  W  Plane  Transfer  Function 


P»r.P<?SB 

Normalize  w  plane  transfer  function  using  an  LCAP2  index.  Normalization  can 
be  either  with  respect  to  the  low  order  non-zero  coefficient  or  the  high  order 
coefficient  of  the  denominator. 

Usage 


CALL  HNORM(I) 


I  input  -  Index  of  the  w  plane  transfer  function 

1.  Normalization  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
HNORM  is  called.  These  parameters  are  defined  below t 


parameter  preset  description 

KNORM  1.  Value  used  for  normalizing  the  transfer  function 

NRMFG  0  If  .EQ.O,  the  low  order  non-zero  coefficient  of  the 

denominator  is  set  equal  to  the  value  of  KNORM  and 
all  other  coefficients  are  normalized  to  this  value. 

If  KN0RM=1.,  the  low  order  non-zero  coefficient  of  the 
numerator  is  the  Bode  gain. 

If  .NE.O,  the  high  order  coefficient  of  the  denominator 
is  set  equal  to  the  value  of  KNORM  and  all  other  coeffi¬ 
cients  are  normalized  to  this  value.  If  KN0RM=1.,  the 
high  order  coefficient  of  the  numerator  is  the  root 
locus  gain. 

R«rtri<rti9P9 


KNORM  cannot  be  zero . 


SUBROUTINE  NPADD  -  LCAP2  Operator,  H  Plana  Transfer  Function  Add 


PUTPgBt 

Add  two  w  plana  transfer  functions  using  LCAP2  indices. 


CALL  HPADDCI, J,K) 

I  input  -  Index  of  resultant  transfer  function  sum 
J  input  -  Index  of  first  transfer  function  to  be  added 
K  input  -  Index  of  second  transfer  function  to  be  added 


NPDIV 


IdtnttficitiQn 

SUBROUTINE  NPDIV  -  LCAP2  Operator,  N  Plane  Transfer  Function  Divide 

P-tlTEMt 

Divide  two  w  plane  transfer  functions  using  LCAP2  indices. 


Usage 

CALL  HPDIVd, J,K) 

I  input  -  Index  of  resultant  transfer  function 
J  input  -  Index  of  dividend  transfer  function 
K  input  -  Index  of  divisor  transfer  function 


SUBROUTINE  HPEQU  -  LCAP2  Operator,  H  Plana  Equal 
Purpose 

Equata  w  plana  transfar  functions  using  LCAP2  indicas. 


Usaoa 

CALL  HPEQU(I,J) 

I  input  -  Index  of  resultant  transfer  function 
J  input  -  Index  of  transfar  function  to  ba  equated  with 


WPLDC 

Identification 

SUBROUTINE  HPLDC  -  LCAP2  Operator,  Load  Coefficients  Into  H  Plana  Transfar 

Function 

Purpose 

Load  coefficients  into  w  plana  transfar  function  using  an  LCAP2  index. 

Usage 

CALL  HPLDC(I) 

I  input  -  Index  where  transfar  function  is  to  ba  stored 


1. 


Transfar  function  coefficients  are  entered  with  polynomial  coefficient 
arrays  POLYN  and  POLYO  ( LCAP2  format)  which  are  in  COMMON/HEADDB/ .  They  ara 
to  ba  sat  before  NPLDC  is  called. 


SUBROUTINE  HPLDR  -  LCAP2  Operator,  Load  H  Plana  Transfer  Function  In 

Root  Form 

Purpose 

Load  roots  into  w  plane  transfer  function  using  an  LCAP2  index. 

Usage 

CALL  HPLDR(I) 

I  input  -  Index  where  transfer  function  is  to  be  stored 

1.  Transfer  function  roots  are  entered  with  polynomial  coefficient  arrays 
ROOTN  and  ROOTD  (LCAP2  format)  which  are  in  COMMON/HEADDB/ .  They  are  to  be 
set  before  NPLDR  is  called. 

HPMPY 


SUBROUTINE  HPMPY  -  LCAP2  Operator,  M  Plane  Transfer  Function  Multiply 

Purpgs* 

Multiply  two  w  plane  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  HPMPYd, J,K) 

I  input  -  Index  of  resultant  transfer  function  product 

J  input  -  Index  of  first  transfer  function  to  be  multiplied 

K  input  -  Index  of  second  transfer  function  to  be  multiplied 

Hiihsd 

If  only  the  coefficients  of  the  j-th  end  k-th  transfer  functions  are  avail¬ 
able,  the  product  is  computed  by  multiplication  of  the  coefficients.  If  the 
roots  of  the  j-th  and  k-th  transfer  functions  are  available,  the  product  is  com¬ 
puted  by  combining  the  roots.  The  coefficients  of  the  product  are  then  formed 
from  these  roots. 


i 


HPPRN 


Identification 

SUBROUTINE  WPPRN  -  LCAP2  Operator,  Print  Out  H  Plane  Transfer  Function 
Purpose 

Print  out  w  plane  transfer  function  using  an  LCAP2  index. 


CALL  HPPRN(I) 

I  input  -  Index  of  transfer  function  to  be  printed  out 
Method 

Roots  of  the  transfer  function  are  printed  out  only  if  they  are  defined 
(previously  computed  or  loaded  in).  Coefficients  of  the  transfer  function  are 
printed  out  in  ascending  order. 


HPRTS 


Idaniifigalian 

SUBROUTINE  HPRTS  -  LCAP2  Operator,  Find  Roots  Of  H  Plane  Transfer  Function 
PurP9SB 


SUBROUTINE  NPSUB  -  LCAP2  Operator*  H  Plana  Transfer  Function  Subtract 
Purpose 

Subtract  two  w  plane  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  HPSUB(I,J,K) 

I  input  -  Index  of  resultant  transfer  function  difference 

J  input  -  Index  of  first  transfer  function  (minuend) 

K  input  -  Index  of  second  transfer  function  (subtrahend) 


SUBROUTINE  HSXFM  -  LCAP2  Operator,  Transform  H  Plana  Roots  Into  S  Plana 


purposa 

Transform  w  plana  roots  into  "equivalent”  s  plana  roots  using  an  LCAP2 
index.  The  transformation  of  the  w  plana  roots  to  the  s  plane  is  not  unique.  The 
"equivalent"  s  plana  roots  are  provided  solely  to  aid  the  analyst  in  identifying 
and  correlating  w  plane  roots.  The  computed  s  plane  roots  are  not  saved. 

Usage 


CALL  HSXFM(I) 

I  input  -  Index  of  w  plana  transfer  function 

1.  Sampling  period,  SAMPT,  of  COMMON/HEADDB/  must  be  set  before  calling  this 
subroutine. 

Method 

Transformation  of  tha  roots  from  w  to  the  s  plane  is  defined  by 
s  =  ln(  (l+w)/(l-w)  )  /  SAMPT 

Hhen  w  =  -1.0  or  +1.0  the  "equivalent"  s  plana  root  is  undefined.  If  ABS(l.-w) 
is  lass  than  l.R-5,  the  equivalent  root  is  printed  out  as  999999.99.  If 
ABS(w+l.)  is  less  than  l.E-S,  the  equivalent  root  is  printed  out  as  -999999.99. 


b&XFM 


SUBROUTINE  WZXFM  -  LCAP2  Operator,  N  to  Z  Plane  Transformation 
Purpose 

Compute  w  to  z  plane  bilinear  transformation  using  an  LCAP2  index,  (notei 
the  w  is  not  the  w'  defined  by  the  Tustin's  bilinear  rule) 

Usage 

CALL  HZXFMCI,  J) 

I  input  -  Index  of  computed  z  plane  transfer  function 
J  input  -  Index  of  w  plane  transfer  function  to  be  transformed. 

1.  The  sampling  period,  SANPT,  of  CONMON/HEADDB/  must  be  set  before  calling 
this  subroutine. 

Method 

Bilinear  transformation  is  implemented  by  transformation  of  the  w  plane 
roots. 
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ZELCR 


Identification 

SUBROUTINE  ZELCR  -  LCAP2  Operator*  Eliminate  Common  Roots  Of  Z  Plane 

Transfer  Function 


Purpose 

Eliminate  common  roots  from  a  z  plane  transfer  function  using  an  LCAP2 
index. 

Usage 


CALL  ZELCR(I) 

I  input  -  Index  of  z  plane  transfer  function 

1.  Common  root  elimination  parameters  ECRE1  (preset=2 . E-4)  and  ECRE2 
( preset-1. E-S)  are  in  COMMON/HEADDB/ . 

Methgd 

If  a  numerator  root  nrt  and  a  denominator  root  drt  are  found  such  that 
ABS(drt/nrt  -  ( 1 . , 0 . ) ) . LT . ECRE1  for  nrt.NE.O  or  ABS(drt) . LT . ECRE2  for  nrt.EQ.O, 
roots  nrt  and  drt  are  considered  to  be  common  roots  and  will  be  eliminated  from 
the  transfer  function. 
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SUIROUTINE  ZFREQ  -  LCAP2  Operator*  Z  Plano  Proquancy  Rasponsa 


Purpaaa 

Evaluate  z  plana  frequency  rasponsa  using  an  LCAP2  index.  Automatic  fre¬ 
quency  mode  available  to  allow  program  to  dynamically  choose  its  own  frequency 
points  to  yield  a  smooth  plot  of  the  rasponsa. 

Usama 

CALL  ZPREQ(I) 

I  input  -  Index  of  z  plana  transfer  function 

1.  Frequency  rasponsa  parameters  are  in  COMHON/HEADDB .  They  ara  to  be  sat 
before  ZFREQ  is  called.  See  description  of  SFREQ  for  the  complete  list  of 
definitions  of  these  parameters.  The  parameter  SAMPT  is  described  below t 

parameter  preset  description 

SAMPT  1  Sampling  period 


Httlffld 


Same  as  that  described  in  datail  in  description  of  SFREQ 


SUBROUTINE  ZLOCI  -  LCAP2  Operator,  Z  Plana  Root  Locus 


Puroosa 

Evaluata  z  plana  root  locus  using  an  LCAP2  indax.  Automatic  gain  salaction 
avai labia  to  supplamant  usar-salactod  gains. 

Usaga 


CALL  ZLOCI(I) 

I  input  -  Indax  of  z  plana  transfar  function  to  ba  avaluatad 

1.  Root  locus  parameters  are  in  COMMON/HEADDB/.  They  are  to  ba  sat  before  ZLO¬ 
CI  is  called.  Saa  description  of  SLOCI  for  a  complete  list  of  definitions 
of  these  parameters. 

Hathad 

Same  as  that  described  in  detail  in  description  of  SLOCI. 


SUBROUTINE  ZMRFQ  -  LCAP2  Operator,  Z  Plane  Multirate  Frequency  Response 
Purpose 

Evaluate  multirate  (fast  input,  slow  output)  frequency  response  of  a  z 
plane  transfer  function  using  an  LCAP2  index. 

Usage 

CALL  ZMRFQ(I,M) 

I  input  -  Index  of  z  plane  transfer  function 
M  input  -  Integer  ratio  of  output/input  sampling  periods 

1.  The  input  z  plane  transfer  function  is  at  the  faster  sampling  rate. 

2.  Frequency  response  parameters  are  in  COMMON/HEADDB/ .  See  description  of 
SFREQ . 

3.  The  sampling  period,  SAMPT,  is  for  the  slower  output  sampler. 

Hgthftd 

The  frequency  response  is  evaluated  by  direct  application  of  Sklansky's 
frequency  decomposition.  No  explicit  rational  representation  of  the  slower  out¬ 
put  transform  is  computed.  If  an  explicit  representation  of  the  slower  output 
transfer  function  is  desired,  see  LCAP2  operator  ZMRXFM. 


ZMRXFM 


Identification 

SUBROUTINE  ZMRXFM  -  LCAP2  Operator,  Multirate  (fast  input,  slow  output 

sampler)  Z  Transform 


Purpose 

Compute  the  output  z  transform  of  a  fast  to  slow  sampler  using  LCAP2 
indices.  This  operation  will  yield  a  rational  transfer  function  at  the  slower 
sampling  rate. 

Usage 


CALL  ZMRXFMd, J) 

I  input  -  Index  of  resultant  slower  output  z  transfer  function 
J  input  -  Index  of  faster  input  z  transfer  function 

1.  The  integer  ratio,  output/input  sampling  periods,  NTGER,  of  COMMON/HEADDB/ 
must  be  set  before  this  subroutine  is  called. 

Method 

The  output  transform  of  a  fast  to  slow  rate  sampler  is  given  by  Sklansky's 
frequency  decomposition  method  as 

nzl 

1  \ 

\  T/n  j2XpiXk/n 
n  /  G  (z  e  ) 

/ _  n 

k=l 

T/n 

where  G  (z  )  is  the  z  plane  transfer  function  at  the  faster  sampling  rate 
n 

and  n  is  the  integer  ratio  of  the  output/input  sampling  periods.  Using  the 
root  form  representation  of  the  input  transform,  a  rational  representation 
of  the  slower  rate  output  transform  is  computed. 


SUBROUTINE  ZNORM  -  LCAP2  Operator,  Normalize  Z  Plana  Transfer  Function 


Purpose 

Normalize  z  plana  transfer  function  using  an  LCAP2  index.  Normalization  can 
ba  either  with  respect  to  the  low  order  non-zaro  coefficient  or  the  high  order 
coefficient  of  the  denominator. 

Usage 

CALL  ZNORH(I) 

I  input  -  Index  of  the  z  plane  transfer  function 

1.  Normalization  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
ZNORM  is  called.  These  parameters  are  defined  belowi 

parameter  preset  description 

KNORN  1.  Value  used  for  normalizing  the  transfer  function 

NRMFO  0  If  .EQ.O,  the  low  order  non-zero  coefficient  of  the 

denominator  is  set  equal  to  the  value  of  KNORM  and 
all  other  coefficients  are  normalized  to  this  value. 

If  .NE.O,  the  high  order  coefficient  of  the  denominator 
is  set  equal  to  the  value  of  KNORM  and  all  other  coeffi¬ 
cients  are  normalized  to  this  value. 

Rwtrifctivna 


KNORM  cannot  be  zero. 


SUBROUTINE  ZPADO  -  LCAP2  Operator,  Z  Plana  Transfar  Function  Add 
Puraoaa 

Add  two  z  plana  transfar  functions  using  LCAP2  indices. 

Usaoa 

CALL  ZPADD(I,J,K) 

I  input  -  Index  of  resultant  transfar  function  sun 
J  input  -  Index  of  first  transfar  function  to  bo  added 
K  input  -  Index  of  second  transfer  function  to  be  added 


ZPDIV 

Idtnttficitign 

SUBROUTINE  ZPDIV  -  LCAP2  Operator,  Z  Plana  Transfar  Function  Divide 
PurPORt 

Divide  two  z  plana  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  ZPDIVCI, J#K) 

I  input  -  Index  of  resultant  transfar  function 
J  input  -  Index  of  dividend  transfar  function 
K  input  -  Index  of  divisor  transfer  function 


SUBROUTINE  ZPEQU  -  LCAP2  Operator,  Z  Plana  Equal 


Purpoaa 

Equata  z  plana  transfar  functions  using  LCAP2  indicas. 


Ussaa 

CALL  ZPEQUCI,J) 

I  input  -  Indax  of  rasultant  transfar  function 
J  input  -  Indax  of  transfar  function  to  ba  aquatad  with 


ZPLDC 


IdtntifistUgn 

SUBROUTINE  ZPLDC  -  LCAP2  Operator,  Load  Coafficiants  Into  Z  Plano  Transfar 

Function 


PilCPORt 

Load  coafficiants  into  z  plana  transfar  function  using  an  LCAP2  indax. 


Ussaa 

CALL  ZPLDC(I) 

1  input  -  Indax  whara  transfar  function  is  to  ba  storad 

1.  Transfar  function  coafficiants  ara  ontsrod  with  polynomial  coafficiant 
arrays  POLYN  and  POLYD  (LCAP2  format)  which  ara  in  COMMON/HEADDB/ .  They  ara 
to  bo  sat  bafora  ZPLDC  is  callad. 


SUBROUTINE  ZPLDR  -  LCAP2  Operator*  Load  Z  Plane  Transfer  Function  In 

Root  Form 

Purpose 

Load  roots  into  z  plane  transfer  function  using  an  LCAP2  index. 

Usage 

CALL  ZPLDR(I) 

I  input  -  Index  where  transfer  function  is  to  be  stored 

1.  Transfer  function  roots  are  entered  with  polynomial  coefficient  arrays 
ROOTN  and  ROOTD  (LCAP2  format)  which  are  in  COMMON/HEADDB/ .  They  are  to  be 
set  before  ZPLDR  is  called. 


SUBROUTINE  ZPMPY  -  LCAP2  Operator,  Z  Plan*  Transfer  Function  Multiply 
Purpose 

Multiply  two  x  plan*  tranafar  function*  uains  LCAP2  indicaa. 

CALL  ZPMPYCI* J»K> 

X  input  -  Index  of  resultant  tranafar  function  product 
J  input  -  Xndax  of  first  transfer  function  to  ba  Multiplied 
K.  input  -  Indax  of  second  transfer  function  to  ba  Multiplied 

Method 

If  only  tha  coofficiants  of  th*  j-th  and  k-th  tranafar  functions  are  avail¬ 
able*  the  product  is  conputod  by  Multiplication  of  the  coefficients.  If  th* 
roots  of  tha  j-th  and  k-th  tranafar  functions  are  available*  th*  product  is  coa- 
puted  by  coabining  the  roots.  Tha  coefficients  of  tha  product  are  than  fornad 
fro*  thasa  roots. 


SUBROUTINE  ZPPRN  -  LCAP2  Operator,  Print  Out  Z  Plana  Transfer  Function 


Purpose 

Print  out  z  plana  transfer  function  using  an  LCAPZ  index. 


Usage 

CALL  ZPPRN(I) 

I  input  -  Index  of  transfer  function  to  be  printed  out 

Hftlhfid 

Roots  of  the  transfer  function  aro  printed  out  only  if  they  are  defined 
(previously  computed  or  loaded  in).  Coefficients  of  the  transfer  function  are 
printed  out  in  ascending  order. 


ZPRTS 

Identification 

SUBROUTINE  ZPRTS  -  LCAP2  Operator,  Find  Roots  Of  Z  Plane  Transfer  Function 
PUCRgRl 

Find  roots  of  a  z  plana  transfer  function  using  an  LCAP2  index. 

Useoe 

CALL  ZPRTS(I) 

I  input  -  Index  of  z  plane  transfer  function 

tkibad 

Roots  of  the  numerator  and  denominator  are  computed  by  subroutine  PROOT. 

The  code  for  this  routine  is  in  subroutine  SPRTS. 

Raaicisliaim 

If  the  roots  of  ZPTFi  were  previously  computed  or  loaded  in,  the  program 
will  not  recompute  the  roots  from  the  coefficients.  A  message  to  this  effect 
will  be  printed. 


ZPSUB 


IdiDtifiotign 

SUBROUTINE  ZPSUB  -  LCAP2  Operator,  Z  Plana  Transfer  Function  Subtract 
Purpose 

Subtract  two  z  plane  transfer  functions  using  LCAP2  indices. 

Usage 

CALL  ZPSUB(I, J,K) 

1  input  -  Index  of  resultant  transfer  function  difference 

J  input  -  Index  of  first  transfer  function  (minuend) 

K  input  ~  Index  of  second  transfer  function  (subtrahend) 


ZSXFM 


Identification 

SUBROUTINE  ZSXFM  >  LCAP2  Operator,  Transfors  Z  Plane  Roots  Into  S  Plane 

Eurjp«« 

Transfors  z  plane  roots  into  "equivalent"  s  plane  roots  using  an  LCAP2 
index.  The  transforsation  of  the  z  plane  roots  to  the  s  plane  is  not  unique.  The 
"equivalent"  s  plane  roots  are  provided  solely  to  aid  the  analyst  in  identifying 
and  correlating  z  plane  roots.  The  computed  s  plane  roots  are  not  saved. 


CALL  ZSXFM(I) 

I  input  -  Index  of  z  plane  transfer  function 

1.  Sampling  period,  SAMPT,  of  COMMON/HEADDB/  must  be  set  before  calling  this 
subroutine . 

Method 

Transformation  of  the  roots  from  z  to  the  s  plane  is  defined  by 


a  ■  InC  z  )  /  SAMPT 


SUBROUTINE  ZTIME  -  LCAP2  Operator,  Inverse  Z  Transform  And  Tima  Response 


Purpqsa 

Compute  inverse  z  transform  and  the  time  response  using  an  LCAP2  index. 


Usaae 

CALL  ZTIME(I) 

I  input  -  Index  of  z  plane  transfer  function 

1.  Time  response  parameters  are  in  COMMON/HEADDB/ .  They  are  to  be  set  before 
ZTIME  is  called.  These  parameters  are  defined  belowi 


parameter 

preset 

description 

TSTEP 

1 

.NE.O  for  step  response;  .EQ.O  for  impulse  response 

TMA6N 

1. 

Magnitude  of  input  for  time  response 

TEND 

1. 

End  time  for  evaluating  time  response 

TMAGN 

1. 

Magnitude  of  input  for  time  response 

SAMPT 

1. 

Sampling  period 

GRAFP 

1 

.NE.O  for  printer  (low  resolution)  plot 

FILM 

0 

.NE.O  for  hardcopy  (high  resolution)  plot 

TXMIN 

0 

Minimum  x  axis  for  plot 

TXMAX 

0 

Maximum  x  axis  for  plot 

(Auto,  scaling  of  x  axis  if  TXMIN=TXMAX) 

TYMIN 

0 

Minimum  y  axis  for  plot 

TYMAX 

0 

Mamimum  y  axis  for  plot 

(Auto,  scaling  of  y  axis  if  TYMIN=TYMAX) 

Method 

The  inverse  z  transform  is  computed  by  the  power  series  (long  division) 
method.  Hhile  this  method  of  computing  the  time  response  is  inherently  less 
sccurete  then  the  pertiel  frection  method,  results  for  typical  transfer  func¬ 
tions  are  excellent.  To  provide  a  measure  of  the  accuracy  of  the  response,  the 
results  are  computed  in  double  precision  and  compared  with  single  precision 
results. 


SUBROUTINE  ZVCNO  -  LCAPZ  Operator,  Z  to  <ZXXn>  Transformation 
Purpose 

Computa  transformation  of  z  plana  transfar  function  to  a  faster  z  variable 
by  replacement  of  variables. 

Usage 

CALL  ZVCNG( I  *  J ) 

I  input  -  Index  of  resultant  z  plana  transfer  function  expressed 
in  terms  of  the  faster  z  variable 
J  input  -  Index  of  z  plana  transfar  function  to  be  operated  upon 

1.  The  integer  ratio  of  faster/slowar  sampling  rata,  NTGER,  of  COMMON/HEADDB/ 
must  ba  sat  before  ZVCNG  is  called. 

Utihfid 

The  z  variable  of  tha  slower  sampled  transfer  function  is  replaced  by  zXXn 
and  stored  into  tha  faster  sampled  transfer  function.  This  subroutine  calls 
ZVCHNG1 . 

Restrictions 

Since  tha  format  used  to  represent  transfar  function  arrays  in  LCAPZ  limits 
tha  degree  of  the  polynomials  to  lass  that  50,  the  degree  of  the  j-th  z  plane 
transfer  function  times  NTGER  must  be  less  than  50. 


SUBROUTINE  ZHXFM  -  LCAP2  Operator,  Z  to  H  Plano  Transformation 
Purooaa 

Compute  z  to  w  plana  bilinear  transformation  using  an  LCAP2  index,  (notei 
the  w  is  not  the  w'  defined  by  the  Tastin' s  bilinear  rule) 

Usaoa 

CALL  ZMXFM( I ,  J  ) 

I  input  -  Index  of  computed  w  plana  transfer  function 
J  input  -  Index  of  z  plana  transfer  function  to  be  transformed. 

1.  The  sampling  period,  SAMPT,  of  COMMON/HEADDB/  must  bo  set  before  calling 
this  subroutine. 

Method 

Bilinear  transformation  is  implemented  by  transformation  of  the  z  plana 


All  of  the  following  parameters  are  initialized  by  subroutines  INITO  and 
MINITO  which  should  be  called  as  the  first  and  second  executable  statements  in 
the  main  program. 


Paran.  Preset 


Description 


CONTP  0 


CYCLE 

0 

DBMAX 

0 

DBMIN 

0 

DEGMN 

-360. 

DELAY 

0 

ECRE1 

2.E-4 

ECRE2 

l.E-8 

EDB1 

1. 

EDB2 

2. 

EDEG1 

4. 

EDEG2 

10. 

EPAD1 

l.E-10 

EP1 

l.E-0 

EP2 

l.E-10 

EP3 

.01 

EP4 

l.E-8 

EP5 

l.E-5 

ERCNJ 

l.E-4 

ERCX 

l.E-4 

ERCZ 

l.E-5 

FAUTO 

1 

Film  plot  continuation  flag  (cannot  be  used  for  Bode  plots) 

. EQ. 0  -  Single  plot 

.EQ.l  -  First  curve  of  a  plot 

.  EQ.2  -  Continuation  of  a  plot 

.  EQ.3  -  Final  curve  of  a  plot 

.EQ.O  for  automatic  selection  of  2  or  3  cycle  seal  e 
for  BODE  plots.  (1  cycle  is  not  available) 

Max.  db  for  plotting  frequency  response 
Min.  db  for  plotting  frequency  response 
(auto,  scaling  if  DBMAX. EQ. DBMIN) 

Min.  defined  phase  in  frequency  response 
(phase  defined  from  DEGMN  to  DEGMN+360.) 

Delay  time  for  sampled-data  transform  -  (sec.) 

Tolerance  for  common  root  in  common  root  elimination 
subroutine  CRELIM 

Tolerance  for  zero  root  in  subroutine  CRELIM 
Min.  delta  db  in  freq.  response  for  plotting 
Min.  delta  db  in  freq.  response  for  print  out 
Min.  delta  deg.  in  freq.  response  for  plotting 
Min.  delta  deg.  in  freq.  response  for  print  out 
Tolerance  for  negligible  higher  order  coefficients 
in  subroutine  ADDP 

An  input  quantity  for  subroutine  MULE 

An  input  quantity  for  subroutine  MULE 

An  input  quantity  for  subroutine  MULE 

An  input  quantity  for  subroutine  MULE 

An  input  quantity  for  subroutine  MULE 

Tolerance  to  determine  if  complex  root  pairs  are 

conjugates 

Roots  are  considered  to  be  complex  if  imag.  part 
.GT.  ERCX 

Roots  are  considered  to  be  zero  if  CABS.LT.ERCZ 
. NE.Q  for  automatic  frequency  mode.  Uses  NOMEG  and 
OMEGA  arrays. 

.EQ.O  for  user  supplied  frequency  points.  Uses 
FREQlf  FREQ2,  ...»  FREQ5  arrays. 

. NE.  0  for  Bode  plots  with  frequency  response 
Time  delay  for  s-plane  frequency  response 


FBODE 

FDLAY 


1 

0 


FILM 

0 

FNICO 

0 

FNYQS 

0 

FREQl(l) 

1. 

FREQK2) 

10. 

FREQ3(1) 

1. 

FREQk(l) 

0 

FREQk(2) 

0 

FREQk(3) 

0 

(k=2, 5) 

FXYDL 

.5 

FXYMN 

-2.5 

GRAFP 

1 

ITLOC 

50 

KDELT 

l.E+4 

KFLG 

1 

KGAIN 

.5 

2. 

KNORM 

1 

MAXIT 

80 

MDEG 

0 

MNDW 

.0005 

MTGER 

1 

MX  DM 

.22 

MXITF 

3000 

MXM 

1 

MO 

0 

Ml 

0 

M2 

0 

M3 

0 

M4 

0 

. NE.O  for  hardcopy  (high  resolution)  plots 
.NE.O  for  Nichols  plot  with  frequency  response 
.NE.O  for  Nyquist  plot  with  frequency  response 
Starting  frequency  point  for  first  segment  of  user 
specified  values. 

End  frequency  point  for  second  segment  of  user 
specified  values. 

Delta  frequency  for  third  segment  of  user 
specified  values. 

Starting  frequency  point  for  k-th  segment,  * 

End  frequency  point  for  k-th  segment,  * 

Delta  frequency  for  k-th  segment,  X 
(X  -  only  if  FAUTO.EQ.O) 

Nyquist  plot  scale  in  units  per  inch 
(Auto,  scaling  if  FXYDL=0) 

Nyquist  plot  parameter  -  minimum  real  and  imag.  value 
plotted 

.NE.O  for  printer  (low  resolution)  plots 

Max.  no.  of  different  gains  computed  in  root  locus 

Value  for  changing  root  locus  gains 

(preset  to  large  value  so  that  no  additional  gains 

are  computed  by  the  program) 

Flag  for  computing  root  locus  gains 
.EQ.O  for  computing  gains  by  ratio  KDELT 
.NE.O  for  computing  gains  by  increment  KDELT 
Array  of  values  used  for  computing  root  locus  gains 
KGAIN(l)=first  user  specified  root  locus  gain 
KGAIN(2)=second  user  specified  root  locus  gain 


KGAIN(NLOCI)=last  root  locus  gain 

(Gains  computed  and  used  only  if  they  are  between 

KGAIN(l)  and  KGAIN(NLOCI)  ) 

Gain  used  for  normalizing  polynomial  or  transfer  function 
Max.  no.  of  iterations  allowed  per  root  in  subroutine 
MULE 

Degree  of  highest  order  polynomial  in  matrix  (0-4) 

Min.  relative  freq.  step  size  in  freqy.  response  when 
FAUTO.NE.O 

Integer  m  for  multirate  configuration 

Max.  relative  freq.  step  size  in  freqy.  response  when 

FAUTO.NE.O 

Max.  no.  of  iterations  in  auto,  mode  of  freq.  response 
Dimension  of  matrices  (1-30) 

Matrix  for  coefficients  of  sXXQ 
Matrix  for  coefficients  of  sxxi 
Matrix  for  coefficients  of  sxx2 
Matrix  for  coefficients  of  sXX3 
Matrix  for  coefficients  of  sXX4 


HANOT  4  -  N unbar  of  linos  of  annotations  on  hardcopy  (high 

rasolution  aloctrostatic)  plot  (0-4) 

NLOCI  2  -  Nuabar  of  root  locus  gains  antarod  in  array  KGAIN 

(sax* 25) 

NOMEG  2  -  Nunbar  of  frequency  points  antarod  in  OMEGA  array 

for  usa  in  auto,  fraquancy  soda  (maxs20) 

NP  2  Flag  for  dataraining  output  froa  subroutina  MULE 

.E4.0  -  print  all  itarants  and  BCI  output  for  special 
procedures 

.EQ.l  -  print  only  tha  final  iteration  of  each  root 
.EQ.2  -  suppress  all  internal  printing 
(note  -  program  always  print  final  iteration  if  max 
iteration  obtained) 

NQDB  0  -  . NE.O  for  hardcopy  Nyquist  plot  in  db 

NRMFG  0  -  Polynoaial  And  Transfer  Function  Normalization  Flag 

If  .EQ.O,  tha  low  ordar  non-zaro  coefficient  of  tha 
danoainator  is  sat  equal  to  tha  value  of  KNORM  and 
all  other  coefficients  are  normalized  to  this  value. 

For  tha  a  planar  if  KN0RM=1.,  tha  low  order  non-zaro 
cooffieiant  of  tha  numerator  is  the  Bode  gain. 

If  . NE.O,  tha  high  ordar  coefficient  of  tha  denominator 
is  sat  equal  to  the  value  of  KNORM  and  all  other  coeffi' 
ciants  are  normalized  to  this  value.  For  the  s  plane, 
if  KN0RM=1 . ,  tha  high  ordar  coefficient  of  tha  numer¬ 
ator  is  tha  root  locus  gain. 

NTGER  1  -  Intagar  n  for  multi rata  configuration 

OMEGA  -  Array  of  fraquancy  points  for  auto,  frequency  moda 

(units  defined  by  RAD) 

1.  OMEGA(l)  *  first  fraquancy  point  used  in  auto,  mode 

10.  0MEGA(2)  »  second  fraquancy  point  used  in  auto,  mode 


0MEGA( NOMEG) • a  last  freq.  point  used  in  auto,  mode 


PMARG 

0 

. NE.  0  for  plotting  phasa  margin  instead  of  phase  for 
tha  Nichols  plot 

POLY 

0 

- 

Array  used  to  input  coefficients  for  polynomials 

P01YN 

0 

* 

Array  used  to  input  numerator  coefficients  for  transfer 
functions 

POL  YD 

0 

Array  used  to  input  denominator  coefficients  for 
transfer  functions 

PRNFLG1 

1 

- 

(not  usad) 

PRNFLG2 

1 

.EQ.O  to  suppress  print  out  of  arguments  of  an  LCAP2 
operation 

PRNFLG3 

1 

- 

.EQ.O  to  suppress  print  out  of  LCAP2  operation  executed 

PRNFLG4 

1 

* 

.EQ.O  to  suppress  print  out  of  results  of  an  LCAP2 
operation 

PRNFLG5 

1 

- 

(not  used) 

PRN1 

0 

- 

.EQ.O  for  suppressing  print  out  from  CRELIM 

PRN2 

0 

- 

.EQ.O  for  suppressing  print  out  from  FREQS 

PRN3  0  -  .£9.0  for  suppressing  print  out  fron  RCLAS 

PRN4  0  -  .EQ.O  for  suppressing  print  out  fron  RESDU 

PRN5  0  -  .EQ.O  for  suppressing  print  out  fron  STIME 

PRN6  0  -  .EQ.O  for  suppressing  print  out  fron  MTRANS 

PRN7  0  -  .EQ.O  for  suppressing  print  out  fron  PR00T  end  MR00T1 

PRN10  0  -  . NE.O  for  print  out  of  psrtiel  frection  coefficients 

fron  subroutine  RESDU 

RAD  1  -  .NE.O  For  frequency  in  red/sec 

.EQ.O  For  frequency  in  Hz 

RLF61  -1  -  Fleg  for  numbering  root  locus  points  on  the  herdcopy 

(higher  resolution)  plots 
=  1  for  numbering,  =  -1  for  no  numbering 
RLXMN  -450.  -  Root  locus  plot  peram.,  min.  x  scale 

RLXMX  50.  -  Root  locus  plot  perem.,  max.  x  scale 

(euto.  scaling  of  x  axis  if  RLXMN=RLYMX) 

RLYMN  -50.  -  Root  locus  plot  param.,  min.  y  scale 

RLYMX  450.  -  Root  locus  plot  param.,  max.  y  scala 

(auto,  scaling  of  y  axis  if  RLYMN=RLYMX) 

ROOT  0  -  Array  used  to  input  roots  of  polynomials 

ROOTN  0  -  Array  used  to  input  numerator  roots  for  transfer 

functions 

R00TD  0  -  Array  used  to  input  denominator  roots  for  transfer 

functions 

RTMAX  l.E+7  -  Max.  root  to  be  found  by  subroutine  MULE 

RZERO  l.E-5  -  Roots  .LT.  RZERO  returned  fron  subroutine  MULE  ere  set 

identicelly  to  zero 

SAMPT  1  -  Sampling  period  -  (sec.) 

SHADE  16  -  Pen  intensity  for  electrostetic  plot  (Z-Z8) 

TDELT  1.  -  Incremental  time  for;  s-plene  tine  response 

TEND  1.  -  End  tine  for  tine  response  celculetion 

TMAGN  1.  -  Megnitude  of  input  for  tine  response  celculetion 

TSTEP  1  -  .NE.O  for  stop  response  when  computing  time  response 

.EQ.O  for  impulse  response  when  computing  time  response 
TXMAX  0  -  Meximun  x  exis  for  time  response  plot 

TXMIN  0  -  Minimum  x  axis  for  time  response  plot 

(euto.  seeling  of  time  exis  if  TXMIN*TXMAX) 

TYMAX  0  -  Meximun  y  exis  for  time  response  plot 

TYMIN  0  -  Minimum  y  exis  for  tine  response  plot 

(euto.  seeling  of  y  exis  if  TYMINSTYMAX) 

TZEP1  1.  -  Inverse  z  transform  computational  error  criterie. 

Percent  difference  between  single  end  double  precision 
for  which  time  of  first  occurrence  will  be  printed  out 
TZERO  0  -  Stert  time  for  eveluating  s  plane  time  response 

TZFLG  0  -  .EQ.O  for  step  plot  in  inverse  z  transform  response 

plot  (hardcopy  only) 

.NE.O  for  symbol  plot  in  inverse  z  transform  response 
plot  (hardcopy  only) 

XGAP  50.  -  Plot  pen  lifts  up  when  x  distance  exceeds  XGAPX(.Ol) 

inches 

Y  position  for  second  line  of  herdcopy  plot  ennotetion 


YANOT 


9.6 


-  (in.)  (range  is  0-10) 

Plot  pen  lifts  up  when  y  distance  exceeds  YGAPK(.Ol) 
inches 

Pen  intensity  for  accenting  zero  line  for  hardcopy  plots 
. NE.0>  include  zero  order  hold  in  computation  of  sampled 
data  transform 


.EQ.0>  do  not  include  zero  order  hold 


SUBROUTINE  HEADINi  -  Heading  Statement  For  Entering  Plot  Titles 

(i=l,5) 

Purpose 

Simple  FORTRAN  statement  for  entering  Hollerith  data  into  the  plot  array 
HEAD  of  COMMON/HEADDB/  used  for  labeling  plot  titles. 

USflflfi 

CALL  HEADINi (INDX, WORD) 

INDX  input  -  Pointer  to  array  HEAD  of  COMMON/HEADDB/  where  array 
WORD  will  be  copied  into 

HORD  input  -  Hollerith  data  with  format  10H...  if  i=l 

20H. . .  if  i=2 

50H . . .  if  i=5 

1.  First  line  of  plot  title  is  in  HEAD(i),i=l,7 
Second  line  of  plot  title  is  in  HEAD(i ) , i=8 , 1A 
Third  line  of  plot  title  is  in  HEAD(i ) , i  =  15» 21 
Fourth  line  of  plot  title  is  in  HEAD(i ) , i  =  22, 28 

2.  Example i  CALL  HEADIN2(8,20HTHIS  IS  AN  EXAMPLE  )  will  yield, 

HEAD(3)=10HTHIS  IS  AN 
HEAD( 4 )  =  1 OHEXAMPL  E 

which  will  appear  as  the  second  line  of  the  plot  title. 

3.  First  line  of  plot  title  will  appear  at  the  top  of  the  plot. 

Second  line  will  begin  a  YANOT  units  from  the  bottom  of  page 
(full  scale  defined  from  0-10  units).  YANOT  (preset=9.6)  is 
in  COMMON/HEADDB/ . 


SUBROUTINE  REMARKi  -  Print  Out  Remarks,  (i=l,5) 

Purpose 

Single  FORTRAN  statement  for  printing  out  Hollerith  data. 

Usaflfl 

CALL  REMARKi(A) 

A  input  -  Hollerith  data  with  format  10H. . .  if  i=l 

20H...  if  is2 
30H. . .  if  i=3 
40H. . .  if  i=4 
50H. . .  if  i=5 

1.  Example i  CALL  REMARK2C20HTHIS  IS  AN  EXAMPLE  )  will  print  out. 


THIS  IS  AN  EXAMPLE 


SFAUX 


COMPLEX  FUNCTION  SFAUX  -  Evaluate  S  Plane  Transfer  Function  Coefficient  Array 
PUTBOSfl 

Evaluate  s  plane  transfer  function  coefficient  array  (LCAP2  format)  for  use 
in  computing  the  frequency  response.  This  complex  function  can  be  be  used  by 
subroutine  FREQS1  or  FREQS2  to  evaluate  the  transfer  function  specified  by  its 
first  argument.  It  can  also  be  used  by  user-supplied  subroutines  similar  to 
SFAUX1 . 

Usage 

SFAUX(TFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

SFAUX  output  -  Complex  value  of  response 

1.  Independent  s  plane  frequency  to  be  used  in  evaluating  the  response  is 
determined  in  subroutine  FREQS1  or  FREQS2. 


SFAUX1 


Identification 

COMPLEX  FUNCTION  SFAUX1  -  Evaluate  S  Plane  Transfer  Function  Coefficient  Array 
Purpose 

This  complex  function  is  similar  to  SFAUX  except  that  it  is  written  so  that 
it  can  be  easily  modified  by  the  user  to  allow  creation  of  a  user — defined  s  plane 
transfer  function. 

Usage 

SFAUX1 CTFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

SFAUX1  output  -  Complex  value  of  response 

Method 

This  complex  function  has  only  one  line  of  code 
SFAUX1=SFAUX(TFC) 

so  that  it  will  yield  the  same  results  as  SFAUX. 

To  create  a  user  defined  s  plane  transfer  function?  a  different  value  is 
returned  for  SFAUX1 .  For  example?  if  the  function  is 

SPTF2  +  SPTF4/2. 

the  user  would  change  the  FORTRAN  code  to 

SFAUX1=SFAUX(SPTF2)  +  SFAUXC SPTF4)/2 . 


The  argument  TFC  must  be  an  array  in  memory.  Since  only  the  first  five  s 
plane  transfer  functions  are  in  COMMON/SCMBLKB/  and  all  others  are  on  a  disk 
file?  only  SPTF1?  SPTF2?  SPTF3?  SPTF4  and  SPTF5  can  be  used  to  create  a 
user — defined  s  plane  transfer  function.  However?  the  user  can  define  additional 
transfer  function  coefficient  arrays  in  a  separate  labeled  common  block  to  be 
accessible  by  SFAUX.  Subroutine  FETSTF  can  be  used  to  copy  transfer  functions 
from  the  disk  file  to  the  transfer  function  in  this  common  block. 

The  code  for  this  routine  is  included  so  that  Example  16?  which  uses  this 
routine?  can  be  better  understood.  The  code  for  is  routine  is 


COMPLEX  FUNCTION  SFAUXl(TFC) 

C  THIS  FUNCTION  MILL  YIELD  RESUITS  WHICH  ARE  IDENTICAL  TO 
C  FUNCTION  SFAUX  IF  STATEMENT  100  BELOW  IS  NOT  CHANGED. 

C 

C  THE  INTENT  OF  THIS  FUNCTION  IS  TO  PROVIDE  THE  LCAP2  USER  WITH 

C  A  ROUTINE  WHICH  CAN  BE  EASILY  MODIFIED  TO  ALLOW  CREATION  OF  A  USER 

C  DEFINED  S-PLANE  FUNCTION  TO  BE  USED  WITH  SUBROUTINE  FREQS.  THE 
C  PROCEDURE  IS  TO  REPLACE  STATEMENT  100  WITH  A  USER  DEFINED 

C  FUNCTION. 

C 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


C  *  EXAMPLE  -  IF  THE  USER  DEFINED  FUNCTION  IS  TO  BE  * 
C  x  x 
C  X  SPTF2  +  SPTF4/2.  * 
C  *  x 
C  X  REPLACE  STATEMENT  100  BY  X 
C  X  x 
C  X  100  SFAUX1  =  SFAUX(SPTF2)  +  SFAUX(SPTF4)/2.  X 


C  XXXXXXXXKXXXXXXXXXXXXXXKXKXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C  XXXXXXKXXXXXXKKXXXXXXXXXXXK  NOTE  XXKXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

C  ONLY  THE  5  S-PLANE  TRANSFER  FUNCTIONS  SPTF1 , SPTF2, SPTF3, SPTF4, SFTF5 
C  CAN  BE  USED  DIRECTLY  WITH  FUNCTION  SFAUX.  ALL  OTHER  STORED  S-PLANE 
C  TRANSFER  FUNCTIONS  ARE  ON  DISK  STORAGE  AND  THEREFORE  NOT  DIRECTLY 
C  ACCESSIBLE  WITH  FUNCTION  SFAUX. 

C  IN  ORDER  TO  USE  ANY  OF  THESE  OTHER  S-PLANE  TRANSFER  FUNCTIONS, 

C  THE  USER  MUST  (1)  DECLARE  AN  SCM  COMMON  BLOCK  IN  THE  MAIN  PROGRAM 

C  AND  IN  SFAUX1 ,  AND  (2)  TRANSFER  THE  DESIRED  TRANSFER  FUNCTIONS 
C  INTO  THIS  COMMON  BLOCK.  SEE  E.A.  LEE  FOR  DETAILS. 

C 

COMMON/ FRQBLK/U,X,TWOPI,MMTGER 
COMPLEX  X 

COMMON/HEADDB/HEAD(70),DB(900) 

COMMON/SCMBLK/XTFSL 1520 ) , XTFW( 1520 ) , XTFZ( 1520) , XPYL760) 

DIMENSION  SPTF1(102),SPTF2(102),SPTF3(102),SPTF4(102),SPTF5(102) 
EQUIVALENCE  (XTFS( 1 ) , SPTF1U ) ) , (XTFS( 305) , SPTF2( 1 ) ) 

♦, CXTFS(609) ,SPTF3(1 ) ),<XTFS(913), SPTF4< 1 )),(XTFS(1217), SPTF5C 1 ) ) 
DIMENSION  WPTFK102) , WPTF2(  102) , WPTF3L 102) , WPTF4(  102) , WPTF5C 102) 
EQUI VAL ENCE  ( XTFWt 1 ) , WPTF1 ( 1 ) ) , < XTFWt  305 ) , WPTF2( 1 ) ) 

♦  , CXTFWC609) ,WPTF3( 1 ) ) , t XTFWt 913) ,WPTF4t 1 ) ) , <  XTFW( 1217 ) , WPTF5( 1 ) ) 
DIMENSION  ZPTF1( 102) ,ZPTF2( 102) , ZPTF3( 102 ) , ZPTF4L 1 02) ,ZPTF5( 102) 
EQUI  VAL  ENCE  (  XTFZU ) ,  ZPTF1  ( 1 )  ) ,  (  XTFZ  C  3  05  ) ,  ZPTF2  ( 1 ) ) 

♦, (XTFZ(609) ,ZPTF3( 1 ) ) , (XTFZL913) ,ZPTF4( 1 )),(XTFZ(1217 ) ,ZPTF5( 1 ) ) 
COMPLEX  SFAUX 
100  SFAUX1*SFAUX(TFC) 

RETURN 

END 


SFREQY 


If  LCAP2  defined  transfer  function  coefficient  arrays  are  to  be  used,  only 
the  first  five  transfer  functions  for  each  plane  are  available,  since  the  others 
are  on  disk  files.  However,  a  user  common  block  can  be  defined  so  that  these  oth¬ 
er  transfer  functions  can  be  first  transferred  from  disk  file  to  memory  with 
subroutine  FETSTF  so  that  SFREQY  can  be  used. 


COMPLEX  FUNCTION  WFAUX  -  Evaluate  W  Plane  Transfer  Function  Coefficient  Array 

PMrpggg 

Evaluate  w  plane  transfer  function  coefficient  array  (LCAP2  format)  for  use 
in  computing  the  frequency  response.  This  complex  function  can  be  used  by  sub¬ 
routine  FREQW1  or  FREQW2  to  evaluate  the  transfer  function  specified  by  its 
first  argument.  It  can  also  be  used  by  usei — supplied  subroutines  similar  to 
HFAUX1 . 

This  subroutine  can  also  evaluate  the  multirate  (fast  input*  slow  output) 
response  of  the  transfer  function. 

Usage 

WFAUX(TFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

MFAUX  output  -  Complex  value  of  response 

1.  Independent  w  plane  frequency  used  in  evaluation  of  the  response  is  com¬ 
puted  by  the  program  using  real  frequency  X  of  COMMON/FRQBLK/  and  sampling 
period  SAMPT  of  COMMON/HEADDB/ . 

2.  If  MMTGER  of  COMMON/ FRQBLK/  is  .GT.O,  the  multirate  response  is  computed  by 
using  Sklansky's  frequency  decomposition  method.  MMTGER  is  the  ratio  of  the 
(output/input)  sampling  periods  and  SAMPT  is  the  sampling  period  of  the 
faster  input  sampler. 


WFAUX1 


Identification 

COMPLEX  FUNCTION  HFAUX1  -  Evaluate  W  Plane  Transfer  Function  Coefficient  Arrav 
£MT-PgS8 

This  complex  function  is  similar  to  NFAUX  except  that  it  is  written  so  that 
it  can  be  easily  modified  by  the  user  to  allow  creation  of  a  user — defined  w  plane 
transfer  function. 

Usage 

HFAUX1CTFC) 

TFC  input  -  Transfer  function  coefficient  array  CLCAP2  format) 

HFAUX1  output  -  Complex  value  of  the  response 

Method 

This  complex  function  has  only  one  line  of  code> 

WFAUX1=WFAUX(TFC) 

so  that  it  will  yield  the  same  results  as  WFAUX. 

To  create  a  user  defined  w  plane  transfer  function,  a  different  value  is 
returned  for  MFAUX1 .  For  example,  if  the  function  is> 

WPTF2  +  HPTF4/2. 

the  user  would  change  the  FORTRAN  code  to 

HFAUX1 = WFAUX ( WPTF2 )  +  WFAUXCWPTF4V2. 


HFREQY 


Identifi cation 

SUBROUTINE  HFREQY  -  Evaluate  Frequency  Response  Of  A  W  Plane  Transfer 

Function  Coefficient  Array 


Purpose 

Evaluate  frequency  response  of  a  w  plane  transfer  function  coefficient 
array.  User  supplies  name  of  the  array. 

Usage 


CALL  HFREQY(TFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

1.  Frequency  response  parameters  are  in  COMMON/HEADDB/ .  See  description  of 
subroutine  SFREQ  for  definition. 

B«gtriQUfl,ng 

If  LCAP2  defined  transfer  function  coefficient  arrays  are  to  be  used,  only 
the  first  five  transfer  functions  for  each  plane  are  available,  since  the  others 
are  on  disk  files.  However,  a  user  common  block  can  be  defined  so  that  these  oth¬ 
er  transfer  functions  can  be  first  transferred  from  disk  file  to  memory  so  that 
HFREQY  can  be  used. 


COMPLEX  FUNCTION  ZFAUX  -  Evaluate  Z  Plane  Transfer  Function  Coefficient  Array 


Purpose 

Evaluate  z  plana  transfer  function  coefficient  array  (LCAPZ  format)  for  use 
in  computing  the  frequency  response.  This  complex  function  can  be  used  by  sub¬ 
routine  FREQZ1  or  FREQZ2  to  evaluate  the  transfer  function  specified  by  its 
first  argument.  It  can  also  be  used  by  user-supplied  subroutines  similar  to 
ZFAUX1 . 

This  subroutine  can  also  evaluate  the  multirate  (fast  input*  slow  output) 
response  of  the  transfer  function. 

Uaaju 


ZFAUX(TFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

ZFAUX  output  -  Complex  value  of  response 

1.  Independent  z  plane  frequency  used  in  evaluation  of  the  response  is  com¬ 
puted  by  the  program  using  real  frequency  X  of  COMMON/FRQBLK/  and  sampling 
period  SAMPT  of  COMMON/HEADDB/ . 

2.  If  MMTGER  of  COMMON/FRQBLK/  is  .GT.O*  the  multirate  response  is  computed  by 
using  Sklansky's  frequency  decomposition  method.  MMTGER  is  the  ratio  of  the 
(output/input)  sampling  periods  and  SAMPT  is  the  sampling  period  of  the 
faster  input  sampler. 


COMPLEX  FUNCTION  ZFAUX1  -  Evaluate  Z  Plane  Transfer  Function  Coefficient  Array 
Purpose 

This  complex  function  is  similar  to  ZFAUX  except  that  it  is  written  so  that 
it  can  be  easily  modified  by  the  user  to  allow  creation  of  a  user — defined  z  plane 
transfer  function. 

Usage 

ZFAUXl(TFC) 

TFC  input  -  Transfer  function  coefficient  array  (LCAP2  format) 

ZFAUX1  output  -  Complex  value  of  the  response 

Method 

This  complex  function  has  only  one  line  of  code 
ZFAUX1=ZFAUX(TFC) 

so  that  it  will  yield  the  same  results  as  ZFAUX. 

To  create  a  user  defined  z  plane  transfer  function,  a  different  value  is 
returned  for  ZFAUX1 .  For  example,  if  the  function  is 

ZPTF2  +  ZPTF4/2. 

the  user  would  change  the  FORTRAN  code  to 

ZFAUX! =ZFAUX(ZPTF2)  +  ZFAUX(ZPTF4)/2 . 


Alphanumeric  information  for  labeling  plots  are  in  array  HEAD  of 
COMMON/HEADDB/ .  For  hardcopy  (high  resolution  electostatic)  plots,  up  to  four 
lines  of  annotation  are  available.  For  printer  plots,  only  one  line  of  anno¬ 
tation  is  available.  The  data  in  this  array  is  used  as  follows: 

HEAD(i),  i=l,7  for  1st  line,  hardcopy  and  printer  plots 

HEAD(i),  i=8,14  for  2nd  line,  hardcopy  plots  only 

HEAD(i),  i=15,21  for  3rd  line,  hardcopy  plots  only 

HEAD(i),  i=22,28  for  4th  line,  hardcopy  plots  only 

This  HEAD  array  is  preset  to  blanks.  Whenever  a  hardcopy  plot  is  made  the  con¬ 
tents  in  the  HEAD  array  will  be  printed  out.  The  user  must  enter  or  change  the 
contents  in  this  array  prior  to  a  FORTRAN  call  (an  LCAP2  operator)  which  produc¬ 
es  a  plot. 

For  example,  if  the  first  line  is  to  be  "EXAMPLE  1  S  PLANE  FREQUENCY 
RESPONSE,"  either  one  of  the  following  can  be  used: 

Method  1 

HEAD(1 )=10HEXAMPLE  1 
HEAD(2)=10HS  PLANE  FR 
HEAD(3)=10HEQUENCY  RE 
HEAD(4)=10HSPONSE 

Mtihqd-2 

CALL  HEADIN4( 1 , 40HEXAMPL E  1  S  PLANE  FREQUENCY  RESPONSE  ) 

Method  i 

ENCODE! 40, 100, HEAD( 1 ) ) 

100  FORMAT ( 40HEXAMPL E  1  S  PLANE  FREQUENCY  RESPONSE  ) 

Method  2  is  the  simplest  to  use  since  only  a  single  line  of  FORTRAN  code  is 
required.  This  subroutine  is  of  the  form  HEADINi, (i-1,5)  where  i  designates  the 
multiples  of  10  characters  to  be  entered  and  the  first  argument  of  the  subrou¬ 
tine  designates  the  starting  location  of  the  HEAD  array  where  the  alphanumeric 
information  is  to  be  stored.  For  more  details,  see  description  of  this  subrou¬ 
tine  in  Appendix  C. 

Method  3  has  the  advantage  of  allowing  the  user  to  annotate  the  plots  with 
data  if  the  ENCODE  statement  were  to  include  a  list.  For  example,  if  the  second 


line  of  annotation  on  the  hardcopy  plot  were  to  be  labeled  with  the  values  of  the 
variables  FLEXW  and  ZETA,  the  following  FORTRAN  statements  could  be  usedt 


ENCODEC 31 , 200 , HEADC  8 ) ) FL  EXW, ZETA 
200  F0RMATU1HFL EXMODE  =,F7.3,8H,  ZETA  =,F5.3) 

If  the  values  of  FLEXW  and  ZETA  were  14.523  and  .707,  respectively,  HEADC8) 
through  HEADC 11)  would  have  the  following  characters  in  it. 


FLEX  MODE 

=  14.523, 

7 

HEADC8) 

HEAD(9) 

HEADC 10) 

HEADC 11 ) 

The  first  line  of  annotation  will  always  be  at  the  top  of  the  plot.  The  sec¬ 
ond,  third  and  fourth  lines  (hardcopy  plots  only)  normally  begin  just  below  the 
first  line.  These  last  three  lines,  however,  can  be  placed  lower  on  the  plot  if 
they  should  interfere  with  the  plot  data.  The  parameter  YANOT  (preset  to  9.6)  of 
COMMON/HEADDB/,  which  positions  the  y  coordinate  for  the  second  line,  can  be 
changed  to  a  smaller  value.  The  range  for  this  parameter  is  0-10,  which  corre¬ 
sponds  to  the  length  of  the  y  axis. 

For  the  hardcopy  plots  the  number  of  lines  of  annotation  to  be  produced  is 
determined  by  the  parameter  AN0TAT  of  COMMON/HEADDB/,  which  is  preset  to  4. 


The  source  code  for  this  program  is  available  to  agencies  supporting  DOD 
projects  and  studies.  The  requester,  however,  should  be  aware  that  some  non-ANSI 
FORTRAN  code  and  one  assembly  language  subroutine  are  utilized.  If  the  program 
is  be  be  run  on  a  CDC  176  or  7600  computer  under  the  SCOPE  2.1  operating  system, 
no  problems  should  be  encountered.  If  the  program  is  to  be  run  on  any  other  com¬ 
puter,  modifications  to  the  program  most  likely  will  have  to  be  made.  The 
following  facts  will  be  of  interest  if  modifications  are  to  be  made: 

(a)  A  FORTRAN  version  of  CXMTX1,  which  is  written  in  assembly  language,  is 
available. 

(b)  Non-ANSI  COC  FORTRAN  4  ENCODE  and  DECODE  statements  must  be  replaced  with 
ANSI  standard  internal  write  and  read  statements.  Even  if  the  target  com¬ 
puter  and  operating  system  supports  the  ENCODE  and  DECODE  statements  mod¬ 
ifications  might  still  be  necessary  since  the  ENCODE  statements  used  in 
LCAP2  makes  use  of  the  fact  that  the  CDC  word  length  is  60  bits  long. 

(c)  Work  has  been  initiated  to  convert  this  program  to  FORTRAN  5  and  have  it 
operational  for  the  IBM  3033  computer  as  well. 

This  program  as  well  as  this  user's  guide  is  in  a  continuous  process  of 
evolution  and  development.  For  these  reasons,  this  program  and  related  materials 
will  be  made  available  under  the  understanding  that  no  warranty,  express  or 
implied,  is  made  by  the  Aerospace  Corporation  as  to  the  accuracy  and  functioning 
of  the  program  and  related  materials  and  that  no  responsibility  for  program 
maintenance  is  implied. 

The  current  reproduction  and  handling  fee  is  $220.00.  Request  for  a  copy  of 
this  program,  which  also  includes  the  Interactive  version  of  LCAP2  as  well, 
should  be  addressed  to: 

Administrator 

Information  Processing  Division 

The  Aerospace  Corporation 

2350  E.  El  Segundo  Blvd. 

El  Segundo,  California  90245 
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